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PREFACE 


This  report  Is  the  ninth  issue  of  the  Air  Force  Academy  Aeronautics 
Digest  *  Our  policy  is  to  print  articles  which  represent  recent 
scholarly  work  by  students  and  faculty  of  the  Departnent  of  Aeronautics, 
members  of  other  departments  of  the  Academy  and  the  Frank  J.  Seiler 
Research  Laboratory,  researchers  directly  or  indirectly  involved  with 
USAFA-sponsored  projects,  and  authors  in  fields  of  Interest  to  the  USAFA. 

In  addition  to  complete  papers,  the  Digest  Includes,  when 
appropriate,  abstracts  of  lengthier  reports  and  articles  published  In 
other  formats.  The  editors  will  consider  for  publication  contributions  in 
the  general  field  of  Aeronautics,  including: 


Aeronautical  Engineering 
Aerodynamic  s 
Flight  Mechanics 
Propul  si  on 
S  true  tur e s 
Instrumentation 
Fluid  Dynamics 

Thermodynamics  and  Heat  Transfer 
Biomechanics 
Engineering  Education 
Aeronautical  History 


Papers  on  other  topics  will  be  considered  on  an  individual  basis. 
Contributions  should  be  sent  to: 


Editor,  Aeronautics  Digest 
DFAN 

US  Air  Force  Academy,  Cf  80840 


The  Aeronautics  Digest  is  presently  edited  by  MaJ  Jay  DeJongh,  PhD; 
Maj  A.M.  Higgins,  PhD;  Maj  E.J.  Jumper,  PhD;  and  Capt  J.M.  Kempf,  PhD, 
Department  of  English,  who  provided  the  final  editorial  review.  Our 
thanks  also  to  Associate  Editor,  Martha  Arends,  and  Production  Artist, 
Deborah  Ross,  of  Contract  Techical  Services,  Inc. 


*The  first  eight  Issues  of  the  Digest  can  be  ordered  from  the  Defense 
Technical  Information  Center  (DTIC),  Cameron  Station,  Alexandria,  VA 
22324.  Use  the  following  AD  numbers:  Aeronautics  Digest  -  Spring  1978, 
ADA060207;  Aeronautics  Digest  -  Fall  1978,  ADA069044;  Aeronautics 
Digest  -  Spring  1979,  ADA075419;  Aeronautics  Digest  -  Fall  197$, 
ADA085770;  Aeronautics  Digest  -  Spring/Summer  1980,  ADA093i78; 
Aeronautics  Digest  -  Fall/Wlnter  I960,  ADA108338;' Aeronautics  Digest  - 
Spr  lng/Summer  1981,  ADA1  12421  ;  Aeronautics  Digest  ~  Fall/Winter  19181. 
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THE  SUPERSONIC  AREA  RULE  AND  A  PROPOSED  SIMPLIFICATION 
Eric  J.  Jumper* 

Abstrac t 

This  paper  outlines  the  development  of  the  supersonic  area  rule  from 
theory  to  application.  Then  a  simplification  to  the  rule  is  Introduced 
which  greatly  reduces  the  complexity  Involved  in  using  the  full  supersonic 
area  rule  to  predict  wave  drag  increments.  Finally,  the  results  of 
applying  the  modified  rule  to  an  F-15  aircraft  configured  with  near-axis 
conformal  pallets  are  described. 

I .  Introduction 

In  supersonic  flight  drag  is  due  to  three  separate  phenomena.  The 
first  two  are  the  same  as  in  subsonic  flight:  skin-friction  drag  (l.e., 
viscous  drag)  and  pressure  drag  due  to  separated  flow.  In  supersonic 
flight  the  third  phenomenon  is  known  as  wave  drag,  which  is  due  to  the 
difference  in  pressure  between  regions  of  the  flow  separated  by  shock  or 
expansion  waves.  These  waves  emanate  from  flow  boundaries  (for  example, 
body  surfaces)  and  may  Intersect  with  waves  from  other  boundaries  causing 
interference  effects.  These  interference  effects  further  complicate  an 
already  complicated  analysis  problem. 

In  1956  Jones  (Ref.  1)  published  a  method  of  predicting  wave  drag 
that  included  the  Interference  effectB.  This  method  was  similar  in  form 
to  the  transonic  area  rule  and  was,  therefore,  christened  the  supersonic 
area  rule.  Since  the  appearance  of  Jones'  report,  aerodynamlclsts  have 
exploited  the  supersonic  area  rule  to  determine  wave-drag  effects  for 
aircraft  of  various  geometries  and  to  determine  the  wave-drag  effects  of 
changes  in  geometry  such  as  fuselage  shape  and  wing  location.  In  1957 
Nelson  and  Velah  (Ref.  2)  used  a  Fourler-serles  method  for  applying  the 
supersonic  area  rule  in  an  extensive  comparison  between  predicted  results 
obtained  by  using  the  rule  and  experimental  data  obtained  for  the  same 
configurations  lV  wind-tunnel  tests.  They  found  that  the  rule  provided 

•Major,  USAF,  Associate  Professor,  Department  of  Aeronautics  and 
Astronautics,  Air  Force  Institute  of  Technology,  Wright  Patterson  AFB, 
Ohio 
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predictions  of  wave  drag  which  were  reasonably  close  to  those  obtained  In 
the  wind  tunnel.  Then  In  1963  Harris  (Ref.  3)  published  a  report 
describing  a  computer  program  written  by  Boeing  Company  that  made  possible 
the  general  application  of  the  supersonic  area  rule  for  predicting  wave 
drag.  With  the  availability  of  high-speed  computers,  the  Boeing  program, 
or  derivatives  of  It,  came  into  standard  use  In  larger  and  more  extensive 
aerodynamic  prediction  and  design  programs  such  as  that  by  Baals,  Robins, 
and  Harris  (Ref.  4). 

Aerodynamic  prediction  programs  are  useful  not  only  in  the  design 
phase  of  planning  a  new  airplane  but  also  In  understanding  the  effect  of 
changes  to  existing  aircraft  configurations,  such  as  adding  wing  stores, 
or  enlarging  a  turret.  While  programs  like  the  one  written  by  Boeing  are 
generally  available  and  are.  In  fact,  used  in  conjunction  with  wind-tunnel 
tests  In  the  final  stages  of  the  design  process,  they  are  seldom  used  for 
early-management  decisions  or  In  system  design  studies.  The  main  reason 
for  this  Is  that  these  computer  programs  require  large  computers  and 
complicated  Input  parameters.  Further,  the  configurations  to  be 
considered  might  be  quite  numerous.  These  factors  add  up  to  an  extremely 
expensive  exercise. 

If  the  supersonic  area  rule  could  be  adapted  for  use  on  a  small 
computer.  It  would  greatly  benefit  early-management  decisions  and  system 
design  studies.  However,  the  supersonic  area  rule,  In  Its  present  form. 

Is  not  suited  for  use  on  smaller  computers  because  of  large  core  storage 
requirements  and  the  long  computation  time  necessary.  Looking  at  some 
specific  applications  of  the  area  rule  might  indicate  ways  In  which  the 
rule  could  be  streamlined  to  make  It  compatible  with  small  computers. 

One  example  of  a  program-management  decision  that  could  be  aided  by 
using  the  supersonic  area  rule  Is  the  placement  of  a  near-fuselage-axis 
store  or  equipment  pod.  This  decision  might  Involve  numerous  candidate 
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locations  chosen  for  other  than  aerodynamic  reasons,  such  as  structural 
considerations.  The  use  of  the  area  rule  would  limit  the  number  of 
configurations  to  be  studied  further  in  wind-tunnel  tests  and/or  in  more 
extensive  analytic  studies.  Because  this  application  involves  changing  an 
existing  configuration,  much  is  already  known  about  the  original  airplane 
and  only  drag  increments,  due  to  locating  the  pod  or  to  minor  changes  in 
the  pod  shape  at  the  same  location,  are  of  interest  to  the  analyst 

The  supersonic  area  rule  would  also  be  helpful  in  system  design 
studies;  however,  these  studies  often  determine  volume  requirements  rather 
than  specific  configurations.  For  example,  the  study  might  examine  the 
aerodynamic  performance  trade-off  of  adding  some  new  avionics  package 
which  can  be  accommodated  only  by  adding  volume  to  the  fuselage  of  an 
existing  aircraft.  If  the  envelope  of  the  aircraft  Includes  supersonic 
flight,  the  supersonic  area  rule  would  clearly  be  helpful.  However,  these 
system  studies  include  many  considerations  other  than  predicting  the 
wave-drag  effect,  making  the  large  computer-time  requirement  of  the 
supersonic  area  rule  prohibitive.  Further,  the  supersonic  area  rule 
requires  detailed  information  about  a  specific  configuration  and,  as 
mentioned  above,  the  only  configuration  information  which  is  likely  to  be 
available  is  the  volume  requirement. 

In  both  the  management-decision  example  and  the  system-design-study 
application,  the  requirement  for  precision  in  the  wave-drag  prediction  is 
relaxed,  and  the  theoretician  is  given  considerable  license  in  modifying 
the  supersonic  area  rule  to  streamline  its  use.  From  a  user's  point  of 
view,  the  most  important  modifications  to  the  supersonic  area  rule  would 
be  those  which  would  simplify  the  formulation  of  the  input  data  describing 
the  aircraft  configuration;  however,  modifications  which  reduce  process 
time  and  core  storage  requirements  of  the  associated  computer  program  are 
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also  important  and,  In  the  end,  may  be  the  deciding  factor  in  whether  the 
method  is  used. 

This  paper  describes  the  supersonic  area  rule  and  a  modification  to 
it  which  greatly  simplifies  its  use  and  computer  requirements.  A  program 
that  incorporates  the  modified  area  rule  is  also  described.  The  results 
of  using  the  supersonic  area  rule  and  the  modified  rule  are  compared  for  a 
number  of  aircraft  configurations  including  an  F-15  for  which  wind-tunnel 
data  had  been  collected  with  and  without  a  near-fuselage-axis  protuberance 
in  place. 


II .  Supersonic  Area  Rule 
A.  Theory 

Suppose  we  consider  an  airplane  in  a  steady  supersonic  flow  whoae 
surfaces  are  Inclined  to  the  freestream  by  only  a  small  amount.  The 
ideal-fluid  flow  field  around  the  airplane  is  then  amply  described  by  the 
linearised  potential  equation 


6  4> 


yy 


(i) 


where  B2  “  M2 -l ,  is  the  freestream  Mach  number,  ♦  is  the  velocity 

potential,  and  the  subscripts  refer  to  differentiation  in  the  streamvlse 

direction,  x,  and  the  two  other  orthonormal  coordinates,  y  and  a,  of 

cartesian  three-space.  Lomax  (Kef.  5)  showed  that  the  general  solution  of 

Eqn .  (1)  for  the  wave-drag  coefficient,  ,  for  a  lifting  or  non-lifting 

w 

airplane  is  given  by 


2w 
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where 


In  Eqns.  (2)  and  (3)  6  is  measured  from  a  fixed  coordinate  relative  to  the 
freestream  and  indicates  the  orientation  of  a  plane  tangent  to  a  Mach  cone 
for  M  whose  angle  of  inclination,  p  ,  is  sin  ^(l/M-,,).  S,  is  some 

<X>  D 

reference  area  and  q  is  the  reference  pressure,  one-half  the  density  times 
the  freestream  velocity  squared  (®5pV^).  To  define  the  area  S(x),  let  the 
x-axis  of  a  cartesian  coordinate  system  be  fixed  parallel  to  the 
freestream  with  its  origin  at  the  nose  (or  most  forward  point)  of  the 
aircraft,  and  let  a  Mach  plane  of  orientation  6  intersect  the  x-axis  at  a 
point  x.  Then  S(x)  is  defined  as  the  projection,  normal  to  the  x-axis,  of 
the  area  formed  by  the  intersection  of  the  Mach  plane  and  the  aircraft. 
Figure  la  shows  S(x)  at  two  orientation  angles,  tt  / 2  and  an  arbitrary 
angle,  9  .  In  Eqn.  (3),  f  is  the  net  force  conormal  to  the  stream 
direction  and  the  9  direction  (see  Figure  lb),  computed  by  a  closed 
integration  of  the  pressure  around  the  surface,  o  ,  cut  by  the  Mach  plane. 

One  might  think  of  the  coefficient  of  drag  given  in  Eqn.  (2)  as 
consisting  of  the  two  parts  given  in  Eqn.  (4): 


2  IT 


2n 

(9)de  +  J  cD 


(e)de  | 

J 


(4) 


where  the  (g)  consists  of  the  first  term  on  the  right  side  of  Eqn.  (3) 

o 

and  C  (9)  consists  of  the  last  two  terms  of  Eqn.  (3).  It  is  clear  that 


6 


USAFA-TR-83-2 


only  Cn  (9)  contains  terms  requiring  knowledge  o f  the  pressure 
°L 

distribution  on  the  surface  of  the  airplane.  Although  the  subscript  L  on 
the  term  is  slightly  misleading,  Lomax  (Ref.  5)  points  out  that  when  the 
orientation  of  e  Is  such  that  the  positive  net  conormal  force  f  la 
pointing  up  at  90  degrees  to  the  freestream,  f  represents  the  net  section 

lift.  Thus,  the  terms  of  c  (0)  do  in  a  way  resemble  drag  due  to  a  "lift” 

DL 

of  sorts.* 

The  supersonic  area  rule  as  described  by  Jones  (Ref.  1)  gives  the 
formula  for  the  approximate  coefficient  of  drag  as 


2n 


4w2qS 


-/ 

b  o 


r  l  l 

a  /* 

j  J  s"(0  S"  (x)  In 

0  0 

x  5 

L  L 

d(;dx 

dO 


(5) 


where  S(x)  Is  defined  In  a  similar  way  as  for  Eqn.  (3).**  Comparing  Eqn. 

(5)  to  Eqn.  (4)  shows  that  the  supersonic  area  rule  Is  simply  the  Cp  term 

o 

of  Lomax's  more  general  expression  which,  you  recall,  was  subject  only  to 
the  limitations  of  linearised  theory.  As  Nelson  and  Welsh  (Ref.  2) 
demonstrated  and  Lomax  (Ref.  S)  himself  pointed  out,  the  term  of  Eqn. 

(4)  represents  the  limiting  factor  to  the  correctness  of  the  supersonic 
area  rule,  Eqn.  (S),  even  within  the  framework  of  linearised  theory. 

One  cannot  help  but  notice  the  striking  resemblance  of  Eqn.  (5)  to 
von  Karmen' a  slender-body  result  for  a  closed  (or  constant  radius  at  x»L) 

body  of  revolution  (Ref.  6)  that  gives  the  wave-drag  coefficient  as 

L  L  .. 

S  '  sk  // ‘"(e>  Hr  '  £|  JM’‘  <‘> 

b  0  0 

*In  general,  even  nonllftlng  configurations  will  have  non-sero 
Cdt(  9)  '  • • 

**Jones  suggests  using  the  areas  cut  by  the  Mach  planes;  however.  Nelson 
and  Welsh  (Ref.  2)  found  that  using  the  same  definition  for  S(x)  as  In 
Eqn.  (3)  gave  better  results. 
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where  A(x)  Is  the  cross-sectional  area  of  the  body  of  revolution.  Ue  may 
write  the  integrand  of  Eqn.  (S)  as 

L  L 

CD  (9)  *  2^S 7  ff SM(e)  S"(X)  ln  ll  '  l|  d5dx  (7) 

W  b  0  0 

A  comparison  of  Eqn.  (7)  to  Eqn.  (6)  suggests  an  Interpretation  of  Eqn. 

(7)  as  the  von  Karman  coefficient  of  drag  for  an  "equivalent"  body  of 

revolution  for  a  given  family  of  Mach  planes  of  e  orientation.  In  this 

Interpretation  S  of  the  e  Mach  planes  appears  to  simply  replace  the  A  of 

Eqn.  (6).  It  is  in  this  light  that  Jones  (Ref.  1)  presented  the 

supersonic  area  rule.  It  Is  important  to  remember,  however,  that  it  Is  CQ 

of  Eqn.  (4)  which  forms  the  limitation  of  the  area  rule  and  not  the 

restriction  to  von  Karman's  Integral  formula,  which  requires  that  the 

body  be  closed  (or  constant  radius  at  x-L),  for  example. 

If  we  were  to  examine  a  closed  body  of  revolution  at  zero  angle  of 

attack,  it  is  Interesting  that  the  slender-body  theory,  Eqn.  (6),  leads  to 

a  coefficient  of  drag  which  is  independent  of  Mach  number.  In  this  case, 

CD  (0)  of  Eqn.  (7)  1 8  independent  of  e  so  that  equals  (0). 

O  o  o 

Additionally,  as  Lomax  pointed  out  (Ref.  5),  this  case  would  lead  to  the 

symmetry  properties  that 

S"(x,6)  -  S"  (x ,  -0)  (8) 

and 

df  (6)  _  _  df  (-9)  (9) 

dx  dx 

so  that  the  first  term  of  c  would  become  Identically  zero.  Further, 

DL 

since  there  is  no  net  str eamwl se-nor mal  force,  the  second  term  must  also 
be  zero.  For  the  body-of-revolut ion  problem,  then.  It  would  appear  that 
the  supersonic  area  rule  takes  Into  account  wave-drag  changes  with  Mach 
number  due  to  thickness  distribution  that  are  not  accounted  for  ln 
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B.  Application 

In  order  to  obtain  a  wave-drag  computation  within  the  framework  of 

linearized  theory,  one  would  need  a  complete  description  of  the  pressure 

distribution  on  the  surface  of  a  given  airplane;  having  that  pressure 

distribution,  however,  would  lead  directly  to  the  drag  by  means  of  a 

surface  integration.  With  the  exception  of  some  simple  configurations, 

pressure  information  this  complete  requires  extensive  computation,  perhaps 

based  on  paneling  techniques,  and  having  that  information  obviates  the 

need  for  employing  Eqn.  (3).  Lacking  this  information  the  supersonic  area 

rule,  Eqn.  (5),  has  been  employed  to  obtain  an  approximation  to  •  As 

w 

Lomax  (Ref.  1)  and  Kelson  and  Welsh  (Ref.  2)  pointed  out,  there  are 

aircraft  configurations  for  which  Cn  (the  term  which  limits  the  accuracy 

L 

of  the  supersonic  area  rule)  is  not  small.  But  for  supersonic  aircraft  of 
conventional  form  Baals,  Robins,  and  Harris  (Ref.  4)  have  shown  that  the 
area  rule  gives  good  results. 

In  practice,  numerical  computer  algorithms  are  used  to  carry  out  the 
integrations  in  Eqn.  (5).  One  of  the  first  widely  known  algorithms  was 
that  written  by  Boeing  Company  in  the  1960s.  I  believe  that  it  is  fair  to 
say  that  by  today's  standards  the  Boeing  program  seems  crude. 

Nonetheless,  the  Boeing  program  served  as  the  point  of  departure  for  many 
of  today’s  algorithms.  As  described  by  Harris  (Ref.  3)  the  Boeing  method 
calculates  CD(9)  for  each  of  a  finite  number  of  9’s  by  numerically 
integrating  Eqn.  (7)  using  the  equivalent  body  of  revolution  determined  by 
the  projection  on  an  axis-normal  plane  of  the  section  area  cut  by  the  Mach 
plane  oriented  at  the  roll  angle  9  .  Sufficient  6's  are  chosen  &6  apart 
so  that  9  goes  through  360  degrees.  is  then  computed  by  averaging  all 

cD(e)'s. 

A  major  problem  in  building  such  a  computational  algorithm  is  in 
describing  the  geometry  of  the  airplane.  Computer-interactive-graphics 
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methods  have  made  this  problem  a  alapler  one  than  It  once  via  from  the 
operator'e  point  of  view,  but  thia  amplification  waa  paid  for  at  the 
price  of  both  algorithm  and  coaputer  hardware  aophiatlcatlon  (Ref.  4). 

Ill .  Modification  of  the  Superaonic  Area  Rule 

It  night  at  firat  seem  ludicroua  to  auggeat  a  modification  to  the 
area  rule  that  is  less  aophiatlcated  than  the  published  works  of  the 
1960a.  The  modification  proposed  here,  however,  is  for  the  user  who  la 
without  access  to  large  computers  and  operates  with  little  or  no 
interactive  graphics  capability.  If  we  add  to  these  constraints  the 
demand  for  a  simple  input  capability,  it  is  clear  that  some  modification 
to  the  supersonic  area  rule  is  needed  if  it  is  to  be  used. 

The  modification  I  adopted  was  to  construct  a  simple  equivalent  body 
of  revolution  based  on  planes  cut  normal  to  the  fuselage  axis.  It  is  this 
equivalent  body  of  revolution  to  which  the  area  rule  is  then  applied. 

This  means  that  Eqn.  (7)  need  only  be  applied  and  numerically  Integrated 
once  for  each  Mach  number  of  Interest.  Further,  the  algorithm  for 
constructing  new  equivalent  bodies  of  revolution  at  each  Mach  number 
becomes  almost  trivial  when  applied  to  a  body  of  revolution.  In  addition, 
the  input  requirement  is  satisfied  by  a  simple  area  distribution  which,  if 
not  readily  available,  is  simple  to  construct.  And  finally,  any 
protuberance  to  be  added  to  an  existing  aircraft  may  be  specified  as  a 
simple  additive  area  at  applicable  axial  fuselage  locations,  thus  making 
it  ideal  for  systems  studies  where  these  areas  may  be  the  extent  of  the 
information  available. 

A.  Previous  Work 

This  modification  is  not  without  precedent.  Harris  (Ref.  3) 
described  a  similar  technique  for  simplifying  the  description  of  the 
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fuselage  In  the  early  Boeing  program;  however,  authors,  when  speaking  of 
the  supersonic  area  rule,  have  not  been  so  bold  as  to  collapse  the  entire 
wing  structure,  tall  structure,  etc.  onto  a  single  body  of  revolution.  N 
doubt  Investigators  along  the  way  have  done  so,  but  I  found  no  mention  of 
this  In  the  literature. 

Hall  (Ref.  7)  performed  a  beautiful  piece  of  experimental  work 
dealing  with  the  transonic  area  rule.  In  this  work  he  compared  drag  data 
taken  from  actual  aircraft  configurations  to  drag  data  taken  from 
equivalent  body-of-r ev olut ion  models  for  Mach  numbers  up  to  1.3.  The 
results  are  quite  marvelous.  Hall's  conclusion  was  that  the  simple 
equivalent  bodies  of  revolution  were  able  to  predict  drag  rise  of  the 
airplane  configurations  to  within  15  percent.  Selected  results  of  Hall's 
study  are  given  in  Figures  2,  3,  and  4. 

Since  the  proposed  modification  Co  the  supersonic  area  rule  is  the 
exact  (within  the  framework  of  linearized  theory)  computation  for  a  body 
of  revolution,  Hall's  work  tends  to  justify  the  use  of  the  simplified 
supersonic  area  rule.  Further,  because  the  simplified  area  rule  is 
proposed  for  use  within  the  context  of  the  introduction,  a  qualitative 
result,  based  on  Its  use  for  analysis  of  near-axls  protuberances,  se«  i.s 
strengthened  by  Hall's  results. 

IV .  Computer  Program 

A.  Numerical  Integration 

A  computer  program  was  written  to  numerically  Integrate  Eqn.  (7) 
by  dividing  the  axis  Into  a  number  of  small  segments  of  length  Let 

us  assume  that  we  have  an  axial  distribution  of  Mach-plane  areas  S(x,M) 
(the  methods  I  used  to  arrive  at  this  distribution  will  be  discussed 
later)  at  the  appropriate  Mach  number,  K.  The  numerical  Integration  of 
Eqn.  (7)  was  carried  out  using  the  following  simple  algorithm: 


Figure  2.  Experimental  Drag  Rise  Comparison  for  the  Aircraft 

Configuration  and  the  Body  of  Revolution.  Drag  Based  on  Wing 
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The  summation  was  broken  up  and  an  analytic  expression  used  at  i-j  to 
avoid  the  singularity  In  the  natural  logarithm  as  x  approaches  £  (i.e., 
i-j)  • 

Because  the  goal  of  writing  the  program  Included  making  the  Input 
format  as  simple  as  possible,  the  expected  Input  format  was  a  simple 
tabulation  of  the  axial  location  and  cross-sectional  area  of  the  airplane 
configuration.  Although  the  number  of  locations  chosen  would  affect  the 
accuracy  of  the  numerical  Integration  (this  point  will  be  discussed 
later),  a  more  complex  problem  arises  from  the  fact  that  Eqn.  (10) 
requires  the  second  derivative  of  the  area  schedule  at  each  axial 
location,  be  It  the  original  area  schedule  or  the  Mach-plane  area 
schedule.  Several  schemes  for  dealing  with  this  problem  were  tried. 
Including  fitting  the  area  schedule  with  a  cubic  spline  which  provided  a 
continuous  area  distribution  as  well  as  continuous  first  and  second 
derivatives.  The  cubic  spline,  however,  turned  out  to  be  the  source  of 
large  errors.  I  finally  settled  on  a  five-point  nested  average  of  simple 
differences  to  obtain  the  approximation  for  the  required  derivatives. 

Starting  with  the  ares  distribution,  a  range  of  five  slopes 
surrounding  a  given  axis  location  was  averaged  to  arrive  at  a  first 
derivative  value  for  that  point.  The  same  procedure  was  again  used  to 
calculate  the  second  derivative  from  the  calculated  first  derivative 
table.  Each  axial  location  was  handled  In  this  manner  except  the  extreme 
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ends  of  the  body  where  a  five-point  range  was  not  available.  In  these  two 
regions  progressively  smaller  slope-averaging  ranges  were  used. 

B.  Mach-Plane  Area  Distribution 

By  far  the  greatest  computational  simplification  cornea  as  a 
consequence  of  the  method  of  initially  describing  the  airplane 
configuration  as  a  simple  tabulation  of  cross-sectional  area.  Any 
algorithm  for  determining  equivalent  Mach-plane  area  need  only  deal  with  a 
body  of  revolution,  and  this  at  zero  incidence.  The  first  consequence  is 
that  the  Mach-plane  area  need  only  be  calculated  once  for  each  Mach 
number,  because  the  axial  symmetry  leads  to  the  same  area  for  every  roll 
angle  e  (taken  with  respect  to  the  body  axis).  Secondly,  again  because  of 
symmetry,  only  half  the  area  need  be  investigated,  l.e.,  0£<t>£l80  degrees, 
where  <P  defines  the  radius  location  of  the  forward-normal  projection  of 
the  Mach-plane  area  (see  Figure  5). 


Figure  5.  Geometric  Description  for  Computation  of  the  Radii  of  the  New 
Mach-plane  Body  of  Revolution  Applied  to  a  Simple  Body  of 
Revolut ion 

One  hundred  and  eighty-one  such  radii  were  determined  corresponding  to 
equal  &  ♦  of  one  degree. 

It  is  possible  to  perform  this  operation  in  two-dimensional  space. 

It  can  be  shown  that  an  angle  $  for  the  Mach-plane,  as  shown  in  Figure  5, 
is  related  to  r  in  the  two-dimensional  plane  of  Figure  6  by  the  following 
relationship: 
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tany  «  tan(y)  /slr»<t> 


(11) 


Figure  6.  Relation  between  ^  and  y 

By  using  the  two-dimensional  plane  technique  d>  was  chosen  at  one-degree 
Increments  through  a  total  of  it  radians  (180  degrees)  and  the 
corresponding  y’s  were  used  to  determine  the  appropriate  radii  of  the 
forward-normal  Mach-plane  projection  (i.e.,  rj,  ,  etc.  of  Figure  5)  at 
each  value  of  x.  The  area  for  each  pie  slice  between  adjacent  radii  was 
determined  and  added  to  the  other  slice-shaped  areas.  This  gave  the  area 
of  half  the  Mach-plane  projection.  The  total  area  was  obtained  by 
multiplying  by  two. 

The  intersection  of  the  line  at  sn  angle  y  and  the  radius  of  the  body 
of  revolution  was  determined  by  marching  to  the  right  when  the  axis  system 
was  laid  out  as  in  Plgure  5  (or  to  the  left  If  41  was  less  than  0)  in 
Increments  of  6x;  computing  the  height  of  the  y-leg  of  the  triangle 
determined  by  the  Y-angle  ray,  x  marching  distance  and  the  height,  y;  and 
comparing  the  height  with  the  radius  at  that  location.  Once  the  radius 
had  been  exceeded,  an  iterative  proceaB  was  used  to  let  y  approach  the 
linear  Interpolated  r  to  within  some  small  predetermined  amount. 


V .  Results 

A.  Computer  Storage  and  Process  Time 
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The  storage  requirement  is  dependent  upon  the  method  of  lnputlng 
the  original  area  distribution  and  the  manner  of  manipulating  the  various 
quantities  in  the  program;  however,  with  little  effort  a  program 
incorporating  the  algorithm  of  the  last  section  can  be  made  to  fit  into  a 
computer  with  32K  words  of  storage  (8-blt  words).  The  driving  Influence 
is  the  number  of  axial  locations  chosen.  One  hundred  axial  locations  gave 
sufficient  resolution  for  most  area  configurations. 

Process  time  varied  with  the  number  of  axial  points  chosen.  For  a 
program  written  in  Fortran  and  executed  on  a  Harris  500  computer,  the 
times  were  less  than  half  a  minute  for  50  points,  approximately  a  minute 
and  a  half  for  100  points,  and  approximately  6  minutes  for  200  points. 
These  are  the  process  times  to  determine  the  wave-drag  coefficient  of  a 
given  configuration  for  one  Mach  number.  Clearly,  there  is  an  advantage 
to  staying  near  100  points  in  terms  of  both  accuracy  and  execution  time 
considerations. 

B.  Controls 

As  mentioned  earlier,  Figures  2,  3,  and  A  represent  a  set  of 
controls,  in  an  experimental  rather  than  an  analytic  sense,  which  serves 
as  a  Justification  for  the  use  of  the  method;  however,  it  is  of  interest 
to  explore  the  method  further.  In  order  to  investigate  the  sensitivity  of 
the  method  for  determining  the  best  location  to  place  a  protuberance  of 
known  volume  and  shape,  1  performed  a  wind-tunnel  test  of  a  basic  body  of 
revolution  which  could  change  configurations  so  as  to  place  a  known 
protuberance  in  two  locations  on  the  body.  Then  the  same  configurations 
were  analyzed  using  the  computer  program,  and  the  analytical  results 
compared  to  the  wind-tunnel  data. 
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A  brass  model,  the  pieces  of  which  are  shown  in  Figure  7,  was 
constructed  so  that  it  allowed  for  interchanging  the  protuberance 
locat ion . 


Figure  7.  Wind-Tunnel  Model  Showing  the  Various  Parts  Which  Made  the 
Configuration  Changes  Possible 

The  basic  body  was  a  parabolic  (single-arch)  body  of  revolution,  ten 
inches  long,  whose  radius  in  inches,  r,  versus  axial  position  in  inches, 
x,  is  described  by  the  equation 

r  -  (0.225  -  0.01875x)x  (12) 

A  two-inch  long,  f oo t b a  1 1 - sha pe d  protuberance  was  added  to  the  body 
of  revolution,  in  one  case  at  x  =  2  inches  and  in  the  second  case  at  X“7 
inches.  The  radius  of  the  protuberance  in  inches,  rp,  versus  axial 
location  along  the  protuberance  in  inches,  f.  ,  is  given  by  the  equation 

r  =  (0.8  -  0.4CK  (13) 

P 

Each  of  these  configurations  was  tested  at  four  Mach  numbers  (1.2, 
1.67,  2.45,  and  3.02)  in  the  one-foot  square  test  section  of  the  trisonic 
wind  tunnel  at  the  United  States  Air  Force  Academy.  This  tunnel  has  been 
described  in  detail  previously  (Ref.  8).  Schlieren  photographs  of  the 
three  configurations  typical  of  ail  runs  are  shown  in  Figure  8.  These 
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photographs  show  that  the  flow  remained  attached  and  Indicate  that 
differences  in  drag  should  be  due  to  wave-drag  differences. 

The  results  are  shown  in  Figure  9  along  with  the  predicted 
performance,  predicted  using  the  s irapl i f ied-ar ea-r ul e  program. 

The  method  does  seem  to  be  able  to  judge  the  best  location  for  a  given 
protuberance,  as  demonstrated  for  these  configurations. 

The  demonstrated  ability  of  the  method  to  discriminate  the  best 
location  of  a  given  protuberance  must  be  attributed  to  the  supersonic  area 
rule,  since  the  test  involved  a  body  of  revolution  for  which  the 
simplification  is  identical  to  the  rule. 

In  order  to  further  exercise  the  method  and  the  program,  I  tested 
three  model  configurations  originally  studied  by  Nelson  and  Welsh  (Ref. 

2).  1  chose  models  1,  4,  and  5  from  Ref.  2  for  the  following  reasons:  (a) 

model  1  was  a  simple  body  of  revolution  for  which  the  program  should  give 
results  identical  to  the  standard  supersonic  area  rule;  (b)  model  4 
represented  an  example  of  moderately  large  areas  located  far  off  axis  for 
which  the  supersonic  area  rule  did  a  good  job  of  predicting  the 
experimental  wave  drag;  and  (c)  model  5  represented  an  extreme  case  of 
large  areas  located  far  off  axis  for  which  even  the  supersonic  area  rule 
began  to  fail.  These  models  are  described  in  detail  in  Ref.  9  and 
summarized  briefly  in  Table  1. 

The  results  of  the  comparison  are  shown  in  Figure  10.  Notice  that 
both  the  simplified  and  complete  area  rule  give  identical  results  for  the 
body  of  revolution  --  as  they  should.  Model  4  shows  that  the  area  rule 
predicts  the  data  very  closely,  while  the  simplified  area  rule  is  off  by 
a  8  much  as  20  percent  at  M-l.l,  but  is  within  7  percent  by  M“1 . 5 .  For 
model  5  it  is  clear  that  even  the  area  rule  is  unacceptable,  and  the 
simplified  area  rule  is  worse. 
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It  should  be  noted  that  In  those  cases  where  the  supersonic  area  rule 
gives  good  results,  the  simpl 1 f ic at  ion  does  a  fair  job.  Further  note  that 
this  Is  the  total  drag,  and  the  Intended  use  of  the  simplification  la  to 
give  an  estimate  of  drag  Increments  due  to  minor  configuration  changes. 

C.  Application  to  the  F-15 

The  preceding  section  represented  a  set  of  what  might  be  termed 
controls  but  did  not  enable  me  to  decide  the  usefulness  of  the  simplified 
method  for  applications  which  fall  under  the  stated  purpose  of  the 
Introduction.  In  this  section  the  actual  and  predicted  effect  of  adding 
near-axls  protuberances  to  the  basic  clean  F-15  will  be  examined.  Since 
some  wind-tunnel  data  exists  for  the  McDonnel 1-Dougla s  FAST  PACK  fuel 
pallets,  these  have  been  selected  for  examination  (Ref.  10).  Figure  11 
shows  the  location  of  the  left  pallet  on  the  F-15.  Figure  12  gives  the 
area  distributions  (with  and  without  the  pallets)  used  for  the 
computations. 

First,  the  area  distribution  for  the  clean  F-15  was  used  to  obtain 
the  coefficient  of  drag.  Then  the  distribution  was  modified  to  account 
for  the  additional  area  due  to  the  conformal  pallets.  (Conformal  means 
that  the  added  body  Is  designed  to  blend  with  the  main  body  so  aa  to 
reduce  undesirable  aerodynamic  characteristics  such  as  flow  separation, 
which  can  Increase  drag.)  The  curve  In  Figure  13  was  obtained  by 
subtracting  the  wlthout-pallet  configuration  results  from  the  wlth-pallet 
results.  In  this  case  the  simplified  area  rule  more  closely  predicted  the 
data  than  did  the  area  rule. 

VI .  Conclusions 

In  this  paper  I  have  presented  the  theory  and  application  of  a 
modified  supersonic  area  rule.  Although  the  modification  Is  applied  to  a 
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•Ingle  example.  It  appears  to  lead  to  acceptable  results  when  used  for  the 
express  purpose  of  obtaining  drag  Increments  due  to  adding  near-axls 
protuberances  to  existing  aircraft  configurations.  Further,  the  cases 
examined  for  the  purpose  of  providing  a  set  of  controls  seem  to  Indicate 
that  the  modification  to  the  supersonic  area  rule  yields  relatively  good 
predictions  of  the  total  wave  drag  in  those  cases  where  the  unmodified 
area  rule  also  yields  good  results.  This  result  exceeds  the  expectations 
of  the  modification. 

There  Is  no  doubt  that  the  merit  of  the  modification  presented  here 
needs  further  study  before  receiving  general  wide-scale  approval.  The 
results  of  this  study,  however,  do  seem  to  indicate  that  this  modification 

to  the  supersonic  area  rule  would  prove  very  useful  in  systems  studies  and 

ear  1 y-management  decisions.  Therefore,  I  recommend  its  use,  since 
employing  an  approximation  to  the  supersonic  area  rule  is  better  than  not 
exploiting  the  area  rule  at  all. 

VI I  Epl 1 ogue 

In  preparing  Figure  13  for  this  paper  I  was  struck  by  how  well  the 
modified  supersonic  area  rule  predicted  the  data  when  dealing  with  drag 
increments.  I  thought  that  there  might  be  something  Inherent  in  the 

method  which  uncoupled  the  area  increments  from  the  base  area  distribution 

and  that  this  might  might  be  demonstrated  mathematically.  To  this  end  I 
expanded  Eqvi.  (10),  including  the  difference  schemes,  to  see  If  there 
might  be  a  way  in  which  nonl 1 nea r 1 t i e s  were  somehow  compensating  for  one 
another.  I  found  that  there  was  still  a  strong  nonlinear  coupling  term 
which  linked  the  additional  protuberance  area  to  the  original 
configuration  area.  Based  on  this  cursory  analysis  1  must  conclude  that 
the  agreement  demonstrated  In  Figure  13  Is  only  by  chance  better  than  that 
determined  by  the  full  supersonic  area  rule.  That  Is  to  say,  the 
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aodlf lcat Ion  approximates  the  full  area  rule  and  in  general  one  would 
expect  that  the  superaonlc  area  rule  foraa  the  lialtatlon  of  how  accurate 
one  would  expect  reaulta  to  be,  and  the  modification  will  approxlaate  this 
solution.  This  question,  however,  aight  warrant  further  study. 
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J69-T-25  TURBOJET  CYCLE  ANALYSIS 
Robert  E.  Boyle,  Jr.* 

Abstract 

Numerous  test  runs  of  the  J69-T-25  turbojet  engine  have  been  made  at 
the  Aeronautics  Laboratory  of  the  United  States  Air  Force  Academy  over  the 
past  25  years.  However,  the  test  data  —  specifically  the  nominal  values 
of  thrust,  specific  fuel  consumption,  and  compressor  and  turbine 
efficiencies  —  have  not  been  readily  accessible.  To  correct  this 
situation,  this  paper  makes  the  data  available  in  a  comprehensive 
reference.  Since  the  instructor  or  researcher  will  be  provided  with 
sufficient  information  for  planning  purposes,  needless  engine  runs  will  be 
prevented  . 

Two  analyses  of  the  data  were  produced.  The  first  is  a  simplified 
analysis  based  on  assumptions  that  make  the  de  t  er  mi  na  t  ijpn  of  engine 
performance  parameters  less  complicated.  This  provides  a  basic 
understanding  of  the  engine  operation  while  reducing  the  amount  of  data 
reduction  required.  Only  two  simplifying  assumptions  are  made  in  the 
second  analysis;  thus  the  determination  of  engine  operating  parameters  is 
more  accurate  but  more  complicated. 

The  simplified  analysis  yields  results  that  are  within  10  percent  of 
the  experimental  values,  demonstrating  that  nominal  values  can  be 
predicted  with  a  less  sophisticated  approach.  This  analysis  lends  itself 
nicely  to  demonstrating  the  principles  of  the  turbojet  engine  to  a  basic 
engineering  class.  The  rigorous  analysis,  on  the  other  hand,  provides  an 
adequate  description  of  what  actually  occurs  during  engine  operation  and 
could  be  used  in  more  advanced  courses  on  jet  engine  design. 


I .  I nt  roduc  t 1  on 

Teledyne  J69-T-25  turbojet  engines  are  used  by  the  United  States  Air 
Force  as  power  plants  on  T-37  trainers,  as  the  booster  engines  on  the 
C-123J  cargo  aircraft,  and  in  the  MA1A  air  cart  to  supply  air  for  starting 
certain  Jet  aircraft.  Two  of  these  enggines  are  mounted  in  test  cells  at 
the  Aeronautics  Laboratory  of  the  USAF  Academy.  Their  primary  purpose  is 
to  demonstrate  to  Aeronautical  Engineering  students  the  principles  of 
turbojet  engines  and  the  operation  of  specific  engine  components.  By 
indicating  nominal  values  of  thrust,  fuel  flow,  and  turbine  and  compressor 
efficiencies,  this  paper  will  aid  instructors  who  demonstrate  the 
performance  of  these  engines  in  classroom  situations. 

Figure  1  shows  a  schematic  of  the  J69-T-25  turbojet  engine  and  the 

*Captaln,  USAF,  Associate  Professor  of  Aeronautics,  DFAN 
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Two  different  analyses  of  the  engine  performance  are  produced.  The  first 
analysis,  which  is  used  In  the  Department  of  Aeronautics  core  course, 
makes  several  simplifying  assumptions.  The  first  of  these  assumptions  is 
that  the  mass  flow  into  the  engine  is  equal  to  the  mass  flow  exiting  the 
engine.  This  is  true  only  *  *he  mass  flow  of  the  fuel  is  ignored.  It  is 
also  assumed  that  the  mechanical  efficiency  of  the  engine  is  unity.  This 
assumption  implies  that  there  are  no  frictional  or  mechanical  linkage 
losses  involved  in  the  operation  of  the  c ompr e s sor- t ur bi ne  combination. 

It  also  ignores  any  power  loss  due  to  the  operation  of  accessories  such  as 
the  fuel  pump,  the  oil  pump  or  the  generator.  The  pressure  loss  in  the 
combustor  is  assumed  to  be  negligible.  The  validity  of  this  assumption, 
which  affects  the  calculation  of  turbine  efficiencies,  will  be  examined 
using  the  experimental  data.  The  nozzle  is  assumed  to  be  adiabatic;  this 
assumption  will  also  be  examined  in  the  experimental  analysis.  The  values 
of  the  specific  heat  at  constant  pressure,  C  ,  and  ratio  of  specific 
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heatstY,  are  assumed  Co  be  constant  throughout  the  engine.  The  high 
temperatures  attained  as  well  as  change  In  composition  of  the  gas  in  the 
combustor  and  turbine  will  affect  the  accuracy  of  this  assumption.  Since 
the  velocity  of  the  flow  at  the  turbine  exit  is  assumed  to  be  negligible 
compared  to  nozzle  exit  velocity,  the  static  temperature  at  station  5  Is 
equal  to  the  total  temperature  there  or  T5  =  tq5‘  It  is  also  assumed  that 
the  exit  nozzle  is  perfectly  expanded,  i.e.,  exit  pressure  Is  just  equal 
to  atmospheric  pressure.  Since  the  Mach  number  at  the  exit  can  be  shown 
to  be  subsonic  at  all  RPM  settings,  this  assumption  is  accurate.  Also 
Inherent  In  a  static  thrust  analysis  is  the  fact  that  the  total  pressure 
and  temperature  at  the  entrance  of  the  compressor  Is  just  equal  to 
atmospheric  values.  With  the  low  entrained  velocities  in  the  inlet,  the 
changes  in  total  properties  are,  in  fact,  negligible.  The  assumptions  are 
therefore  : 


(1)  Hass  flow  rate  into  the  engine  equals  the  mass  flow  out 

of  the  engine  (m.  «  m  ),  i.e.,  the  mass  flow  rate  of 

in  out 

•  'V 

the  fuel  is  assumed  to  be  negligible  (  m  =  0) 

(2)  Mechanical  efficiency  of  unity  Or]  n|  *  1002) 

(3)  No  pressure  loss  in  combustor  (1>  =  I’  ,) 

U<4  U  i 

(4)  Adiabatic  nozzle  (T_,  =  T'  .) 

Ob  Ob 

(5)  C  and  Y  are  constant  throughout  the  engine 

P 

<C  -  .24  Y  =  1.4) 

p  lbm  R . 

(6)  Velocity  at  the  turbine  exit  is  much  less  than  nozzle 

exit  velocity  ( V.  <<  V  ) 

3  6 

(7)  Perfectly  expanded  nozzle  (Pft  =  P.jtn)) 

(8)  Total  pressure  and  total  temperature  in  the  inlet  are 
assumed  constant,  i.e.,  PQ7  =  patm>  T02  =  Tatm  *  This  is 
only  true  for  a  stationary  engine  with  no  losses  in  the  inlet. 
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(9)  The  air  and  fuel-air  mixture  are  assumed  to  be  perfect 
gases  (  P  =  .  RT  ) 

The  more  rigorous  approach,  which  is  taken  in  the  Aero  majors'  courses, 
uses  only  the  static  thrust  assumption  (no.  8)  and  the  perfect  gas 
assumption  (no.  9).  This  approach  better  describes  what  actually  occurs 
in  the  engine  and,  therefore,  should  produce  a  more  accurate  analysis  of 
component  performance.  Throughout  this  paper  the  engine  parameters  for 


both  analyses  will  be  compared. 


II.  Test  Apparatus 


The  Teledyne  J69-T-25  engine,  designed  in  1956,  is  very  similar  to 
the  original  turbojet  designs.  It  consists  of  a  centrifugal  compressor, 
an  annular  combustor,  an  axial  flow  turbine,  and  an  exhaust  notzle  of 
fixed  exit  area.  A  Bell-mouth  nozzle  has  been  added  at  the  inlet  to  allow 


steady,  one-dimensional  air  flow  to  be  delivered  to  the  compressor  and  to 
permit  an  accurate  determination  of  the  mass  flow  entering  the  engine. 

See  Figure  2  for  the  location  of  the  Bell-mouth  on  the  engine. 


Figure  2.  Picture  of  J69-T-25  Engine,  Bell-mouth  inlet, 
and  Instrumentation 
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PARAMETER 

SYMBOL 

DIMENSIONS 

Engine  Thrust 

F 

lbf 

Fuel  Flow 

■f 

lbm/hr 

Atmospheric  Pressure 

P 

atm 

in.  Hg 

Dynamic  Pressure 

in  Bell-mouth  (P  -  P  ) 

atm  2 

AP 

in.  H 20 

Total  Pressure  at 

Compressor  Entrance 

P02 

p  si 

Total  Pressure  at 

Compressor  Exit 

P03 

psl 

Total  Pressure  at 

Turbine  Exit 

P04 

psi 

Total  Pressure  at 

Nozzle  Exit 

P06 

in .  Hg 

Total  Temperature  at 
Compressor  Inlet 

T02 

°F 

Total  Temperature  at 
Compressor  Exit 

T03 

°F 

Total  Temperature  at 

Turbine  Exit 

T05 

°F 

Total  Temperature  at 

Nozzle  Exit 

T06 

°F 

Percent  of  Maximum 

Engine  RPM 

Z  RPM 

Z 

Actual  Engine  RPM 

RPM 

revolutions 
per  minute 

Table  1  show*  the  parameters  neasured  during  our  test  runs.  The 
teaperatures  in  the  engine  were  neasured  using  thermocouples  Implanted  in 
pitot  tubes.  Tq^  was  neasured  by  a  single  probe  at  the  inlet.  T^  was 
neasured  by  a  single  probe  located  at  the  compressor  exit. 


was  not 
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■  enured  ■  Although  measurement  of  thl*  temperature  by  aeane  of 
theraocouplee  is  possible,  the  failure  rate,  as  well  as  the  nuaber  of 


probes  required  to  accurately  aeaaure  t  .  make  It  econoalcally  attractive 


04 


to  calculate  the  turbine  Inlet  teaperature  analytically.  was  aeasured 


by  a  total  of  24  probes  mounted  radially  In  a  6  x  4  configuration  as  shown 
In  Figure  3. 


•  PROBES 


Figure  3.  Probe  Arrangement  for  Measurement  of  T 
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Tn.  Is  aeasured  at  the  noitle  exit  by  eight  probes  mounted  radially. 
Pressures  were  aeasured  at  different  locations  using  static  ports  or  pitot 
tubes.  In  soae  locations  aultlple  pressure  ports  were  used  and  the 
aeasured  pressures  were  averaged.  ?2  was  aeasured  using  four  static 
pressure  ports  aounted  radially  at  the  throat  of  the  Ball-aouth  and  Is 
used  to  determine  AP,  the  difference  between  the  total  pressure  and  the 
static  pressure  at  the  throat.  p  was  measured  by  one  probe  located  at 
the  compressor  exit.  Pq^  was  aeasured  by  two  probes  located  at  the 
turbine  entrance.  p  was  aeasured  by  six  probes  aounted  radially  at  the 
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nozzle  exit.  The  engine's  thrust  vas  measured  with  a  lead  cell  mounted 
at  the  base  of  the  test  stand. 

An  oil  pump,  a  fuel  pump,  and  an  engine  generator  are  the  only 
accessories  driven  by  the  engine.  The  power  required  to  operate  these 
engine  accessories  Is  Included  In  the  mechanical  efficiency,  >jm ,  used  In 
the  majors'  courses.  The  value  of  T)m  Is  .96  for  the  majors’  analysis 
and  1.0  in  the  simplified  core  course  analysis. 


Ill .  Analytical  Approach 

The  thrust  of  a  turbojet  engine  can  be  determined  analytically  by 
using  a  momentum  balance  on  a  control  volume  surrounding  the  engine.  If  a 
one-dimensional,  steady  flow  is  assumed  to  enter  and  leave  the  engine,  the 
thrust  can  be  represented  as 

P  '  \\  -  4>V,  +  *6  <P6  '  Pl>  <‘> 

where  F  Is  the  thrust  of  the  engine.  If  the  nozzle  Is  perfectly  expanded 
for  all  RPM  settings,  then  P&  -  -  Patm-  Also,  this  analysis  Is  a 

determination  of  thrust  for  a  stationary  turbojet  engine,  l.e.,  static 
thrust,  so  the  Inlet  velocity  is  essentially  zero  ( •'  0  ) ■  The  resulting 
thrust  equation,  then.  Is  simply 

F  -  m6V6  (2) 


where  m.  Is  the  mass  flow  out  of  the  engine  and  V,  is  the  velocity  at  the 

6  o 

nozzle's  exhaust.  The  mass  flow  and  velocity  at  the  nozzle  exit  are  not 
easily  measured  during  engine  operation,  so  Eqn.  (2)  must  be  altered  to 
Include  known  parameters.  One  form  of  this  equation  Is 


F 


m | ( 1+f ) 
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where  f  is  the  fuel-air  ratio  (mf/m  ).  For  the  core  course  thrust 
calculation,  Eqn.  (3)  reduces  to 


F  = 


2C 

P 


2 


(4) 


since  the  fuel  flow  is  neglected  ( “  0)  and  the  mass  flow  out  of  the 
engine  is  assuned  to  be  equal  to  the  inlet  mass  flow. 

It  is  clear  from  Eqn.  (4)  that  both  C  and  y  affect  the  calculation 

P 

of  engine  thrust.  No  appreciable  change  in  these  two  parameters  is 
encountered  in  the  inlet  or  compressor  due  to  the  relatively  low 
temperature  changes  encountered  there.  Thus,  they  are  assumed  to  be 
constant  up  to  the  combustor.  These  values  are:  Yc  *  1.4,  Cpc  “  •  24  y^-jj  • 

In  the  simplified  analysis  these  parameters  are  assumed  to  be  constant 
throughout  the  engine  regardless  of  the  temperature.  In  the  majors' 
course  analysis,  the  variation  with  temperature  is  taken  into  account  by 
assuming  that  Y  and  Cp  are  constant  after  the  combustion  process  but  with 


.  The  subscript  "t"  is  usually 


new  values  Y,  ■  1.35  and  C  -  .264  ■  - 

1  pt  lbm  R 

used  to  denote  the  values  of  and  C  downstream  of  the  combustor.  Later, 
another  calculation  of  theoretical  thrust  is  discussed  in  which  the  values 
of  and  C  are  varied  in  the  turbine  and  nozzle. 

Isentroplc  and  polytropic  efficiencies  of  the  compressor  and  turbine 
are  very  important  to  any  Jet  engine  analysis.  The  isentroplc  efficiency 
is  one  way  of  measuring  the  amount  of  power  lost  due  to  the  entropy 
produced  during  the  compression  or  expansion  of  the  flow.  In  other  words, 
we  compare  an  actual  device  to  an  ideal  or  the  best  possible  adiabatic 
device.  For  example. 


ideal  power  required 
'c  m  actual  power  required 


(5) 
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where  As  is  entropy  produced  and  R  is  the  gas  constant  (Ref.  5).  If  this 
pressure  change  is  analyzed  differentially,  i.e.  the  total  change  is 
divided  into  many  small  increments  of  pressure  change,  Eqn.  (11)  can  be 
represented  as 


03 

02 


(12) 


In  other  words,  we  are  replacing  the  entropy  production  term,  which  can 
vary  significantly  from  point  to  point  in  the  compressor,  with  a  term 
which  is  easily  determined  and  remains  practically  constant  throughout  the 
compressor.  The  polytropic  efficiency  of  the  turbine,  ,  can  be 
determined  from  the  following  relation: 


IV .  Results 

Two  test  runs  of  the  J69-T-25  turbojet  engine  on  8  July  1982  produced 
the  data  used  in  this  report.  Table  2  shows  the  average  values  of  the 
data  taken  for  the  two  runs,  while  Table  3  lists  the  absolute  values  of 
these  same  temperatures  and  pressures,  and  Table  4  gives  operating 
parameters  calculated  from  the  variables  in  Table  3. 

The  actual  mass  flow  into  the  engine  inlet  was  determined  from  the  P 
measured  in  the  Bell-mouth  (Ref.  7)  and  is  tabulated  in  Table  4  for 
various  engine  RPM  (Ref.  1).  The  fuel-air  ratio,  determined  from  the  air 
mass  flow  calculated  at  the  inlet  and  the  measured  fuel  flow,  is  also 
listed.  The  Mach  number  at  the  engine’s  exit,  M  ,  is  tabulated  in  Table 

D 

4  along  with  the  calculated  value  of  turbine  inlet  temperature,  TQ^  .  The 

determination  of  T-.  will  be  discussed  later.  The  variation  of  the 

04 
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lsentroplc  and  polytropic  efficiencies  of  the  turbine  and  compressor  as  a 
function  of  RPM  tabulated  In  Table  4  are  plotted  In  Figures  4  and  5. 


%  RPM 

Figure  4.  Isentroplc  Efficiency  of  Compressor  and  Turbine  vs.  RPM 


Figure  5.  Polytropic  Efficiency  of  Compressor  and  Turbine  vs.  RPM 
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Notice  the  fairly  constant  values  of  ec  above  70  percent  RPM.  Thus  ec 
behaves  as  expected.  The  isentropic  efficiency  of  the  compressor  peaked 
at  about  90  percent  RPM  indicating  the  desired  cruise  operating  point  for 
the  engine.  Figure  6  shows  the  different  turbine  efficiencies  obtained 
from  the  two  analyses.  The  relatively  large  difference  occurs  because 


Figure  6.  Comparison  of  Isentropic  Turbine  Efficiencies  from  the  Two 
Analyse  s 


the  pressure  drop  in  the  burner  yields  a  turbine  inlet  pressure  lower  than 
that  assumed  in  the  core-course  analysis.  Because  of  this  lower  pressure 
the  related  ideal  temperature  drop  is  lower  than  originally  predicted, 
resulting  in  higher  turbine  efficiencies  than  those  obtained  from  the  core 
analysis.  T_^  was  not  measured  but  was  calculated  using  a  power  balance 
between  the  compressor  and  turbine.  A  mechanical  efficiency  of  96  percent 
was  assumed  in  the  majors'  course  analysis.  This  efficiency  takes  into 
account  not  only  frictional  losses  but  also  the  power  required  to  drive 
the  engine  accessories,  as  previously  mentioned.  For  example. 
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These  parameters  are  also  corrected  to  sea  level  as  shown  in  the  figures. 
Figure  9  shows  the  comparison  of  measured  thrust  to  those  values 
determined  by  the  two  analyses.  The  corrected  sea-level  thrust  is 
approximately  8.5  percent  less  than  the  original  engine  data  (Ref.  6). 
This  is  probably  due  to  the  deterioration  of  the  engine  components  over 
the  past  25  years.  Figure  10  shows  the  thrust  specific  fuel  consumption 

(TSFC)  as  a  function  of  RPM. 

12001 


%  RPM 

Figure  9.  Comparison  of  Calculated  to  Measured  Thrust 

2.2s| 


2.0 


1.75 


I  5 


1.5 


“  1-25 


1.0 


O  MEASURED  TSFC 


G 


60 


-L. 


70 


80 


90 


100 


%  RPM 

Figure  10.  Thrust  Specific  Fuel  Consumption  vs.  RPM 
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For  the  core  course  T  calculation,  Eqn.  (15)  reduces  to 

*04  =  f03  +  *05  “  *02 


(15) 


since  c  =  C  ,  V  ■  1 ,  and  m  =  m,  .  As  previously  mentioned,  these 

pc  pc  m  lb 

calculated  values  of  T  for  both  core  and  majors'  calculations  are 

(M 

listed  in  Table  4.  Thrust  and  fuel  flow  are  shown  as  a  function  of  engine 
RPM  in  Figures  7  and  8. 


%  RPM 

Figure  7.  Measured  and  Corrected  Sea  Level  Thrust  vs.  RPM 


%  RPM 

Figure  8.  Measured  and  Corrected  Sea  Level  Fuel  Flow  vs.  RPM 
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In  the  present  engine,  at  the  altitude  of  the  Academy,  the  minimum  TSFC 
occurs  at  100  percent.  When  corrected  to  sea  level  the  minimum  occurs  at 
approximately  95  percent.  The  sea  level  TSFC  indicates  almost  no  shift  in 
the  original  sea  level  data  produced  by  the  engine  manufacturer. 

V .  Conclusions 

The  simplified  analysis  gives  a  thrust  within  4  percent  of  the 

measured  value  as  shown  in  Figure  9.  This  apparently  good  prediction  Is 

due  to  compensating  errors  introduced  by  the  assumptions.  For  example, 

the  adiabatic  nozzle  assumption  requires  the  use  of  Instead  of  the 

lower  T  ,  •  This  causes  a  higher  T  ,  a  higher  V.  ,  and,  therefore,  a 
06  6  6 

predicted  thrust  higher  than  that  in  the  majors'  analysis.  On  the  other 
hand,  the  core  analysis  does  not  consider  the  mass  flow  of  the  fuel  in  the 
thrust  calculation.  This  results  in  a  predicted  thrust  lower  than  that  of 
the  majors'  analysis-  These  two  assumptions  result  in  an  accurate  thrust 
prediction  (see  Table  4). 

The  two  analyses  show  large  differences  in  their  predictions  of 
turbine  efficiencies.  This  Is  due  to  the  core  analysis  assumption  that 
there  is  no  pressure  loss  in  the  combustor.  The  lsentropic  efficiency  is 
calculated  using  and  ideal,  lsentropic  temperature  at  the  turbine  exit. 
This  ideal  temperature  depends  on  the  pressure  loss  in  the  turbine.  The 
pressure  at  the  turbine  inlet  is  assumed  to  be  higher  than  it  actually  is, 
resulting  in  a  T  gs  lower  than  it  should  be,  as  shown  by  Eqn.  (10).  Using 
Eqn.  (9)  one  can  see  that  a  lower  T7t  results.  The  turbine  lsentropic 
efficiency  should  be  inherently  higher  than  that  in  the  compressor  due  to 
their  modes  of  operation.  The  compressor  is  centrifugal,  causing  large 
entropy  changes,  while  the  axial  flow  turbine  accomplishes  its  task  more 
efficiently.  The  polytropic  and  lsentropic  efficiencies  in  the  turbine, 
as  expected,  are  very  similar  because  the  turbine  is  single  stage.  One 
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must  also  consider  that  the  polytropic  efficiency  of  the  compressor  may  be 
meaningless  in  that  ec  is  used  to  determine  stages  of  compression  during 
engine  design.  The  centrifugal  compressor,  by  its  very  nature.  Is  not 
staged  . 

An  analysis  was  completed  to  determine  the  effect  on  thrust  of 

variation  of  Y  and  C  with  temperature  at  the  turbine  inlet,  the  turbine 

t  Pt 

exit,  and  the  nozzle  exit.  This  resulted  in  no  more  than  0.1  percent 

difference  in  thrust.  The  assumption  that  y  =  1.35  and  C  =  .264  -l'~- 

t  pt  lbm  R 

and  that  they  are  constant  downstream  of  the  burner  is,  therefore,  quite 
adequate  for  this  analysis. 

The  turbine  inlet  temperature  was  almost  constant  below  80  percent 

RPM-  This  is  due  to  the  fact  that  the  fuel  air  ratio  is  nearly  constant 

in  this  range.  As  f  then  increased  toward  stoichiometric  (f  “  0.066) 

st 

the  turbine  inlet  temperature  also  increased.  The  turbine  exit 
temperature  showed  a  decrease  and  then  an  Increase  at  higher  RPM  because, 
as  the  engine  RPM  increased,  the  amount  of  work  required  by  the  turbine 
continually  Increased.  Since  T was  constant  at  first,  TQ,_  must 
decrease  (work  of  the  turbine  is  proportional  to  T ^  -  T ) .  As  T rose, 
T  also  increased  at  higher  RPM. 

Despite  its  somewhat  inaccurate  calculations,  caused  by  the  many 
restrictive  assumptions,  the  simplified  core-course  analysis  provides  a 
basic  understanding  of  the  turbojet  engine  for  the  non-Eng inee r i ng  general 
student.  The  majors'  course  analysis  is  much  less  restrictive  and 
provides  a  more  accurate  description  of  the  actual  component  operating 
parameters.  This  approach  is  necessary  for  the  Aero  Propulsion  student  to 
gain  the  more  sophisticated  background  necessary  for  a  propulsion  design 
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pc 

Cpt 
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P02 

P03 

P04 

P05 

P06 

T02 

T03 

T04 

T05 

T06 

\ 

Yc 


't 
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Symbol « 

inlet  frontal  area 
exhaust  nozzle  exit  area 

specific  heat  at  constant  pressure  ahead  of  combustor 

specific  heat  at  constant  pressure  aft  of  conbustor 

mass  flow  rate  Into  compressor 

mass  flow  rate  from  nozzle  exit 

mass  flow  rate  of  fuel 

Mach  number  at  the  exit 

total  pressure  at  compressor  entrance 

total  pressure  at  compressor  exit 

total  pressure  at  turbine  entrance 

total  pressure  at  turbine  exit 

total  pressure  at  nozzle  exit 

total  temperature  at  compressor  entrance 

total  temperature  at  compressor  exit 

total  temperature  at  turbine  entrance 

total  temperature  at  turbine  exit 

total  temperature  at  nozzle  exit 

mechanical  efficiency  of  compre ssor-turbine  combination 

ratio  of  specific  heats  ahead  of  combustor 

ratio  of  specific  heats  aft  of  combustor 

lsentropic  efficiency  of  compressor 

lsentropic  efficiency  of  turbine 

polytropic  efficiency  of  compressor 

polytropic  efficiency  of  turbine 
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ROUGHNESS  EFFECTS  IN  AXIAL  FLOW  COMPRESSORS: 

AN  EMPIRICAL  MODEL 

R.J.  Stiles* 

Abstract 

This  paper  presents  an  empirical  model,  taken  from  recent  technical 
literature,  which  can  be  used  to  compute  both  the  compressor  airfoil 
surface  finish  required  to  create  a  hy d r od y na m lc a  1 1 y  smooth  fluid  flow 
through  an  axial  flowjet  engine  compressor  and  the  effects  of  changes  In 
airfoil  surface  finish  on  compressor  efficiency.  Standard  surface 
roughness  measurements  and  the  problems  associated  with  using  those 
measurements  in  the  model  are  also  discussed.  To  illustrate  its  use,  the 
model  has  been  applied  to  the  compression  sections  of  three  Air  Force 
gas-turbine  engines:  the  TF33-PW-7,  the  F100-PW-100,  and  the  J85-GE-5H. 

I .  IntroJuctl. n 

In  recent  yea.s  rapidly  escalating  fuel  costs  have  created  renewed 
interest  among  aircraft  gas-turbine  engine  designers  in  the  effects  of 
compressor  airfoil  surface  finish  on  turbine  engine  performance.  There  is 
some  experimental  data  which  shows  that  reduced  roughness  in  compressor 
components  can  result  in  improved  compressor  efficiency  which,  in  turn, 
will  yield  lower  thrust  specific  fuel  consumption  (TSFC)  and  lower  exhaust 
gas  temperature  (EGT)  for  the  engine  (Refs.  1  through  3).  Since  current 
Air  Force  fuel  consumption  is  approximately  3.6  billion  gallons  per  year, 
even  a  small  reduction  In  TSFC  for  the  engines  which  consume  most  of  that 
fuel  can  result  In  substantial  cost  savings.  One  manufacturer  of  smooth 
coatings  for  compressor  components  has  suggested  that  the  application  of 
his  product  to  J57,  TF33,  and  J79  engines  could  yield  an  average  1.5 
percent  reduction  in  TSFC  for  those  engines;  that  would  save  as  much  as 
200  million  gallons  of  fuel  over  a  six  year  period  (Ref.  4).  Since  high 
EGT  is  a  primary  life-limiting  factor  for  the  materials  in  a  jet  engine, 
lower  EGT  should  result  in  lower  engine  maintenance  and  replacement  costs 
over  a  period  of  time. 

The  objective  of  this  paper  is  to  describe  and  illustrate  an 
empirical  model  developed  by  Schaffler  (Ref.  3)  which  can  be  used  to 

•Major,  USAF,  Aero  Propulsion  Laboratory 
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estimate  the  effect  of  airfoil  surface  finish  on  axial  compressor 
performance.  The  model  permits  designers  to  calculate  the  airfoil  surface 
finish  required  to  obtain  maximum  available  compressor  efficiency  at  a 
given  compressor  operating  point  (specified  by  pressure  ratio  and 
corrected  speed)  and  aircraft  mission  point  (specified  by  Mach  number  and 
altitude).  Also,  the  model  can  be  used  to  estimate  changes  In  compressor 
efficiency  associated  with  changes  in  airfoil  surface  finish.  This 
information  can  then  be  used  In  an  engine  cycle  analysis  to  determine 
changes  In  TSFC  and  EGT. 

The  discussion  that  follows  begins  with  a  definition  of  surface 
roughness  and  provides  some  background  information  concerning  the 
difficulty  of  measuring  roughness  in  a  given  material.  Then,  after  the 
empirical  model  is  presented,  its  use  is  Illustrated  by  application  to  the 
compression  sections  of  three  Air  Force  engines:  the  TF33-PW-7,  the 
F100-PW-100,  and  the  J85-GE-5H. 

I I .  Definition  and  Measurement  of  Roughness 

Surface  roughness  is  defined  as  the  "fine  Irregularities  in  the 
repetitive  or  random  deviations  from  the  Intended  surface  contour"  of  a 
given  material  (Ref.  5).  There  are  a  variety  of  surface  measurements  that 
can  be  made  to  characterize  a  material  surface  In  terms  of  Its  roughness 
(see  Ref.  5  for  examples),  but  a  measurement  called  the  arithmetic  average 
roughness,  Ra,  is  most  commonly  used.  As  Illustrated  In  Figure  1,  Ra  Is 
the  arithmetic  average  of  deviations  from  a  mean  line  through  the  surface 
profile.  It  is  defined  by  the  equation 


X-L 


X-0 
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where  y 
length . 


1 8  the  ordinate  of  the  curve  of  the  profile  and  L  Is  the  sampling 


t 


% 


- A - J^J. 


Ra 


Surface 
'  L 


Figure  1-  Nomenclature  for  Determining  the  Arithmetic  Average  Roughness 

While  It  is  relatively  easy  to  measure  various  physical 
characteristics  of  a  material's  surface,  determination  of  the 
relationship  between  the  measured  characteristic  and  fluid  flow  behavior 
over  the  surface  is  more  difficult.  In  fluid  mechanics  surface  roughness 
Is  often  specified  in  terms  of  "equivalent  roughness,"  a  concept  which  has 
its  basis  in  systematic  and  extensive  experiments  performed  by  Nlkuradse 
(Ref.  6)  and  which  is  symbolized  by  the  parameter  k  g.  Roughness  on  the 
inside  of  a  pipe  will  reduce  the  fluid  flow  rate  through  the  pipe  for  a 
given  pressure  gradient  along  the  length  of  the  pipe.  This  loss  can  be 
quantified  by  a  parameter,!,  called  the  resistance  coefficient  of  pipe 
flow  (see  Ref.  7,  page  79).  In  Nikuradse's  experiments,  sand  grains  of  a 
uniform  height,  ks,  were  applied  (as  closely  packed  as  possible)  to  the 
inside  surface  of  a  circular  pipe  of  a  given  radius,  R,  and  X  was  measured 
for  a  range  of  values,  ks,  R,  and  pipe  Reynolds  numbers.  For  a  particular 
Reynolds  number  range  the  flow  in  the  pipe  can  be  characterized  as 
"completely  rough"  (see  P.ef.  7,  page  5R0  for  a  complete  definition  of  this 
term),  and  1  is  found  to  agree  very  closely  with  values  given  by  Eqn.  (2). 


X  «  - - - 

(2  log  +  1 . 74) 2 
Ks 

A  surface  is  said  to  have  an  equivalent  sand  roughness  of  ks  it  the 
measured  value  of  X  for  flow  over  the  surface  la  equal  to  the  value  of  X 
computed  using  Eqn  (2)  with  the  value  for  ks.  Equivalent  sand  roughness 


(2) 
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for  a  particular  surface  depanda  on  the  shape,  height,  and  distribution  of 
surface  Irregularities  (or  roughness  clesents).  Figure  2,  taken  fro*  Ref. 
7,  page  587 ,  provides  a  suaiary  of  so*e  equivalent  sand  roughneas 
aeasuresents  aade  hy  Schllchtlng. 


Figure  2.  Equivalent  Sand  Roughness  for  Various  Roughness  Types  (from 
Boundary-Layer  Theory  by  Dr.  Hermann  Schllchtlng,  translated 
by  Dr.  J.  Kestin.  Copyright  (c)  1968,  McGraw-Hill,  Inc.  Used 
with  the  permission  of  McGraw-Hill  Book  Company) 

Note  from  Figure  2  that  the  ratio  of  roughness  element  spacing,  D,  to 

roughness  element  height,  k,  appears  to  be  an  important  factor  Influencing 

equivalent  sand  roughness;  that  Is,  Increasing  the  ratio  D/k  for  a  given 

[  surface  reduces  ka  for  that  surface. 
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If  the  ratio  D/k  Is  large,  the  roughness  Is  said  to  be  k-type  or 
protuberance  roughness,  while  If  the  ratio  D/k  is  snail,  the  roughness  Is 
said  to  be  d-type  or  cavity  roughness.  These  types  of  roughnesses  are 
illustrated  in  Figure  3.  An  experimental  study  by  Perry  et  al  (Ref.  8) 
showed  that  turbulent  boundary  layer  flow  is  significantly  different  for 
surfaces  with  cavity  roughness  than  it  is  for  surfaces  with  protuberance 


HbMwh 

nxy;a 


m  3i 


d-type  (cavity)  roughness 


Figure  3.  k-type  and  d-type  Roughness 
Surfaces  of  compressor  blades  and  vanes  contain  a  variety  of  peak  (or 
cavity)  shapes  and  a. broad  spectrum  of  peak  (or  cavity)  densities  and 
heights.  As  noted  by  Schaffler  (Ref.  3),  the  Ra  value  of  a  given  surface 
may  not  be  indicative  of  the  hydrodynamic  character  of  the  surface.  The 
effect  of  surface  roughness  is  felt  in  turbulent  flow  when  the  roughness 
elements  penetrate  the  laminar  sublayer  of  the  turbulent  boundary  layer. 

It  is  conceivable  that  three  surfaces  could  have  equal  Ra  values  and  yet 
could  have  significantly  different  effects  on  the  boundary  layer  flow.  In 
Figure  4,  surface  (a)  has  a  series  of  moderate  deviations  fron  the  mean 
line  with  no  pronounced  peaks  or  cavities.  Surface  (b)  is  quite  amooth 
with  occasional  pronounced  peaks,  and  surface  (c)  is  a  mirror  image  of 
(b).  Suppose  that  all  protuberances  are  contained  in  the  laminar  sublayer 
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of  the  boundary  layer  over  surface  (a)  while  the  peaks  penetrate  the 
laminar  sublayer  in  case  (b).  Then  (b)  should  produce  a  significantly 
greater  effect  on  the  flow  than  (a).  Surface  (c)  probably  would  not 
affect  the  flow  In  the  same  way  that  surface  (b)  would.  Clearly,  Ra  Is  an 
Inadequate  measure  of  "effective  roughness,"  k£.  Effective  roughness  la  a 
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Eqns.  (4)  and  (5)  are  of  the  form 

ke  -  f(Ra)  (6) 

Both  Schaffler  and  Koch  and  Smith  obtained  a  relation  for  effective 
roughness  as  a  function  of  only  Ra ,  as  shown  in  Eqn.  (6).  Relationships 
of  this  form  are  used  in  actual  practice  (and  will  be  used  in  the  model 
presented  below),  because  they  are  currently  the  best  interpretation  of 
available  data.  From  experiments  which  were  discussed  above,  however,  it 
is  known  that  effective  roughness  depends  on  a  number  of  geometric 
characteristics  of  the  surface,  that  is,  effective  roughness  can  be 
written  in  the  form  of  Eqn.  (7).  The  precise  form  of  Eqn.  (7)  is  not 
known  . 

ke  =  f(Ra,  peak  density,  peak  shape,  k  or  d-type,  etc.)  (7) 

A  series  of  systematic  experiments  would  be  required  to  obtain  a  better 
expression  for  effective  roughness  in  terms  of  measurable  surface 
characteristics  and,  as  suggested  by  Schaffler,  a  more  accurate 
specification  of  Eqn.  (7)  would  provide  a  significant  contribution  to  our 
understanding  of  roughness  effects  in  compressors. 

III.  A  Model  for  Calculating  the  Effect  of  Surface  Finish  on 
Axial  Compressor  Performance  in  an  Aircraft  Engine 
A.  Approach 

The  model  is  intended  to  provide  a  means  for  estimating  changes  in 
compressor  performance  which  result  from  changes  in  blade  or  vane  surface 
roughness.  The  flow  of  the  working  fluid  in  a  multistage  compressor  is 
exceedingly  complex,  and  thus  it  would  be  very  difficult  to  base  the  model 
on  the  details  of  the  interactions  between  individual  airfoils  and  the 
viscous  flow  passing  over  them.  It  is  possible,  however,  to  proceed  from 
the  assumption  that  each  stage  of  a  compressor,  or  the  compressor  as  a 
whole,  can  be  treated  as  a  control  system,  so  that  the  Integrated  effect 
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of  detailed  interactions  between  the  flow  and  the  airfoils  la  nodeled. 

That  approach  is  taken  here.  Viscous  losses  in  the  compressor  result  In  a 
temperature  rise  which  la  greater  than  that  for  lsentroplc  compression 
between  the  compressor  Inlet  and  exit  total  pressures.  Since  the  loss 
mechanism  stems  from  the  viscous  nature  of  the  flow,  we  examine  compressor 
efficiency  as  a  function  of  Reynolds  number  (the  ratio  of  Inertia  forces 
to  viscous  forces  In  the  flow)  with  airfoil  surface  finish  as  a  parameter. 
With  proper  definition  compressor  efficiency  can  be  made  Independent  of 
the  overall  compressor  ratio,  so  that  results  of  analyses  for  compressors 
with  different  pressure  ratios  can  be  directly  compared.  The  measure  of 
compressor  efficiency  that  allows  this  direct  comparison  la  polytropic 
efficiency  (this  term  Is  defined  In  Ref.  10).  Polytropic  efficiency  Is 
used  throughout  this  work.  The  discussion  which  follows  Is  based  largely 
on  the  work  of  Schaffler  (Ref.  3). 

B.  Efficiency  versus  Reynolds  Number 

Reynolds  number  effects  in  axial  flow  compressors  have  been 
studied  extensively.  Refs.  11,  12  and  3,  and  13  provide  examples  of  test 
results  of  cascade,  single  stage,  and  multistage  turbomachines, 
respectively.  I  have  drawn  several  conclusions  from  a  review  of  the  data 
in  these  references: 

First,  for  a  given  compressor  operating  point  (specified  by  pressure 
ratio  and  corrected  speed),  as  Reynolds  number  (based  on  first  rotor  blade 
mean  chord)  is  varied  from  low  to  high  values,  two  critical  values  are 

observed:  a  lower  critical  Reynolds  number.  Re  ,  and  an  upper  critical 

crit  L 

Reynolds  number,  Re  ^  y' 

Second,  below  Re  ,  compressor  efficiency  Increases  rapidly  with 

crit  L 

increasing  Reynolds  number. 
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Third,  above  Re  ,  compressor  efficiency  does  not  change  with 

crit  U 

Reynolds  number  and  is  a  function  of  the  effective  roughness  only. 

Fourth,  between  Re  ,  and  Re  the  relationship  between  C  and 

crit  L  crit  U 

Re  can  be  represented  by  the  equation 


where  K  - 
parameter  . 


1-  T)  =  K  Re 

proportionality  constant  and  n  is  a  Reynolds  number  sensitivity 
Figure  5  illustrates  conclusions  1  through  4. 


(8) 


Laminar 

Separation 

A. 


Turbulent  attached  Turbulent  attached 
Hydrodyn.  Smooth  Hydrodyn.  Rough 


R®crit  L 


R«crit  U  At  Re 


Figure  5.  Compressor  Performance  --  r\  versus  Re  (for  Given  PR  and  N /  6  ) 

The  poor  compressor  performance  in  the  region  below  Re  is 

cr  i  t  L 

attributed  to  laminar  separation.  In  this  region  the  laminar  boundary 
layer  on  a  compressor  blade  contains  insufficient  energy  to  overcome  the 
adverse  pressure  gradient  impressed  on  it  by  the  outer  flow,  and  it 
separates.  Compressor  flow  is  thus  reduced,  and  viscous  losses  increase 
sharply.  Additional  information  concerning  this  region  below  Re  crlt  Lmay 


he  found  in  Ref.  3  and  will  not  be  discussed  further  here. 
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In  the  region  he:  <■  n  Re  and  Re  ,  the  flow  is  turbulent  and 

crit  L  crit  U 

remains  attached.  When  all  surface  protuberances  are  contained  within  the 
laminar  sublayer  of  the  turbulent  boundary  layer,  the  flow  is  called 
"hydrodynamical ly  smooth."  Losses  are  independent  of  roughness  so  long  as 
all  roughness  elements  are  contained  in  the  laminar  sublayer.  With 
increasing  Reynolds  number  in  this  region,  the  turbulent  velocity  profile 
becomes  fuller  and  the  thickness  of  the  laminar  sublayer  is  reduced.  As 
roughness  elements  begin  to  penetrate  this  layer,  the  losses  become  a 
function  of  Reynolds  number  and  roughness. 

Above  Re  ,  the  flow  is  still  turbulent  and  attached,  but  most 

crit  U 

roughness  elements  penetrate  the  laminar  sublayer  and  losses  become  a 
function  of  the  roughness  only.  This  region  is  called  “ hy d r od y namic al ly 
rough."  Note  that  T|  increases  with  Reynolds  number  until  the  flow  becomes 
hyd r od ynamic al ly  rough.  Thus,  roughness  provides  a  constraint  on  best 
possible  compressor  efficiency  (at  a  given  compressor  operating  point). 

Schaffler  (Ref.  3)  has  experimentally  demonstrated  the  three  regions 
of  compressor  behavior  described  above  with  five-  and  six-stage 
compressors  (Figure  6,  taken  from  Ref.  3).  Also,  he  has  noted  that  the 
"curl  over"  region  (a  transition  region  between  hy dr odynamic al 1 y  smooth 
and  h yd r od y na mic a  1 1 y  rough  flow  in  which  C  is  a  function  of  both 
roughness  and  Reynolds  number)  occurs  over  a  small  Reynolds  number  range. 

The  effect  of  changing  airfoil  roughness  is  shown  schematically  in 
Figure  7.  Also  note  that  data  showing  this  same  effect  is  presented  in 
Figure  6.  Increasing  surface  roughness  reduces  the  efficiency  level  at 
which  the  flow  becomes  hydr odynamic al 1 y  rough.  Given  the  compressor 
operating  Reynolds  number,  Reop,  and  the  r\  versus  Re  curve  for  various 
surface  finishes,  it  is  possible  to  estimate  compressor  performance 
changes  associated  with  changes  in  surface  finish  and  the  surface  finish 
required  for  best  possible  performance  at  that  operating  point  (i.e. 


63 


U  SAFA-TR- 83-2 


h  yd  r  od  y  na  m  lc  a  1 1  y  smooth  flow).  Consider  the  situation  shown  In  Figure  7. 

For  the  given  operating  Reynolds  number,  surface  finishes  (In  terms  of  Ra) 

better  than  .508pm  (20  pin.),  would  provide  no  performance  Improvement. 

However,  If  surface  finish  Is  .762  pm  (30uln.)  and  can  be  reduced  to  .508pm 

(20pln.),  then  compressor  efficiency  should  Improve  by  two  percent  as  a 

result  of  that  change.  Note  that  the  curve  of  Tf  versus  Re  for 

hyd r odynamlc al 1 y  smooth  and  rough  flow  Is  defined  by  the  parameters  n  and 

Re  .  The  procedure  for  estimating  these  two  parameters  will  be 

crlt  U 

discussed  below. 

C.  Critical  Roughness  Reynolds  Number 

To  determine  Re  ,  the  point  of  transition  from 

crit  U 

hydrodynamlcal  ly  smooth  to  h  yd  r  ody  na  mlc  al  1  y  rough  flow  on  the  Tj  versus  Re 
curve  must  be  found.  This  transition  will  occur  when  roughness  elements 
on  the  surface  protrude  through  the  laminar  sublayer  of  the  turbulent 
boundary  layer  flow  over  the.  surfaces.  The  size  of  the  roughness  elements 
for  which  this  will  occur  has  been  related  to  the  equivalent  sand 
roughness  parameter  by  a  series  of  experiments  with  rough  flat  plates. 

These  experiments  have  led  to  the  definition  of  the  "critical  roughness 

Reynolds  number,"  Re  .  This  is  a  Reynolds  number  based  on  equivalent 

ks 

sand  roughness  that  defines  the  transition  from  smooth  to  rough  flow 

described  above.  Schllchting  (Ref.  7)  presents  the  following  simple 

formula  applicable  to  fully  turbulent  flow  over  a  rough  flat  plate: 

Rek„  *  ks  —  5  100  (9) 

s  v 

where  kg  Is  equivalent  sand  roughness.  W  j  Is  the  relative  velocity  of 
flow  past  the  plate,  and  v  Is  the  kinematic  viscosity.  Eqn.  (9)  means 

that  for  a  rough  flat  plate  with  Re k  less  than  100,  the  flow  Is 

s 

hydrodynamlcal ly  smooth.  As  suggested  by  Koch  and  Smith  (Ref.  9),  flow 
conditions  in  an  operating  compressor  may  be  significantly  different  from 
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those  of  fist  piste  flow  (l.e.  strong  adverse  pressure  gradient,  laminar 
and  turbulent  boundary  layer,  and  unsteady  flow)  and  therefore  Eqn.  (9) 
should  be  modified  for  use  in  coapressors.  They  recoaaend 


ks  <  90 

\> 


(10) 


where  k_  is  6.2  Ra .  Schaffler's  experiaental  results  (Ref.  3)  showed 


W, 

kg  _ L  <  88  (electrochemically  machined  blades) 


(11) 

(12) 


k_  —  <  1 35  (ground  blades) 

O  \) 

where  kg  is  8.9  Ra .  Not  that  Eqns.  (10)  through  (12)  can  be  written  in 
the  form 

ke  =  ke  <  E  (13) 


where  ke  is  the  effective  roughness,  Re/c  is  the  Reynolds  number  per  unit 
length  (based  on  stage  relative  velocity  and  blade  chord),  and  E  is  the 
empirically  determined  constant  as  in  the  right  hand  sides  of  Eqns.  (9)  to 
(12).  Considering  the  equal  sign  in  Eqn.  (13),  rearranging  and  taking 
logarithms  yields 

ln(kg)  =  ln(E)  -  lr.(Re/c)  (14) 

This  result  is  shown  graphically  in  Figure  8  for  a  general  (unspecified) 

value  of  E.  The  equation  is  a  straight  line  on  logarithmic  axes.  The 

region  below  the  diagonal  ^Eqn.  (14)J  ia  hyd  r  ody  namic  al  1  y  smooth  (  rj  “  f(n,Re)J 

f(n,Re)  and  that  above  the  diagonal  Is  hyd r odynamic al ly  rough  (  r)»  f(k  )] 

e 

The  bar  represents  the  distribution  of  stage  Reynolds  number  divided  by 
blade  mean  chord  length,  c,  (Re/c).  This  must  be  determined  for  the 
compressor  with  known  surface  finish  at  a  given  operating  point  (pressure 
ratio  and  corrected  speed)  and  aircraft  mission  point  (Mach  number,  M,  and 
al tltude  ,  h)  . 
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Increase-lnc.  Mach  number- M 
decrease-dec.  altitude -h 


ROUGH 


A=t(ke) 


dac.M 


dac.  ke 


dec.  h 


8MOOTH 
K  =f(n,Raop> 


In,  (Re/c) 

Figure  8.  Hy d r ody namic a  1 ly  Rough  and  Smooth  Regions 

The  location  of  the  Re/c  distribution  changes  (for  the  sane  compressor 
operating  point)  as  the  mission  changes;  high  altitude  and  low  Mach  number 
shift  the  distribution  to  the  left  and  low  altitude  and  high  Mach  number 
shift  It  to  the  right  In  Figure  8.  Changes  In  surface  finish  move  the 
Re/c  distribution  In  the  vertical  direction  in  Figure  8.  Best  possible 
compressor  performance  (at  a  given  operating  point)  occurs  when  the  Re/c 
distribution  Is  In  the  hyd r odyna mlc al 1 y  smooth  region.  It  is  not  known 
with  certainty  how  many  stages  of  a  multistage  compressor  must  be 
hydrodynamlcally  rough  before  the  whole  compressor  behaves  as  rough; 
Schaffler's  data  (Ref.  3)  Indicates  approximately  60  percent  of  the  total 
number  of  stages. 


Using  Eqn.  (13),  It  Is  possible  to  determine  Re  .  Note  that  when 

crit  U 

the  equal  sign  applies  In  Eqn.  (13), 
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where  c/fc e  Is  Che  ratio  of  blade  chord  to  effective  roughness  for  a  stage. 

So,  for  a  given  value  of  E,  and  with  a  knowledge  of  c/kgfor  each 

compressor  stage.  Re  Is  computed  as  follows: 

cri t  U 


Re  .  ,  *  —  l 

crit  U  Z  * 


where  Z  «  number  of  compressor  stages.  This  is  simply  the  average,  over 


all  machine  stages,  of  the  stage  Reynolds  number  corresponding  to  the 
transition  from  hydrodynamic al ly  smooth  to  rough  flow. 

D.  Reynolds  Number  Sensitivity  Parameter 

To  determine  the  Reynolds  number  sensitivity  parameter,  n,  recall 

that  below  Re  ,  T\  and  Re  are  related  by  Eqn  (8),  which  is 

crit  U 

1  -  T]  =  k  Re~n  (8) 

Hassell  (Ref.  14)  has  developed  an  empirically  based  correlation  for  the 

determination  of  n.  He  examined  V  versus  Re  for  20  different  compressors 
over  a  broad  range  of  Reynolds  number.  He  assumed  that 

n  =  p  •  q  (17) 

where  p  depends  on  the  mean  Mach  number  level  in  the  compressor  and  q  is  a 
function  of  compressor  geometry.  Schaffler  (Ref.  3)  compared  experimental 
values  of  n  with  those  found  from  the  Hassell  correlation  and  found  good 
agreement.  Refer  to  Ref.  13  for  details  of  the  use  of  the  Hassell 
correlat ion . 

E.  Calculation  of  Efficiency  Change 

If  the  value  of  n  for  the  compressor  and  values  of  Re  for 

crit  U 

various  surface  roughnesses  are  known  (corresponding  to  a  given  compressor 
operating  point),  then  curves  of  tj  versus  Re  can  be  generated  for  the 
various  roughnesses.  For  hydrodynamlcally  smooth  flow,  we  have  Eqn.  (8): 

1-T?  -  K  Re~n  (8) 
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So,  with  respect  to  reference  values,  *1  f  and  Re  ,  we  may  write 


ref  "  -ref' 
-n 


_LJ_ 

0-n) 


ref 


If  we  take  logarithms  and  then  exponentials  of  each  aide  we  obtain 

1  -T?  =  U-^ref  exP 

Then,  defining  AT)  as 


AT?  =  T?ref  -TJ 


(18) 


(19) 


(20) 


we  obtain 


AT? 


(1-T?) 


ref 


- 

exp 

; &  Jj 

-  i 

(21) 


Eqn  (21)  can  be  used  to  compute  the  compressor  efficiency  change  between 
two  values  of  Re  ^ cor r e spond ing  to  two  different  roughnesses*  Recall 

that  Re  is  computed  by  Eqn.  (16). 

c  r  i  t  U 


F.  Relationship  between  T)  versuB  Re  and  K  versus  Re/c 
Curves . 

Figure  9  shows  how  the  results  of  the  preceding  sections  can  be 
used  to  estimate  the  effects  of  surface  roughness  on  compressor 
e  f  f ic ienc  y . 


Figure  9.  Relationship  between  k  versus  Re/c  and  T)  versus  Re 
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In  part  (a)  of  Figure  9,  Re/c  distributions  corresponding  to  two  operating 
Reynolds  numbers  (Re  j  and  Re2)  and  two  roughnesses  (ke  ^  and  ke^)  are 
shown.  The  bar  represents  the  distribution  of  Re/c  for  the  compressor. 
Corresponding  operating  Reynolds  numbers  and  roughnesses  are  shown  In  part 
(b).  Consider  a  compressor  operating  at  point  B  (roughness  - 
Reducing  the  Re/c  distribution  to  point  A  (by  Increasing  altitude,  for 
example)  produces  no  change  In  T)  because  the  Re/c  distribution  stays  In 
the  hyd r odynamic al ly  rough  region.  Reducing  roughness  from  B  to  D 
produces  a  two  percent  Increase  In  r?  at  the  same  operating  Reynolds  number 
even  though  the  flow  is  still  hyd r odynamic al ly  rough  at  D.  Now  consider 
moving  from  D  to  C  by  reducing  the  operating  Reynolds  number.  As  shown  In 
part  (a),  the  Re/c  distribution  soon  enters  the  smooth  region  and  then  T? 
decreases  according  to  1-  I)  ■  k  Re  "  as  shown  in  part  (b). 


G.  Model  Summary 

Eqns.  (7),  (13),  (16),  (17),  and  (21)  constitute  the  empirical 
model  for  estimating  roughness  effects  in  axial  flow  compressors.  They 
are  repeated  here  for  convenience. 


k  «  f(Ra,  peak  density,  peak  shape,  k-  or  d-type ,  etc.) 


k  *£  1  E 
e  c 


Re 


crlt 


(7) 

(13) 

(16) 


n  »  p  •  q  where  p  and  q  are  determined  by  the  method  of  Ref.  14  (17) 


AT?  -  (l-T?) 


ref 


- 

exp  j 

;•  -  fcJ: 

-  i 

(21) 
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Until  further  date  becoaes  available,  I  recommend  the  use  of  Schaffler's 
reaulta  In  Eqns.  (7)  and  (13).  Thus,  Eqns.  (4)  and  (11)  are  uaed  In  the 

Illustrative  examples  In  the  next  section  so  that 

kg  -  8.9  Ra  (22) 

k  -  Re  <  88  (23) 

e  c 

I V .  Application  of  the  Model 

A.  Choice  of  the  Application 

In  this  section,  results  are  presented  from  the  application  of  the 
empirical  model  to  the  compression  sections  of  three  Air  Force  Engines: 
the  TF33-PW-7,  the  P100-PW-100,  and  the  J85-GE-5H.  The  surface  finishes 
required  for  hydrodynamlcally  smooth  flow  at  various  mission  points  are 
estimated.  Also,  the  changes  In  compressor  efficiency,  which  could  be 
achieved  by  reducing  airfoil  roughness  to  that  for  hydrodynamlcally  smooth 
flow,  are  computed. 

The  TF33-PW-7  was  chosen  for  three  reasons:  first,  It  Is  a  primary 
Air  Force  fuel  user,  second,  other  studies  of  airfoil  surface  finish  have 
been  performed  on  this  or  similar  engines  (Refs.  1  and  2),  and  third.  It 
Is  a  transport  engine  with  a  well-defined  mission  profile.  The 
F100-PV-100  was  chosen  because  It  Is  a  fighter  engine  representing 
relatively  current  technology.  The  J85-GE-5H  was  chosen  because  of  Its 
potential  use  In  a  roughness  effects  test  In  the  Aero  Propulsion 
Laboratory's  Compressor  Research  Facility.  A  J85  test  rig  has  been 
procured  for  facility  checkout  and  It  will  be  used  as  a  vehicle  for 
testing  roughness  effects. 


B.  TF33-PW-7 


The  TF33-PW-7,  7A  Is  used  to  power  the  USAF  C-141  fleet.  The 
engines  were  first  delivered  In  1964  and  are  currently  consuming  about  600 


j 
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million  gallons  of  fuel  per  year  (Ref.  15).  Table  1  shows  an  aggregate 
mission  for  the  C141B;  It  Is  a  statistical  representation  of  23  different 
missions  Including  training,  exercise,  airlift,  and  contingency  flight 
profiles  (Ref.  15). 


TABLE  1 

AGGREGATE  MISSION  FOR  C141B 

Altitude 

Feet 

M 

Power 

Time 

Minutes 

Takeoff 

Sea  Level 

0.15 

Max.  T/O 

6.0 

Climb 

(1) 

<1> 

NR 

16.0 

Cruise 

37,000 

0.74 

75-80%  NR 

203.0 

Descent 

(2) 

(2) 

Idle 

21.0 

Approach  and  Land 

500* 

30%  MR 

18.0 

Taxi  (In  and  Out) 

0 

Hi 

(die 

25.0 

— 

TOTAL 

289.0 

(1)  250  Kts  to  10,000  Ft.,  280  Kts  to  0.7  M 


Fuel  Increased  50%  to  Reflect  Touch  &  GO  &  GO  Around  Training 
(2)  Reverse  Climb  Path 

MRs  Military  Rated 
NRr  Normal  Rated 


Note  that  the  aircraft  Is  In  the  cruise  condition  for  70  percent  of  total 
alsalon  time.  The  empirical  model  for  roughness  effects  can  be  applied 
for  any  given  aircraft  mission  point  (M  and  altitude)  and  compressor 
operating  point.  Only  the  cruise  condition  Is  considered  here  since  It  is 
the  major  portion  of  the  aggregate  mission. 

Figure  10  shows  the  distribution  of  Reynolds  number  per  unit  length 
(Re/c)  for  the  compression  section  of  this  engine  at  M  -  0.75  and  h  - 
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35,000  feet,  with  an  assumed  Ra  of  1.778  ui  (70  p  In).  This  la  a 
conservatively  high  Ra  froa  field  measurements . 


(//In) 

Zk  Ra  (fjm)  4a 

70-44  +  1.46* 

1.778  1.118 

70-20  +1.62* 

1.778  -*  .608 

70  +  13  +1.87* 

1. 778-+. 330 


Figure  10.  Re/c  Distribution  for  TF33-PW-7  —  Cruise  Mission 


Note  that  the  fan,  low  pressure  compressor  (LPC),  and  high  pressure 
compressor  (HPC)  are  all  In  the  hydr ody namlc al 1 y  rough  region.  Also,  note 
the  values  of  n,  which  were  computed  using  the  Hassell  correlation  for 
each  component.  These  components  would  have  to  be  smoothed  to  Ra  -  1.118, 
.508,  and  .330  ua,  respectively.  In  order  for  the  flow  to  be  completely  In 
the  hydrodynaalc al ly  saooth  region.  Component  efficiency  changes 
corresponding  to  a  reduction  froa  Ra  “  1.778pm  to  these  levels  are  1.45 
percent,  1.62  percent,  and  1.87  percent.  An  Aero  Propulsion  Laboratory 
analysis  of  the  TF33-PW-7  for  similar  efficiency  changes  yielded  an 
approximate  1.5  percent  reduction  In  TSFC  at  the  cruise  thrust  level. 

This  result  Is  consistent  with  the  experimental  data  of  Refs.  1  and  2. 
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C.  F 1 OO-PW-1 00 

The  F100-PW-100  Is  used  in  the  F-1S  aircraft.  In  contrast  to  the 
TF33,  the  FIOO  is  required  to  operate  over  an  extensive  range  of  flight 
conditions.  The  aggregate  mission  for  the  F100  in  the  F-15  has  more 
elements  and  is  more  complex  than  that  for  the  TF33.  Figure  11 
illustrates  the  effect  of  Mach  number  and  altitude  on  the  distribution  of 
Reynolds  number  per  unit  length  for  the  compression  section  of  this 
engine  . 


(Re/c) 


Figure  11.  Re/c  Distribution  for  F200-PW-100  —  Various  Missions 

The  M  "  .88,  h  “  45,000  feet  condition  represents  a  ferry  mission,  while 
the  low  altitude  mission  demonstrates  the  high  values  of  Re/c  (and 
corresponding  low  values  of  Ra  required  for  hydrodynamically  smooth  flow) 
which  can  occur  during  F100  operation.  Given  the  amount  of  time  spent  at 
each  altitude/Mach  number  condition  in  the  aggregate  mission,  the 
corresponding  Re  per  unit  length,  and  the  model  described  in  this  paper, 
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it  is  possible  to  estimate  an  aggregate  TSFC  reduction  which  would  result 
from  changes  in  compressor  airfoil  surface  finish.  No  further  analysis  or 
discussion  will  be  presented  here,  since  the  primary  purpose  of  this 
example  is  to  illustrate  Mach  number  and  altitude  effects  on  compressor 
Re/c  distributions. 

D.  J  8  5-GE-5H 

A  J85-GE-5H  compressor  is  being  used  to  verify  the  operation  of 
the  Aero  Propulsion  Laboratory's  Compressor  Research  Facility.  This 
compressor  will  be  used  for  an  investigation  of  roughness  effects  on  J-85 
compressor  efficiency.  The  Compressor  Research  Facility  provides 
compressor  inlet  pressure  control  (atmospheric  to  approximately  2  psia), 
so  performance  for  a  range  of  Reynolds  numbers  can  be  examined. 

Production  J85  blades  and  vanes  have  a  surface  roughness  of  approximately 
Ra  •  .508  urn.  A  parametric  variation  of  surface  finish  will  be  used  in 
this  test  in  an  attempt  to  isolate  and  measure  roughness  effects  on 
efficiency  independent  of  other  influences.  Figure  12  shows  the  Re/c 
distributions  for  a  J-85  compressor  (at  design  conditions)  over  a  range  of 
inlet  pressures  and  surface  roughnesses.  Figure  13  shows  curves  of 
efficiency  versus  Reynolds  number  corresponding  to  the  roughness  levels  in 
Figure  12.  Note  that  the  range  of  surface  roughnesses  and  inlet  pressures 
suggested  in  Figure  13  provides  ample  opportunity  for  investigation  of  the 
transition  from  hyd r ody nam ic a  1 1 y  smooth  to  rough  flow  in  the  J85 
compressor.  This  example  has  been  provided  to  illustrate  the  application 
of  the  empirical  model  in  forecasting  test  results  for  a  planned 
Compressor  Research  Facility  test. 

V .  Conclusion 

This  paper  has  provided  an  xplanation  of  some  of  the  problems 


75 


USAFA-TR-83-2 


associated  with  surface  roughness  measurement  and  a  description  of  an 
empirical  model  which  can  be  used  to  estimate  the  effect  of  airfoil 
surface  finish  changes  on  axial  compressor  performance.  It  must  be 
emphasized  that  compressor  efficiency  changes,  which  are  calculated 
through  the  use  of  this  model,  must  be  treated  as  estimates;  m^ch  more 
experimental  data  is  required  for  "calibration"  of  the  model.  However, 
this  model  is  useful  in  that  it  provides  insight  into  the  problems  of 
estimating  airfoil  surface  finish  effects  on  compressor  efficiency.  The 
relationship  between  measured  surface  properties  and  effective  roughness 
is  of  primary  importance  and  the  understanding  of  that  relationship 
requires  additional  investigation. 


Symbols 


English  Symbols 

b  roughness  element  width  in  k-  and  d-type  roughness 

c  blade/vane  hord  length 

d  roughness  element  width 

D  peak  spacing 

E  empirical  constant  In  equation  (13) 

EGT  exhaust  gas  temperature 

h  altitude 

k  roughness  element  height 

kg  effective  roughness 

ks  equivalent  sand  roughness 

ks  roughness  defined  by  Eqn.  (3) 

K  proportionality  constant 

L  sampling  length  in  Ra  definition 

M  Mach  number 

n  Reynolds  number  sensitivity  parameter 
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N 

P 

PR 

q 

R 

Ra 

Re 


Re 

Re 


k 

s 

op 


Re 


c  r  i  t 


L 


Recrit  U 


TSFC 

“l 

W 

y 

y  y 

peak  groove 

z 


engine  speed 

parameter  used  to  define  n;  n  «  p  •  q 
pressure  ratio 

parameter  used  to  define  n;  n  «  p  •  q 
pipe  radius 

arithmetic  average  roughness 

Reynolds  number  (usually  based  on  midspan  of  first 
rotor  ) 

critical  roughness  Reynolds  number 

Re  corresponding  to  the  compressor  operating  point 
and  aircraft  mission  conditions 

lower  critical  Reynolds  number  based  on  first 
rotor  blade  mean  chord 

upper  critical  Reynolds  number  based  on  first 

rotor  blade  mean  chord 

thrust  specific  fuel  consumption 

roughness  element  spacing  in  k-  and  d-type 

r  oughne  s  s 

relative  velocity 

ordinate  or  curve  of  roughness  profile 
see  definition  of  k  ,  Eqn.  (3) 
number  of  stages  in  turbomachine 


Greek  Symbols 


V 

u 

e 

% 


polytropic  efficiency 
kinematic  viscosity 

temperature  ratio  -  used  in  corrected  speed  N 
flow  loss  coefficient  (pipe  flow) 
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Subscripts 

1  denotes  1th  stage  of  t ur bomachine 

ref  reference  values  of  TJ  and  Re 
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AN  EXAMPLE  OF  UNCERTAINTY  ANALYSIS  IN  COMPRESSOR  TESTING 
R.J.  Stiles*  and  C.A.  Boedicker** 

Abstrac  t 

This  paper  presents  an  overview  of  the  most  fundamental  and  important 
concepts  in  uncertainty  analysis,  a  technique  which  can  be  used  by  a  test 
engineer  to  estimate  the  accuracy  of  measured  data  and  values  that  are 
calculated  from  measured  data.  Some  of  the  results  of  an  uncertainty 
analysis  which  was  accomplished  for  a  turbine  engine  compressor  test  are 
used  to  illustrate  the  analysis  procedure  and  results.  Also,  the  benefits 
of  doing  an  uncertainty  analysis  as  part  of  a  test  are  briefly  discussed. 

I .  Introduct ion 


Performance  indices  and  aer othermodynamic  parameters  that  are  used  to 
assess  the  performance  of  turbine  engine  components  are  based  upon 
measurements  of  basic  quantities  like  pressure,  temperature,  force,  and 
position.  Values  for  these  basic  quantities  are  obtained  by  processes 
that  inevitably  yield  measurement  error  (i.e.,  the  difference  between  the 
true  value  of  the  quantity  of  Interest  and  the  value  indicated  by  a 
particular  measurement  process).  Furthermore,  if  a  performance  index  is 
computed  as  a  function  of  measured  values,  then  the  effect  of  measurement 
error  is  to  introduce  error  in  the  performance  index  itself;  the 
measurement  error  propagates  through  the  performance  calculation. 

Uncertainty  analysis  is  a  procedure  that  can  be  used  before  a  test  to 
estimate  the  maximum  measurement  errors  that  can  reasonably  be  expected 
for  particular  measurement  processes  and  the  effects  of  these  errors  on 
calculated  parameters,  which  are  derived  from  the  measured  values.  In 
order  to  conduct  a  meaningful  evaluation  of  a  turbine  engine  component,  it 
Is  necessary  to  understand,  and  often  to  minimize,  the  sources  of  error  in 
both  measured  and  computed  values. 

Our  objective  in  writing  this  paper  is  to  provide  an  overview  of  the 
fundamentals  of  uncertainty  analysis  and  to  present  some  results  from  an 

♦Major,  USAF,  AFWAL  Aerc  Propulsion  Laboratory 
**C1C,  USAF  Academy 
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uncertainty  analysis  which  was  performed  for  the  test  of  a  centrifugal 
compressor  at  the  AFWAL  (Air  Force  Wright  Aeronautical  Laboratories)  Aero 
Propulsion  Laboratory's  Compressor  Research  Facility  (CRF), 

Wright-Patter son  AFB,  Ohio.  Ref.  1  presents  a  complete  description  of  the 
CRF.  We  have  not  written  a  step-by-step  procedure  for  accomplishing  an 
uncertainty  analysis;  the  interested  reader  should  see  Ref.  2  for  that 
kind  of  information.  This  overview  of  r  specific  application  of 
uncertainty  analysis  and  the  accompanying  examples  are  intended  to 
illustrate  the  utility  of  the  analysis  procedure  and  results. 


II  Compressor  Flowpath  and  Performance  Calculations 

The  Garrett  Turbine  Engine  Company  has  designed  and  manufactured  a 
variable-geometry  centrifugal  compressor  (Figure  1),  which  will  be  tested 


at  the  CRF. 


Figure  1.  Variable-Geometry  Centrifugal  Compressor 
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At  its  design  operating  speed  (23,000  RPM)  ,  this  compressor  produces  a 
pressure  ratio  of  8:1  in  a  single  stage  and  has  an  airflow  rate  of 
approximately  11.34  Kg/ sec  (25  lbm/sec).  The  test  Is  designed  to 
demonstrate  the  effects  of  variable  inlet  guide  vanes  and  variable 
diffuser  area  on  centrifugal  compressor  performance.  Analysis  has  shown 
that  near-deslgn  levels  of  compressor  ratio  and  efficiency  can  be 
maintained  over  a  broad  range  of  flow  rate  by  properly  adjusting  the 
variable  geometry  hardware.  The  ability  to  modulate  airflow  In  a  turbine 
engine  by  simple  geometry  changes  Is  quite  useful.  For  example.  In  a 
turbofan  application  the  v ar 1 ab 1 e- ge ome t r y  compressor,  when  coupled  with  a 
variable-geometry  turbine,  would  provide  variable  engine  bypass  ratio  by 
changing  core  airflow.  Bypass  ratio  could  then  be  adjusted  to  minimize 
fuel  consumption  for  a  range  of  flight  conditions. 

A  cross-section  of  the  compressor  flowpath,  with  measurement  station 
designation,  is  shown  In  Figure  2.  2.9  3.0 


DESWIRL  VANE 


Figure  2.  Compressor  Flowpath  and  Station  Designation 
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More  than  300  Individual  measurement  probes  have  been  Installed  In  this 

compressor  and  20  different  measured  values  M  (J«I,20)  along  the  flowpath 

j 

are  used  in  the  calculation  of  65  performance  and  aer othermodynamic 
parameters  c  (i«l,65).  Measured  values  of  pressures  and  temperatures  at 
a  specific  measurement  station  along  the  flowpath  are  derived  from  an 
arithmetic  average  of  the  measurements  from  individual  probes 
circumferentially  distributed  (with  respect  to  the  flow)  at  that  station. 

A  detailed  description  of  all  of  the  measurements  and  calculations 
associated  with  this  test  la  outside  the  scope  of  this  paper.  Discussions 
of  specific  measurements  and  calculations  are  presented  where  appropriate 
in  the  sections  that  follow. 

Ill .  Uncertainty  Analysis  Overview 

In  this  section  we  overview  the  most  fundamental  and  important 
concepts  in  uncertainty  analysis.  The  Interested  reader  should  refer  to 
the  many  books  and  journal  articles  concerning  uncertainty  analysis  for 
additional  details  (Refs.  2  through  5,  for  example).  We  have  followed  the 
methodology  of  Refs.  2  and  5  in  the  Information  presented  here. 

"Measurement  error"  is  the  difference  between  a  measured  value  and 
the  true  value,  as  defined  by  the  National  Bureau  of  Standards  (NBS),  of 
the  quantity  being  measured  (Ref.  2).  "Measurement  uncertainty"  is  the 
maximum  measurement  error  that  might  reasonably  be  expected;  uncertainty 
is  a  measure  of  accuracy.  The  problem  in  uncertainty  analysis  is  to 
construct  an  uncertainty  Interval  which  limits  measurement  error. 

Assuming  that  there  are  no  gross  errors  in  the  measurement  (l.e., 
human  error  or  instrument  malfunctions),  then  there  are  two  components  of 
measurement  error:  (1)  precision  (random)  error  and  (2)  bias  (fixed) 
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A.  Precision  Error 

If  a  particular  measurement  process  is  repeated  many  times,  one 
will  find  that  the  results  can  be  plotted  In  a  frequency  distribution,  as 
illustrated  in  Figure  3. 


Figure  3.  Precision  Error 


The  frequency  distribution  will  usually  be  a  normal  one,  the  consequence 
of  random  errors  from  a  number  of  Independent  Influences  in  the 
measurement  system.  Precision  error  is  the  variation  between  repeated 
measurements;  the  standard  deviation  (o)  of  the  frequency  distribution 
is  used  as  a  measure  of  this  error.  Also,  for  finite  sample  size,  the 
statistic  S  is  used  to  calculate  an  estimate  of  the  standard  deviation.  S 
is  called  the  precision  index  and  is  calculated  as  follows. 


S 


N 

L  (X  -  X)2 
1-1 


N-l 


(1) 


where  N  -  number  of  messurements  represented  In  the  frequency  distribution 
and  X  -  average  value  of  the  individual  measurements  X^ . 
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It  is  particularly  significant  to  note  that  the  precision  error  can 
be  reduced  by  making  many  measurements  of  the  same  variable.  The  standard 
deviation  of  the  mean  of  "n"  independent  measurements  is  equal  to  1 //n~ 
times  the  standard  deviation  of  the  frequency  distribution  for  a  single 
measurement.  Therefore,  the  use  of  redundant  instrumentation  and/or 
multiple  measurements  over  a  period  of  time  for  values  which  must  be  known 
accurately  is  very  important  in  turbine  engine  testing.  It  is  important 
to  distinguish  between  N  and  n.  N  corresponds  to  the  number  of  data  scans 
one  uses  to  arrive  at  a  single  measured  value  of  the  parameter  of 
interest.  The  symbol  n  refers  to  either  the  number  of  times  the 
measurement  is  obtained  over  a  period  of  time  or  the  number  of  redundant 
measurements  of  the  same  parameter. 

B  .  Bias  Error 

The  second  component  of  measurement  error,  bias  error,  is 
significantly  more  difficult  to  estimate  than  precision  error.  Bias  error 
(sometimes  called  systematic  error)  is  constant  during  a  test  and, 
therefore,  averaging  repeated  measurements  does  not  affect  the  bias.  Bias 
can  only  be  determined  by  comparing  a  measured  value  with  the  true  (NBS) 
value  of  the  quantity  being  measured.  Known  bias  is  normally  determined 
by  a  calibration  procedure  in  which  a  measurement  instrument  is  compared 
(and  corrected)  so  that  it  agrees  with  a  standard  instrument  over  the 
measurement  range.  However,  after  elimination  of  the  known  bias,  unknown 
bias  remains,  and  it  must  be  estimated  as  part  of  the  uncertainty 
analysis.  In  most  cases,  the  data  engineer  must  use  his  judgment  to 
estimate  an  upper  limit  on  the  unknown  bias  error;  i.e.,  he  must  specify 
the  "bias  limit,  B . "  Note  that  the  estimate  of  bias  error  does  not  have 
a  statistical  basis  as  does  the  estimate  of  precision  error.  Figure  4 
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Illustrates  the  bias  and  precision  error  that  combine  too  form  measurement 


Figure  4.  Bias  and  Precision  Error 
C.  Combining  Elemental  Error  Sources 

As  part  of  an  uncertainty  analysis,  one  must  develop  estimates  of 
bias  and  precision  errors  for  each  measured  quantity.  Before  this  can  be 
done,  though,  the  measurement  process  must  be  defined.  Specifically,  the 
test  procedure.  Instrumentation  type,  and  number  of  sensors  must  be  known. 
Then  estimates  of  bias  and  precision  errors  for  measurement  subprocesses 
(l.e.,  calibration,  data  acquisition,  and  data  recoding)  can  be  made.  The 
elemental  estimates  for  the  subprocesses  are  then  combined  by 
r oo t- sum- sq ua r e  (as  In  Eqns.  (2)  and  (3)  to  give  bias  limit  and  precision 
index  estimates  for  individual  measurements. 
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D.  Combining  Bias  and  Precision  --  Measurement  Uncertainty 

It  is  convenient  to  combine  the  estimates  of  bias  and  precision 
error  into  a  single  number  (i.e.,  measurement  uncertainty),  which  can  be 
interpreted  simply  as  the  largest  error  one  should  reasonably  expect  for  a 
given  measurement.  A  more  rigorous  statement  of  the  meaning  of 
uncertainty  is  not  possible  because  it  is  the  combination  of  a  statistic 
(precision)  and  a  judgment  (bias).  At  best,  we  can  say  that  the 
measurement  uncertainty  estimate  provides  a  number  that  is  analogous  to  a 
confidence  interval.  The  author  of  Ref.  5  suggests  calculation  of 
measurement  uncertainty  by  the  equation 

%  =  +  (B  +  t95  (4) 

where  B  is  the  bias  limit,  S  is  the  precision  index,  and  t^,.  is  the  95th 
percentile  point  for  the  two-tailed  Students  t -d i s t r lbut i on .  The  tg5 
value  depends  on  sample  size;  it  decreases  with  increasing  sample  size  and 
has  a  lower  limit  of  1.96.  For  sample  sizes  of  30  or  more,  t^5  “2  is  often 
used  in  Eqn.  (4).  Note  that  the  effect  of  t  is  to  inflate  and  thus 
reduce  the  risk  of  underestimating  S  when  a  small  sample  is  used  to 
calculate  S.  The  subscript  99  on  U  indicates  (loosely)  a  99  percent 
confidence  level  for  the  uncertainty  computed  from  Eqn.  (4).  Figure  5 
illustrates  the  concept  of  measurement  uncertainty. 

E.  Uncertainty  of  Calculated  Values 

Measurement  uncertainties,  which  are  computed  using  Eqns.  (2) 
through  (4)  propagate  through  calculations  involving  those  measurements. 

One  can  estimate  the  effect  of  measurement  uncertainty  on  a  given 
calculated  value  by  expanding  the  calculated  value  (which  is  a  function  of 
j  measurements)  in  a  Taylor  Series.  It  can  be  shown  that  the  precision 
and  bias  errors  of  a  computed  value  are  (for  a  first-order  Taylor  Series 
expansion) 
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(5 


(6 


where  the  subscripts  on  S  and  B  pertain  to  the  i  computed  value  (C  )  and 

th  1 


,  th 


the  j  measured  value  ( Mj ) .  Details  for  the  derivation  of  Eqns.  (5)  and 


(6)  are  given  In  Appendix  B  of  Ref.  2.  The  uncertainty  of  calculated 
values  can  be  determined  using  Eqn.  (4). 

Note  that  In  order  to  determine  the  bias  limit  and  precision  Index 
(and  then  the  uncertainty)  of  a  computed  value,  the  partial  derivatives 

in  Eqns.  (5)  and  (6)  must  be  known.  For  a  few  simple  calculations 
one  may  be  able  to  compute  the  partial  derivatives  analytically.  However, 
In  moat  experiments  In  which  many  calculations  and  measurements  are 
Involved,  It  Is  more  convenient  to  use  a  numerical  procedure  to  estimate 
these  terms.  For  the  test  that  Is  the  subject  of  this  paper,  there  were 
65  calculations  Involving  20  Independent  measured  values.  The  partlals 


3C  /3M.  (sometimes  called  the  in  '  .nee  coefficient  matrix)  were  computed 


1  j 

using  the  first  order  approximation 


aci  c;  -  c* 


'  Mj  Mj  '  Mj 


(7 


where  is  derived  from  the  nominal  set  of  measurements  and  c|  Is 
derived  from  a  measurement  set  with  all  nominal  values  except  Mj  ,  which 
is  and  Is  one  percent  higher  than  its  nominal  value.  Note  that  the 

terms  given  by  Eqn.  (7)  will.  In  general,  have  dimensions.  It  Is 
desirable  to  be  able  to  review  the  influence  coefficient  matrix  to 
determine  the  most  significant  terms,  but  this  is  not  possible  unless  all 
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the  terns  have  the  same  dimensions  or  are  dimensionless.  Therefore,  we 
chose  to  alter  the  influence  coefficient  matrix  as  follows. 

3Ci  M.  (Cj  -  C^/C 

Dimensionless  Influence  Coefficient  =  ^  =  -^7 — ^  (8) 

j  i  j  i’  i 

Since  M'.  is  always  one  percent  greater  than  M.  ,  influence  coefficients 
in  the  form  of  Eqn.  (8)  give  the  percentage  change  in  C  corresponding  to 
a  one  percent  change  in  M  •  Also,  influence  coefficients  in  dimensionless 
form  can  be  used  in  Eqns.  (5)  and  (6)  as  long  as  bias  and  precision  errors 
are  expressed  as  percentages. 

IV .  Results  and  Illustrative  Examples 

We  have  chosen  to  present  the  results  of  the  centrifugal  compressor 
uncertainty  analysis  for  only  10  of  the  65  calculated  parameters  involving 
10  (of  the  total  20)  different  measurements,  since  our  objective  is  only 
to  illustrate  the  analysis  procedure  and  results.  Table  X  lists  these  10 
measurements  and  also  gives  estimated  bias  limits  and  precision  indices 
(percent  of  reading)  for  the  measurements. 

Figure  6  is  the  dimensionless  influence  coefficient  matrix  for  the  10 
measurements  and  10  calculations  we  have  chosen  to  highlight.  Values  in 
each  row  are  the  percentage  change  in  the  calculated  value  (at  the  far 
left  of  each  row)  corresponding  to  a  one  percent  change  in  the  measurement 
directly  above  the  value.  One  can  compare  magnitudes  of  the  values  in 
each  row  to  find  the  most  important  measurements  affecting  a  calculation. 
Comparing  magnitudes  of  the  values  in  each  column  shows  w,,ich  calculated 
parameters  are  most  affected  by  a  given  measurement.  The  influence 
coefficient  matrix  is  a  valuable  tool,  not  only  for  determining  Important 
measurements,  but  also  for  direct  estimates  of  the  effects  of  measurement 
changes.  It  is  Important  to  note  that  the  influence  coefficients  in 
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Figure  6  are  valid  for  only  one  aet  of  aeasured  variable*.  In  thia  caee 
they  correspond  to  condition*  at  100  percent  coapreaeor  dealgn  apeed. 

Table  II  liata  the  10  calculated  paraaetera  along  with  the  expected 
uncertainties  in  these  values  (expressed  as  a  percentage).  The 
uncertainties  were  derived  using  equations  (4)  through  (6),  the 
aeasureaent  errors  shown  in  Table  It  and  the  Influence  coefficients  in 
Figure  6.  These  values  are  only  valid  at  100  percent  design  speed;  at 
other  speeds  we  would  expect  changes  in  the  influence  coefficients  and 
thus  the  aeasureaent  uncertainty  values. 

He  now  present  two  examples  that  serve  to  clarify  and  illustrate  the 
use  of  the  results  shown  in  Tables  I  and  II  and  Figure  6. 

A.  Example  1:  Mach  Nuaber  Calculat Iona 

There  are  four  Mach  number  calculations  represented  in  Figure  6 
(not  including  the  relative  Mach  nuaber  calculation  RM20);  they  are  AM19, 
AM20 ,  AM25,  and  AM29.  Referring  to  the  aeasureaent  station  in  Figure  2, 
note  that  these  are  the  Mach  nuaber  values  before  the  inlet  guide  vanes, 
before  the  Impeller,  before  the  diffuser,  and  after  the  diffuser/90  degree 
bend,  respectively.  The  calculated  Mach  nuaber s  are  of  Interest  both  in 
theaaelves,  to  indicate  coapressibllity  of  the  flow,  and  as  eleaents  used 
in  additional  calculations.  For  exaaple,  the  total  pressure  loss  across 
the  inlet  guide  venes  (IGVs)  depends  on  IGV  angle  and  (to  a  lesser  extent) 
AM19.  The  calculations  of  AM20  and  AM25  are  neceesary  for  the 
deteralnatlon  of  flow  velocity  coaponents  at  atatlons  2.0  and  2.5, 
reepectlvely .  AM29  is  used  to  celculate  static  teaperature  at  station 
2.9,  which  is  used,  in  turn,  to  calculate  diffuser  efficiency. 

The  calculation  of  these  Mach  nuabers  provides  a  useful  exaaple  of 
the  Influence  of  calculation  procedure  on  the  uncertainty  of  calculated 
velues  AM19,  AM20,  and  AM25  are  all  calculated  using  aeasured  or  coaputed 
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total  prasaura  and  total  temperature  at  tha  station  of  intaraat  along  with 
tha  continuity  aquation.  Using  thesa  total  propartias  and  tha  aquation 


(Ref.  6) 


(9) 


where  W  la  airflow,  t  la  total  temperature,  la  total  prasaura,  Cg  la 
discharge  coefficient,  A  la  area,  «  is  flow  swirl  angle,  Y  la  tha  ratio  of 
specific  heats  of  tha  flow,  and  M  la  Mach  number,  one  can  solve  for  M  by 
an  Iterative  procedure  If  «  la  known.  Note  that  at  station  2.5,  "  la  not 
known,  so  that  Euler* a  Turbine  aquation  la  used  together  with  Eqn.  (9)  to 
solve  iteratively  for  AM25  and  «  .  (Euler's  Turbine  equation  is  obtained 
by  applying  conservation  of  energy  and  the  aoaentua  equation  to  a  stream 
tube  past  a  rotor  blade.  For  wore  information  see  Ref.  7,  pages  10  to  13 
or  pages  17  to  22.) 

AM29  la  calculated  from  measured  static  and  total  pressure  values  at 
station  2.9  using  the  lsentroplc  relationship  (Ref.  6) 


H^»2) 


where  P  la  static  pressure.  Solving  for  M  yields 

b 


(11) 


He  note  from  Figure  6  thet  the  Influence  coefficients  for  AM29 
corresponding  to  the  measurements  FS29  and  PT29  are  the  largest 
coefficients  In  the  matrix.  It  la  Interesting  to  examine  the  partial 
derivatives  of  Eqn.  (11)  with  raspect  to  P  and  P  .  Expressions  for 

5  T 

these  pertlal  derivatives  are 
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Dimensionless  Influence  coefficients  corresponding  to  Eqns.  (12)  end  (13) 
sre  derived  by  multiplying  (12)  end  (13)  by  (Ps/M)noBinal  snd  (y«>„omlnalt 
respectively.  Also  note  from  Eqns.  (12)  end  (13)  thet 


T  3M 


M  3P- 


_S  3M_ 
H  3P„ 


number)  Is  shown  In  Figure  7. 


Dimensionless 

Influence 

Coefficients 


JM 


iPT'P« 


PT 


3 

or 


M 


p. 


ps 


(.sae) 


<.782) 


(.574) 


(.397) 


(14) 


The  change  in  the  Influence  coefficients  in  (14)  with  p  /p 

D  4 


(or  Mach 


VPT 


(0)(Msch  Number) 


Figure  7.  Influence  Coefficients  versus  Mach  Number  for  the  Given 
Equation 

We  can  conclude  from  Eqns.  (12)  through  (14)  and  Figure  7  that  a  Mach 
number  calculation  using  Eqn.  (11)  Is  quite  sensitive  to  errors  in  P  or 
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P  ,  particularly  for  low,  aubaoolc  Mach  numbers.  Tharafora,  it  la  not 

S 

aurprlalng  to  observe  In  Tabla  II  that  tha  uncertainty  of  tha  calculated 
value  AM29  la  significantly  larger  than  that  for  AM19,  AM20,  or  AM25.  Tha 
laportant  laaaon  In  thle  axaapla  la  that  calculation  procedure  aay  have  a 
atrong  Influence  on  the  accuracy  of  a  coaputed  value.  If  an  engineer  haa 
aaveral  optlona  for  calculation  procedure,  he  ahould  conalder  the  accuracy 
of  the  reault  In  eelectlng  tha  beat  option. 

B.  Kxaaple  2:  Coapraaaor  Efficiency 

Coapraaaor  efficiency  la  one  of  tha  Boat  laportant  calculated 
valuea  used  to  characterise  the  perforaance  of  the  aachlne;  It  can  have  a 
very  significant  Influence  on  the  perforaance  of  the  turbine  engine.  We 
have  said  In  an  earlier  section  that  the  valuea  In  the  Influence 
coefficient  aatrlx  depend  upon  the  aagnltudes  of  the  aeaaured  valuea 
(l.e.,  the  Influence  coefficient  values  will  vary  with  coapresaor  speed). 
In  this  exaaple  we  show  how  the  uncertainty  of  coapresaor  efficiency 
varies  with  coapresaor  speed,  and  why  It  Is  particularly  difficult  to  Bake 
an  accurate  deteralnat Ion  of  coapresaor  efficiency  at  low  speed. 

Coapreseor  efficiency  la  calculated  froa  aeasured  pressures  and 
teaperatures  at  the  compressor  Inlet  and  exit.  Assuming  an  average  value 
of  the  ratio  of  specific  heats  for  the  coapresaor , 7 >  adiabatic  efficiency 
can  be  calculated  using  the  equation 

T-l 

(PT3p\  Y  _  j 

PT19 /  (15 

ij  _  - 

TT30  _  j 
TT19 

Taking  partial  derivatives  of  Eqn.  (15)  with  respect  to  PT3Q  ,  Pn9  ,  tt30  t 
•«<*  If, o  . 
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a  n 


3P, 


T30 


(r32) "^1  At*1-) 

T19  \  T19/  /  yT19  J 


(16) 


3»l 

w 


T19 


(y-D/7 

PT19 


(17) 


(18) 


(19) 


The  dimensionless  Influence  coefficients  corresponding  to  Eqns>  (16) 


through  (19)  are 


and  3Tt19  '  \  V  )  «om. 

Note  that 


3r>  PT30 

3pt3o  n 


T 

T30 

T 

T19 


(20) 


3>?  TT30 
9TT30  ^ 


3>7  TT19 
3TT19  ** 


XT30 

TT19 


(21) 
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Eqns.  (20)  and  (21)  can  be  used  to  compute  Influence  coefficients  for 
expected  temperature  and  pressure  ratio  combinations  at  various  compressor 
speeds  and  nominal  efficiencies. 

Table  III  contains  pressure  rat los  (  PT3Q  / PT19 ) ,  temperature  ratios 
(TT30 /TT19 ) ,  and  nominal  efficiency  values  for  60  percent,  80  percent,  and 
100  percent  design  speed  for  normal  compressor  operating  conditions  at 
those  speeds.  Influence  coefficients  have  been  derived  for  this  data 
using  Eqns.  (20)  and  (21)  and  assuming  Y  >1.4.  Finally,  the  measurement 
errors  given  in  Table  I  and  the  computed  Influence  coefficients  have  been 
used  in  Eqns.  (4)  through  (6)  to  yield  the  uncertainty  in  the  efficiency 
calculation  for  the  speeds  in  Table  III. 


Table  III 

COMPRESSOR  EFFICIENCY  VERSUS  COMPRESSOR  SPEED 


T  Design  Speed 

100 

60 

/PT3o\ 

Pressure  Ratiolp - 1 

\  T19/ 

8.0 

4.1 

2.3 

- 

Temperature  Ratio  (  T30\ 
\TT19 / 

1.966 

1  .  598 

1 . 332 

Nominal  Efficiency  (%> 

84 

83 

81 

Efficiency 

Uncertainty  (2) 

m 

1  .04 

1  .57 

Note  that  efficiency  uncertainty  for  60  percent  speed  Is  twice  that  for 
100  percent.  The  conclusion  to  be  derived  from  this  example  la  that 
efficiency  calculations  at  low  pressure  and  temperature  ratios  require 
very  careful  measurements  (low  bias  and  precision  errors)  of  Inlet  and 


99 


USAFA-TR-83-2 


discharge  pressures  end  temperature s ,  because  the  Influence  coefficients 
ere  relatively  large  at  these  operating  conditions. 

Also  note  that  we  have  assuaed  that  bias  and  precision  errors 
(percent  of  reading)  are  the  aaae  for  80  and  60  percent  speed  as  they  were 
for  100  percent  speed  (l.e.,  the  aeasureaent  errors  in  Table  I  have  been 
used  for  all  calculation  results  shown  In  Table  III).  This,  in  fact,  is 
probably  not  a  good  asauaptlon,  because  aeasureaent  error  (percent  of 
reading)  la  a  function  of  where,  In  the  transducer's  range,  the 
aeasureaent  occurs.  Measurement  error  can  be  quite  high  as  a  percentage 
of  the  measured  value  at  the  low  end  of  the  transducer  range,  so  the  test 
engineer  may  be  forced  to  change  transducers  to  obtain  the  required  data 
accuracy . 

V.  The  Benefits  of  Using  Uncertainty  Analysis 

We  said  in  the  Introduction  to  this  paper  that  measurement  processes 
associated  with  an  experiment  will  inevitably  yield  measurement  error.  If 
that  Is  true,  then  a  presentation  of  test  data  is  Incomplete  without  a 
statement  concerning  the  accuracy  of  the  data.  Therefore,  one  of  the 
major  responsibilities  of  a  test  engineer  Is  to  establish  data  accuracy  as 
part  of  the  test  requirements  and  then  to  report  the  accuracy  after  the 
test.  Uncertainty  analysis  Is  the  tool  that  the  engineer  uses  to  do  this 
task.  The  test  objectives  and  procedure,  instrumentation  layout,  and 
calculation  procedures  together  with  the  uncertainty  analysis  allow  the 
engineer  to  deteralne  the  best  test  method  before  the  test  occurs. 

Then,  following  the  test,  the  data  can  be  compared  to  pre-test  uncertainty 
estlaates  as  a  validity  check  and  to  Identify  problem  areas.  An  laportant 
benefit  froa  the  effort  expended  In  doing  an  uncertainty  analysis  Is  the 
careful  test  planning  that  Is  an  Integral  part  of  the  analysis.  In  aost 
cases,  this  planning  increases  the  probability  of  a  successful  test. 
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Also,  the  uncertainty  analysis  provides  valuable  insight  Into  the  test 
results . 

VI .  Conclusions 

In  this  paper  wo  have  presented  the  fundamental  concepts  associated 
with  uncertainty  analysis  and  have  used  the  results  of  an  analysis  that 
was  done  for  a  centrifugal  compressor  test  to  Illustrate  the  following 
Important  points: 

1.  The  precision  error  of  a  measured  value  Is  Inversely  proportional  I 
the  square  root  of  the  number  of  repeated  (over  a  period  of  time)  or 
redundant  measurements  of  that  value.  Redundant  Instrumentation  is  very 
Important  for  reducing  measurement  uncertainty. 

2.  The  uncertainty  of  calculated  values  can  be  strongly  Influenced  by 
the  calculation  procedure;  If  various  procedures  are  available,  the 
engineer  should  determine  which  procedure  minimises  data  uncertainty. 

3.  The  uncertainty  of  a  measured  value  can  vary  greatly  depending  on 
where,  In  the  transducer's  range,  the  measurement  occurs.  Normally, 
measurements  at  the  low  end  of  a  transducer's  range  should  be  avoided. 

4.  The  Influence  coefficients  3C^/3M^  depend  on  measured  values. 
Therefore,  the  uncertainty  of  a  calculated  value  depends  not  only  on  the 
blss  and  precision  of  measured  values  used  In  the  calculation,  but  also  on 
the  magnitude  of  the  measured  values.  Influence  coefficients  must  be 
recalculated  for  each  set  of  measured  data. 

finally,  we  have  attempted  to  explain  the  Importance  of  uncertainty 
analysis.  Data  accuracy  should  be  specified  as  part  of  the  test 
requirements,  estimated  using  uncertainty  analysis  methods,  and  then 
reported  after  the  test.  The  careful  test  planning  that  is  an  integral 
part  of  the  analyels  Increases  the  probability  that  test  objectives  will 
be  clearly  defined  and  accomplished;  that,  in  Itself,  is  strong 
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Greek  Symbol  8 

6 

Y 

Y 

n 

o 

Subscripts 


1 9 , 30 , e tc  . 


Airflow 

A  measurement  In  the  calculation  of  S 

The  average  of  measurements  In  the  calculation  of 

S 


Flow  Swirl  Angle 
Bias  Error 

Ratio  of  Specific  Heats 

T  + 

Ratio  of  Specific  Heats  evaluated  at  —  - 

Adiabatic  Efficiency 
Precision  Error 

Pertaining  to  the  l1"*1  computed  value, 
Pertaining  to  the  jth  measured  value,  M^ 
Pertaining  to  measurement  stations  shown  In 
Figure  2 
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INTEGRATION  OF  AN  AIRBORNE  LABORATORY  INTO 
THE  UNITED  STATES  AIR  FORCE  ACADEMY  ACADEMIC  CURRICULUM 

Kent  R.  Crenshaw* 

Abstrac t 

The  Department  of  Aeronautics  at  the  United  States  Air  Force  Academy 
has  introduced  a  new  special  topics  course  —  Aero  49S,  Flight  Test 
Techniques  —  into  its  acadealc  curriculum  for  the  Fall  semester  of  19S2. 
By  using  leased  aircraft  piloted  by  faculty  leakers,  the  course  attempts 
to  show  students  the  relationship  between  the  theory  of  flight  mechanics 
and  actual  aircraft  performance  and  flying  qualities  (specifically, 
stability  and  control).  Each  cadet  enrolled  in  the  course  receives  four 
flights  in  a  light,  single- eng i ne  aircraft:  two  flights  to  determine 
aircraft  performance  (i.e.,  turns,  cruise,  climbs  and  descents)  and  two 
flights  to  evaluate  aircraft  flying  qualities  (i.e.,  static  and  dynamic 
stability  and  stalls).  Special  cockpit  instrumentation  was  designed  and 
built  to  measure  static  stability  parameters.  This  instrumentation 
consists  of  a  force  gauge,  to  determine  stick  force,  and  rudder  pedal 
pressure  plates  using  strain  gauges  to  measure  rudder  forces.  The  flying 
program  has  been  ent husla at  leal 1 y  received  by  cadets  and  faculty,  and  its 
potential  for  use  in  other  courses  and  for  research  is  enormous. 


I .  Introduction 

The  facilities  of  the  Departaent  of  Aeronautics  Laboratory  at  the 
USAF  Academy  are  without  a  doubt  some  of  the  best  in  the  country.  The 
subsonic  and  trlsonlc  wind  tunnel  facilities,  along  with  the  engine  test 
cells,  provide  a  vital  link  between  classroom  theory  and  real-life, 
practical  application  of  the  principles  of  aeronautics.  Vhlle  the 
instructional  method  currently  used  by  the  Department  tries  whenever 
possible  to  link  theory  with  reality,  as  in  the  aerodynamics  and 
propulsion  areas  of  aeronautics  instruction,  a  laboratory  course  has 
always  been  needed  to  complement  the  study  of  flight  mechanics.  The 
"Balsa  Glider  Lab,"  required  in  Aero  356,  Flight  Mechanics  I,  is  currently 
the  only  laboratory  exercise  offered  to  cadets  in  the  field  of  flight 
mechanics.  The  "Glider  Lab,"  accomplished  during  the  last  quarter,  allows 
cadets  to  culminate  their  study  of  flight  mechanics  by  applying  basic 
aircraft  design  principles  to  sstlsfy  fundamental  aircraft  performance  and 
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stabll ity-and-contr ol  (S  &  C)  considerations.  The  balaa  glider  laboratory 
experience  has  always  been  enthusiastically  received  by  both  students  and 
faculty . 

Providing  the  student  at  the  USAF  Academy  with  actual  experience  In 
gathering  perforaance  and  atablllty-and-control  data  in  order  to  define 
these  characteristics  for  a  full-scale  aircraft  has  been  neglected.  To 
aeet  this  need,  the  Department  of  Aeronautics  currlculua  needed  an 
Airborne  Laboratory  to  link  classrooa  theory  with  the  knowledge  of 
aircraft  flight  behavior  gained  froa  actual  flight. 

A.  Background  and  Objectives 

Atteapts  have  been  made  In  the  past  to  introduce  a  course  using  an 
airborne  laboratory  Into  the  curriculum.  These  attempts  have  Involved  an 
assortaent  of  aircraft,  including  the  U-10  Hello-Courier,  the  U-4 , 
sailplanes,  and  the  T-37B.  None  of  these  endeavors  was  successful, 
prlaarlly  because  the  aircraft  were  not  available  when  we  needed  them. 
Other  contributing  factors  Included  conflicts  In  scheduling,  financial 
restraints,  aircraft  accidents,  and.  In  general,  failure  to  plan  a 
comprehensive  flying  program  that  could  be  Integrated  with  appropriate 
classrooa  Instruction.  Even  though  these  atteapts  were  failures,  the 
students,  faculty,  and  graduates  were  overwhelmingly  In  favor  of  such  a 
cour  se . 

In  order  to  find  a  solution  to  this  old  problem,  we  began  In  January 
1982  to  foraulate  an  airborne  laboratory  prograa  and  to  evaluate  different 
aircraft  to  aeet  our  prograa  needs.  Our  objectives  were  to  create  a 
uniquely  cosprehenslve  currlculua,  enhance  student  motivation  and 
engineering  experience,  and  develop  and  lapleaent  a  prograa  In  tlae  for 
the  Fall  1982  seaester. 
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B.  Prograa  Proposal* 

In  planning  our  course  of  action,  we  exaalned  existing 
aeronautical  education  prograas  at  other  Institutions  around  the  country. 
The  United  States  Naval  Acadeay,  Pennsylvania  State  University, 

Mississippi  State  University,  the  University  of  Arltona,  and  the  United 
States  Military  Acadeay,  to  naae  a  few,  were  schools  using  airborne 
laboratories  in  their  currlculuas.  Out  of  all  the  prograas  we  studied, 
the  United  States  Naval  Acadeay  prograa  appeared  to  best  suit  our 
requlreaent s .  In  addition,  the  Navy's  willingness  to  send  us  their 
course  aaterlals.  Including  specific  descriptions  of  how  they  Integrate  a 
flying  laboratory  into  their  Aeronautics  education  prograaa,  was 
particularly  helpful. 

After  reviewing  in  detail  the  Naval  Acadeay  prograa,  coaparlng  It  to 
prograaa  at  other  schools,  and  analysing  our  own,  we  decided  to  foraulate 
a  new  course  using  a  leased,  FAA-cer tif led  passenger  aircraft  with  our  own 
faculty  aeabers  as  pilots.  In  order  to  fully  appreciate  the  context 
within  which  this  decision  was  asde,  soae  discussion  of  the  prograa  and 
aircraft  options  we  looked  at  Is  necessary. 

The  two  prograa  proposals  presented  to  the  head  of  the  Departaent  of 
Aeronautics  and  his  staff  suggested  the  following  alternatives: 

1.  Integration  of  the  airborne  laboratory  Into  an  existing  course 
(Aero  356,  450,  457,  or  464). 

2.  Creation  of  a  new  course.  Aero  495,  entitled  "Flight  Test 
Techniques  for  Plight  Mechanics." 

The  first  option,  which  added  aore  aaterlal  to  courses  already  covering  a 
vast  rang*  of  topics,  was  ellalnated  becaus*  the  tla*  required  to 
accoaaodate  a  flying  prograa  would  detract  froa  the  aaln  objectives  of  the 
existing  courses  as  well  as  those  of  the  airborne  laboratory  itself. 
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In  contract,  option  2  allowed  us  to  develop  a  truly  comprehensive  new 
laboratory  program.  The  potential  scope  of  a  new  495  course  appeared  to 
be  wore  than  adequate  to  wake  a  semester  course  offering. 

Our  attention  then  turned  to  the  problen  of  determining  which 
aircraft  to  use  in  support  of  the  laboratory.  Criteria  for  selecting  an 
aircraft,  based  on  evaluations  of  safety,  cost,  convenience,  and 
practicality,  were  finally  established  as  follows: 

1.  an  engine  size  of  at  least  180  horsepower 

2.  four-seat  capacity 

3.  single  engine 

4.  low  wing 

5.  radio  and  transponder  equipped 

6.  available  spare  aircraft  of  the  same  type 

7.  ability  to  base  the  aircraft  at  the  USAFA  Airfield,  and 

8.  established  fuel  and  maintenance  support 

The  selection  of  an  aircraft  to  meet  these  criteria  was  narrowed  down  to 
four  possible  options: 

1.  use  of  the  Air  Training  Command  CATC)  T-41C  aircraft 

2.  use  of  the  USAFA  Aero  Club  aircraft 

3.  purchase  of  our  own  aircraft,  and 

4.  lease  of  a  commercial  aircraft 

Option  1  Involved  problems  we  have  had  In  the  past  using  somebody 
else' a  airplane  In  the  midst  of  somebody  else's  flying  program.  The 
potential  for  scheduling  conflicts  and  the  amount  of  time  ATC  required  for 
their  attached  pilots  to  check  out  and  remain  current  in  the  T-41C 
eliminated  that  option. 

Option  2,  while  attractive  from  a  convenience  and  cost  standpoint, 
was  eliminated  because  of  potential  maintenance  problems  and  lack  of 
scheduling  priority.  The  Aero  Club,  as  a  "Morale,  Welfare,  and  Kecreatlon 
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( MUR ) "  Prograa,  would  not  give  ue  priority  In  echedullng  their  aircraft, 
particularly  If  the  elrcreft  we  deeired  wee  echeduled  by  e  club  aeaber  for 
e  croee-country  flight.  The  flexibility  thet  would  be  required  on  our 
pert  to  uee  thle  option  wee  thought  to  be  detriaental  to  the  continuity  of 
our  leb  prograa.  In  addition,  the  Aero  Club  had  only  two  available 
aircraft  that  aatiefled  our  perforaance  criteria:  the  Beech  Sierra  and  the 
Piper  Arrow.  The  Sierra,  an  older  aircraft,  waa,  and  atill  la,  plagued  by 
both  ainor  and  aajor  aalntenance  probleaa.  Since  the  Piper  Arrow  la  the 
club'a  prlaary  cr osa-countr y  aircraft,  It  la  frequently  unavailable.  The 
other  club  aircraft,  while  aore  available,  failed  our  perforaance 
criteria. 

Option  3  offered  ua  the  aoat  control  of  an  aircraft,  but  It  alao  gave 
ua  the  moat  r e aponalb 11 1 ty  and  added  to  our  workload.  In  addition,  the 
lead  tlae  neceaaary  to  coordinate  acqulaltlon  of  an  aircraft  with  the  Air 
Staff  la  meaaured  In  yeara.  Since  our  objective  waa  to  wake  the  Airborne 
Laboratory  operational  by  the  Fall  1982  eeaeeter,  thin  option,  for  now  at 
leant,  waa  alao  ellalnated. 

The  aoat  viable  option,  and  the  one  uaed  by  aoat  of  the  achoola  we 
contacted  that  have  current  airborne  lab  prograaa,  waa  to  leaae  a 
coaaerclal  aircraft.  Hedrick  Beechcraft  at  the  Colorado  Sprlnga  Municipal 
Airport  had  the  type  of  aircraft  we  were  looking  for  and  aore  than  juat 
one  of  each  type  of  aircraft.  In  addition,  their  aalntenance  and  fuel 
aupport  operatlona  appeared  to  be  outatandlng.  Hedrick  had  two  Beech 
Bonansaa  rated  at  263  horaepover,  two  Beech  Slcrraa  at  200  horaepower,  and 
four  Beech  Sundownere  at  180  horaepower  available  for  leaae.  While 
leeelng  froa  Hedrick  would  coat,  on  the  average,  $20.00  per  hour  aore  than 
leaving  any  of  the  Acadeay  Aero  Club  aircraft,  all  but  one  of  our  aircraft 
criteria  were  eatlafled  by  their  aervlce.  The  coaproalee  here  waa  not 
being  able  to  operate  out  of  the  USAFA  Airfield.  In  practice,  thle  haa 
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been  lees  of  a  problem  than  anticipated.  Uhlle  the  USAFA  Airfield  may  be 
closed  because  of  Inclement  weather  or  winds,  the  Colorado  Springs  airport 
is  usually  open  (visual  flight  rule  conditions)  and  has  several  more 
available  runways.  In  addition,  conflicts  with  USAF  Academy  T41 
operations,  particularly  in  the  traffic  pattern,  are  avoided. 

The  decision  concerning  the  type  of  aircraft  to  lease  from  Hedrick 
was  based  on  cost.  The  Bonanza,  at  $85.00  per  hour,  was  eliminated  from 
our  planning  because  it  was  considered  too  expensive.  The  Sierra,  at 
$55.00  per  hour,  was  thought  to  offer  an  acceptable  performance  capability 
for  our  required  flight  operations.  The  Sundowner,  at  $40.00  per  hour, 
was  found  to  be  adequate  for  the  flying  qualities  phase  of  the  laboratory. 
Both  aircraft  are  aer odynamlcal 1 y  the  same;  the  Sundowner  Is  the 
f ixed-pl t ch-pr op  and  fixed-gear  version  of  the  Sierra.  Both  are  shown  in 
Figure  1. 

After  receiving  the  approval  of  the  Head  of  the  Department  of 
Aeronautics,  the  program  proposal  for  a  new  Aero  495  course  titled  "Flight 
Test  Techniques"  and  the  proposal  to  lease  aircraft  from  Hedrick 
Beechcraft  was  presented  to  the  Dean  and  the  Superintendent  of  the  USAF 
Academy  during  February  and  March  of  1982.  The  program  was  approved 
without  revision  by  both,  and  the  Dean  supplied  matching  funds  to  finance 
the  program  through  the  end  of  FY82  .  At.  this  point  the  decision  was  also 
made  to  enroll  only  ten  cadets  in  the  course  for  both  the  Fall  and  Spring 
semesters . 

The  following  course  description.  Including  the  assignment  of  credit 
and  contact  hours,  was  submitted  to  the  Cadet  Curriculum  and  Scheduling 
Office  of  the  USAF  Academy: 

Aero  495.  Flight  Test  Techniques.  1(2) 

This  course  will  cover  fundamental  flight  test  methods  for  gathering 
performance  and  flying  qualities  data  for  fixed-wing  aircraft.  The  theory 
behind  each  flight  test  method  will  be  studied  using  fundamental  and 
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advanced  principles  of  flight  mechanics.  Related  topics  will  consist  of 
pltot-statlc  calibration,  1 n s t r use n t a 1 1  on ,  data  reduction,  and 
welght-and-balance .  A  light,  single-engine  aircraft  Is  used  as  an 
airborne  laboratory  for  obtaining  data  In  flight.  Pinal  exam. 
Prerequisite:  Completed  Aero  356.  Semester  Hours:  3. 


With  this  background  on  the  evolution  of  our  program  In  mind,  It 
might  be  useful  to  provide  a  detailed  description  of  the  scope  of  the 
course,  the  instructional  method  that  we  intend  to  use,  our  scheduling 
plans,  our  Instrumentation  support,  and  the  current  status  of  the  flying 
program. 


II .  Aero  495.  Flight  Test  Techniques 
A.  Course  Overview 

Aero  495  Is  divided  into  two  phases,  performance  and  flying 
qualities,  which  are  discussed  more  fully  below.  Each  cadet  enrolled 
receives  four  flights  during  the  semester,  two  In  each  phase.  The  Beech 
dlerra  Is  used  in  the  performance  phase  and  the  Beech  Sundowner  during  the 
flying  qualities  phase. 


1 .  Per f or mance 

The  performance  phase  Is  covered  In  lessons  1  through  23.  The 
following  topics  are  dealt  with  during  this  phase: 

a.  the  basics  of  performance  flight  testing, 

b.  standard  and  non-standard  atmospheric  conditions, 

c.  calibration  and  use  of  the  pltot-statlc  system, 

d.  reciprocating  engine  mechanics  and  propeller  theory, 

e.  Implication  and  analysis  of  aircraft  weight  and  balance, 

f.  the  flight  test  techniques  used  In  evaluation  of  each  phase 
of  an  aircraft's  typical  flight  profile  (takeoff,  climb, 
cruise,  level  turn,  descent,  and  landing). 
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g.  the  analysis  and  importance  of  a  teat  plan, 

h.  aircraft  Instrumentation  and  general  cockpit  layout,  and 
1.  ground  and  ln-flight  safety 

The  primary  text  used  in  this  part  of  the  course  Is  Light  Aircraft 
Performance  for  Test  Pilots  and  Flight  Test  Engineers  by  Sean  C.  Roberts 
of  Flight  Research,  Inc.,  Mojave,  California.  Supplemental  readings 
consist  of  notes  provided  by  the  USAF  Test  Pilot  School  and  by  Mr.  Ralph 
D.  Klmberlln  of  the  University  of  Tennessee  Space  Institute  at  Tullahoma, 
Tennessee.  Excerpts  taken  from  the  Pilots'  Operating  Handbook  for  the 
Beechcraft  Sierra  C24R  are  also  used.  The  Roberts  text  is  particularly 
good  since.  In  addition  to  covering  most  of  the  classroom  topics  In 
detail.  It  provides  formats  for  in-flight  data  cards  and  data  reduction 
pr  ocedure  s . 

Of  particular  Interest  to  the  cadets  Is  the  propulsion  aspect  of  the 
course.  While  the  cadets  have  had  at  least  some  exposure  to  most  of  the 
classroom  topics  In  the  performance  phase  of  Aero  495,  they  have  had 
little  or  no  experience  with  propellers.  All  of  our  flight  mechanics  and 
propulsion  courses  address  gas-turbine  engine  technology  only.  This,  of 
course,  Is  because  most  of  the  aircraft  in  the  active  Air  Force  inventory 
are  powered  by  gas  turbines  (jets). 

The  flying  portion  of  the  performance  phase  Is  conducted  according  to 
a  test  plan,  which  has  a  format  similar  to  that  used  at  the  Air  Force 
Flight  Test  Center.  The  test  plan  defines  specific  performance  objectives 
that  must  be  met  If  the  advertised  performance  of  the  test  aircraft  is  to 
be  verified.  Performance  parameters  thst  are  evaluated  Include  maximum 
speed,  range  and  rate-of-climb  capability,  service  ceiling,  and  glide 
ratio.  The  test  plan  also  directs  that  the  cadets  be  divided  Into  two 
separate  test  teams,  as  shown  In  Figure  2. 
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Teat  Team  A  Teat  Team  B 


Project  Engineera  Project  Engineera 

Figure  2.  Test  Plan  Organization 

The  cadets,  functioning  as  project  engineers  under  the  supervision  of  a 
test  director,  prepare  data  cards,  prioritize  the  data  to  be  collected, 
and  analyze  the  data  after  each  test  flight.  The  test  plan  for  the 
"Sierra  C24R  Limited  Performance  Evaluation"  Is  shown  in  Appendix  A. 

Both  performance  evaluation  flights  last  one  hour,  with  flight  1 
dedicated  to  gathering  aircraft  cruise  and  turn  performance  data  and 
flight  2  to  aircraft  climb  and  descent  data.  A  "Flight  Test  Planning 
Guide,"  modeled  after  the  format  used  at  the  Naval  Academy,  Is  provided  to 
assist  the  cadets  In  their  preparation  for  each  flight.  Mission  events, 
pilot  and  student  responsibilities,  and  post-flight  data  reduction 
requirements  are  clearly  defined.  In  addition,  data  sheets,  data 
reduction  sheets,  and  an  "Initial  Flight  Test  Report"  form  modeled  after 
the  Air  Force  Flight  Test  Center  (AFFTC)  Form  365  are  also  used. 
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After  each  flight,  the  cadets  submit  for  grading  a  flight  report  that 
satisfies  the  requirements  laid  out  in  the  "Flight  Test  Planning  Guide." 
Other  graded  events  in  the  performance  phase  Include  a  graded  review  and  a 
formal  oral  report  presented  by  each  test  team  on  lesson  23. 

2.  Flying  Qualities 

The  flying  qualities  phase  of  Aero  495  consists  of  lessons  24 
through  41.  Classroom  topics  include  the  following: 

a.  a  definition  of  the  term  and  an  explanation  of  the 
importance  of  determining  an  aircraft's  flying  qualities, 

b.  the  study  of  flying  qualities  requirements  for  military 
aircraft  as  defined  in  the  Military  Specification  for 
Flying  Qualities  of  Piloted  Airplanes  (MIL-F-878 5C )  ,  and 

c.  the  flight  test  techniques  needed  to  evaluate  the  flying 
qualities  of  an  aircraft  (lateral-directional  static 
stability,  dynamic  stability,  longitudinal  static  stability, 
stall  characteristics,  and  maneuvering  flight 
characteristics) . 

The  primary  text  for  the  Flying  Qualities  phase  is  another  book  by  Sean  C. 
Roberts,  titled  Flying  Qualities  Flight  Testing  of  Light  Aircraft  for 
Teat  P 1 1 o  1 8  and  Flight  Test  Engineers.  A  supplemental  text  is  also 
provided  which  includes  chapters  4,  5,  6,  and  7  of  the  United  States  Air 
Force  Test  Pilot  School  notes. 

The  flying  portion  of  the  flying  qualities  phase  is  also  conducted 
according  to  a  test  plan.  The  objectives  are  to  qualitatively  and 
quantitatively  evalur.te  the  Beech  Sundowner  180  C23  as  a  primary  trainer 
for  Class  I  as  defined  in  MIL-F-8785C.  The  aircraft  is  evaluated  for 
compliance  with  selected  paragraphs  of  MIt,-F-87 8 5C .  The  evaluation  effort 
is  again  organized  as  shown  in  Figure  2. 
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As  in  the  performance  phase,  both  flights  in  the  tlying  qualities 
phase  are  for  one  hour.  Flight  3  is  dedicated  to  evaluating  longitudinal 
and  lateral-directional  stability  and  control  as  well  as  maneuvering 
flight.  Flight  4  concerns  dynamic  stability  and  stalls.  A  "Flight  Test 
Planning  Guide”  that  covers  each  flight  is  again  provided  to  assist  in 
cadet  preparation  for  flying. 

Graded  events  again  include  two  flight  reports,  a  graded  review,  and 
a  formal  oral  report  by  each  test  team.  The  semester  ends  with  a  review 
and  course  evaluation  on  lesson  42  and  a  final  exam  worth  25  percent  of 
the  cadet's  final  grade. 

For  more  detail  on  the  actual  test  plans,  planning  guides,  and  data 
reduction  techniques  used  in  both  phases  of  the  course  as  well  as  sample 
calculations,  refer  to  USAFA-TN-83-3 . 

B.  X n s t r uc t i ona 1  Method 

The  instructional  method  used  in  Aero  495,  shown  in  Figure  3,  is 
modeled  closely  after  the  method  used  at  the  United  States  Air  Force  Test 
Pilot  School. 


Figure  3.  Instructional  Method 
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Note  that  the  flight  test  techniques  are  covered  during  classroom  lectures 
as  they  apply  to  each  area  of  study,  l.e.,  cruise,  climbs,  longitudinal 
static  stability,  etc.  While  the  techniques  addressed  are  for  light 
aircraft  testing,  other  techniques  for  higher  performance  and  multi-engine 
aircraft  are  also  studied. 

Data  reduction  for  the  first  phase  of  the  course  involves 
standardizing  performance  data  to  standard  aircraft  weights  and 
atmospheric  conditions.  During  the  flying  qualities  phase  little  data 
reduction  is  required.  The  data  reduction  schemes  for  both  phases  of  the 
course  are  covered  In  lesson  readings,  class  lectures,  and,  more 
spec  1 f ic al 1 y ,  In  the  flight  test  planning  guides  mentioned  earlier. 

The  performance  phase  of  the  course  lends  itself  very  well  to  using 
computer  routines  for  data  reduction,  since  doing  it  by  hand  is  a  lengthy 
pro-jss.  However,  in  the  interest  of  the  cadets  learning  the  procedures 
by  doing,  they  have  not  been  provided  with  any  computer  codes  to  assist 
them  In  reducing  data.  Realizing  and  accepting  the  challenge,  two  cadets 
this  semester  have  written  data  reduction  programs  for  the  performance 
phase  of  the  course.  Cadets  will  be  permitted  to  use  these  programs  in 
the  future,  but  only  as  a  basis  for  checking  their  work. 

C.  Instrumentation 

In  the  performance  phase  of  this  course,  the  Beech  Sierra  had  the 
standard  cockpit  gauges  for  gathering  all  necessary  parameters.  The  only 
additional  device  required  was  a  stopwatch  for  timing  turns,  climbs,  and 
descents.  Table  I  shows  the  parameters  required  for  each  performance  data 
flight . 

The  flying  qualities  phase,  for  which  the  Beech  Sundowner  was  used, 
required  more  than  just  the  standard  cockpit  instruments.  Special 
Instrumentation  and  procedures  were  necessary.  Stick  force  and 
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displacement  had  to  be  measured  as  well  as  rudder  force  and  displacement, 
an  accelerometer  or  "g"  meter  had  to  be  acquired,  and  a  means  of  center  of 
gravity  control  was  nneeded. 

A  "g"  meter  was  borrowed  from  the  Academy's  Soaring  Branch  of 
Airmanship  and  mounted  on  top  of  the  aircraft  glare  shield,  as  shown  in 
Figure  h. 


Figure  4.  "G"  Meter  Installation 

The  various  displacements  and  rotations  of  the  controls  that  we  needed 
were  measured  as  follows:  (1)  longitudinal  "stick"  (control  push  rod) 
displacement  was  determined  by  using  a  tape  measure,  as  shown  In  Figure  5, 
(2)  a  yardstick  was  used  to  measure  rudder  displacement,  as  shown  In 
Plgure  6,  and  (3)  the  rotation  angle  of  the  yoke  was  obtained  by  marking  a 
reference  line  on  top  of  the  push  rod  and  placing  a  scale  on  the 
Instrument  panel-  This  reference  scale  was  marked  off  in  Increments  of 
one-fourth  from  full  left  rotation  to  full  right  rotation.  This 
arrangement  Is  shown  In  Figure  7. 
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Figure  7.  Measuring  Yoke  Angle 


These  rather  simple  methods  proved  effective.  However,  we  also  needed  to 
measure  the  force  required  to  hold  the  "stick"  and  rudder  pedals  in 
various  positions.  The  "stick"  force  was  measured  mechanically  by  a 
device  mounted  in  the  cockpit,  as  depicted  in  Figure  8. 
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The  force  gauge  taunted  on  the  "stick"  tree  calibrated  up  to  appr oximately 
IS  pounds.  Rudder  pedal  force  required  more  sophisticated 
lnatruaentat Ion .  Strain  gauge*  Installed  In  pressure  plate*  mounted  on 
the  rudder  pedal*  with  velcro  strap*  *ent  six-volt  signals  to  two  meters 
(one  readout  for  each  pedal)  mounted  In  a  box  held  by  the  cadet  In  the 
back  seat.  The  system  was  powered  by  the  aircraft  electrical  system 
(either  14  or  28  volts)  through  a  voltage  regulator  using  the  cigarette 
lighter  on  the  Instrument  panel*  Rudder  pedal  force  readings  were 
calibrated  up  to  SO  pounds*  Figures  9  and  10  show  the  rudder  pedal 
pressure  plates  and  the  six-volt  readout  system. 


Top  Pressure  Plate 
Slide*  on  Track 


TOP  VIEW 
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Notes: 


r=L 


TWO  WAY  THROW  SWITCHES 


Era 


POTENTIOMETERS  FOR  CALIBRATION 


1.  R  x  1  scale  used  for  smaller,  more  sensitive  force  readings.  R  x  2  scale  used 
for  larger  readings  up  to  50  lbs. 

2.  Ammeters  read  In  milliamps  converted  to  force  In  pounds  from  5  to  50  lbs. 


Figure  10.  Rudder  Force  Readout  System 
Flight  3  required  the  cadets  to  gather  longitudinal  and  maneuvering 
flight  data  with  the  center  of  gravity  located  first  In  a  forward  position 
and  then  later  In  an  aft  position.  Center-of-gravlty  control  was 
accomplished  by  using  two  fifty-pound  lead  weights  secured  In  the  baggage 
compartment.  While  cadets  were  required  to  check  the  center  of  gravity 
prior  to  each  flight  to  Insure  that  It  was  within  flight  manual  limits, 
this  became  particularly  important  when  lead  ballast  was  carried. 

Table  II  shows  the  parameters  required  for  each  flight  In  the  flying 


qual ltles  phase  . 
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Table  II 

PARAMETERS  REQUIRED  FOR  FLYING  QUALITIES  DATA 

Flight  3.  Longitudinal  and  Lateral  Directional  SAC; 
Maneuvering  Flight 


Longitudinal  Static 
Stability  Data: 


Lateral-Directional 
Stability  Data: 


Indicated  Airspeed 
Indicated  Altitude 
Outside  Air  Teaperature 

*Tach  Tlae 
Engine  RPM 
Longitudinal  Control 
Displacement  ("Stick”) 
Longitudinal  Control 
Force  ("Stick") 

Maneuvering  Flight  Data 


Indicated  Airspeed 
Indicated  Altitude 
Outside  Air  Teaperature 
*Tach  Time 
Engine  RPM 

Aircraft  Load  Factor  ("g"  meter) 
Longitudinal  Control 
Displacement  ("Stick”) 
Longitudinal  Control 
Force  ("Stick") 


Indicated  Airspeed 
Indicated  Altitude 
Outside  Air  Teaperature 
*Tach  Tlae 
Engine  RPM 
**Aircraft  Heading 
Bank  Angle 

Longitudinal  Control 
Force  ("Stick") 

Rudder  Force 

Yoke  Angle 

Rudder  Displacement 


*Used  to  extrapolate  aircraft  weight  at  test  points,  since  fuel 
quantity  gauges  are  not  accurate. 

**The  heading  Indicator,  assuming  an  Initial  trimmed  and  coordinated 
flight  condition,  was  used  to  measure  sideslip  angle. 


Flight  A.  Dynamic  Stability 


Stall s 


Dynamic  Stability  Data: 


Indicated  Airspeed 
Indicated  Altitude 
Outside  Air  Temperature 
•Tach  Time 
Engine  RPM 
Tlae  Oscillations 
Roll  Rates 
Bank  Angle  Changes 


Stall  Data: 

Indicated  Airspeed 
Indicated  Altitude 
Outside  Air  Temperature 
*Tach  Time 
Engine  RPM 


•Used  to  extrapolate  aircraft  weight  at  test  points,  since  fuel 
quantity  gauges  are  not  accurate. 
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In  summary,  the  instrumentation  was  crude,  but  it  worked.  Aa  a 
spinoff,  instrumentation  to  measure  stick  and  rudder  pedal  forces  is  being 
designed  by  Engineering  430  students.  The  Sundowner,  in  its  capacity  as 
the  Aero  495  flying  qualities  test  bed,  will  be  used  to  evaluate  that 
equipment  later  this  semester. 

D.  Cadet  Scheduling  and  Transportation 

With  the  ten  cadets  enrolled  in  Aero  495  making  up  two  test  teams 

and  two  cadets  per  flight,  five  data  flights  are  flown  for  each  of  the 

four  missions  just  discussed.  Each  of  the  test  teams  plana  and  flies  two 
data  flights  (two/team)  independently,  while  the  fifth  flight  is  a 
combined  effort  flown  with  both  test  directors  to  fill  in  gaps  in  the  data 
gathered  earlier  by  other  test  team  members  or  to  repeat  selected  data 
points.  (See  Figure  2  for  Test  Plan  Or ga nl za 1 1  on .  )  While  two  cadets  are 

flying,  the  other  eight  meet  in  the  classroom  for  one  hour  to  prepare  for 

future  flights,  reduce  data,  or  work  on  flight  reports.  Aero  495  is  a 
two-contact-hour  course  and  meets  for  two  hours  only  during  regularly 
scheduled  classroom  lecture  periods.  Cadets  who  are  flying  use  the  two 
hours  of  scheduled  class  time  on  the  flight  line,  and  their  early  morning 
military  training  time  for  traveling  to  the  airport.  Figure  11  shows  the 
Daily  Flying  Schedule  as  arranged  during  the  Fall  1982  semester. 

Time  Event 

0700  Military  Taxi  to  Hedrick  Beechcraft 

0740  Military  Taxi  Arrivea  at  Hedrick  Beechcraft 

0800  Takeoff  and  Fly  Lab  Misaion 

0900  Land  and  Debrief  Lab  Miaaion 

0910  Tranaportation  back  to  Fairchild  Hall 

Figure  11.  Daily  Flying  Schedule 
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Approximately  80  minutes  are  spent  on  the  road  traveling  to  and  from  the 
airport.  Transportation  to  the  airport  1b  by  military  taxi,  while  return 
transportation  is  provided  by  the  Aero  495  pilot  with  whom  the  cadets 
flew.  Return  travel  time  is  used  to  debrief  the  flight.  Cadets  are 
prebriefed  in  the  classroom  before  flying  on  aircraft  orientation  and 
safety.  The  safety  and  mission  briefing  is  repeated  again  down  at  the 
flight  line.  I nc id e nt a  1 1 y ,  the  faculty  pilot  normally  arrives  at  the 
flight  line  at  least  an  hour  prior  to  flying  in  order  to  preflight  the 
aircraft,  check  the  weather,  and  Install  any  necessary  instrumentation. 

E.  Flying  Program  to  Data 

As  this  paper  is  being  written,  we  are  more  than  three  quarters  of 
the  way  into  the  Fall  1982  semester,  and  the  first  three  data  missions 
have  been  completed.  Table  Ill  is  a  summary  of  our  flying  record  to  date. 

Table  III 

AERO  495  SORTIE  SUMMARY 


Instructor  Only  Sorties  7 

•Cadet  Sorties  16 

Maintenance  Cancellations  1 

•Weather  Cancellations  _5 

Total  Sorties  2_3 

Total  Flying  Hours  27.5 


•One  flight  air  aborted  due  to  weather. 

Flexibility  in  the  flying  schedule  to  compensate  for  weather  cancellations 
and  the  single  cancellation  for  maintenance  we  experienced  was  achieved  by 
flying  four  cadets  per  day  using  two  airplanes  and  both  course  pilots.  In 
addition,  during  down  time  for  weather,  the  classroom  lecture  schedule  was 
moved  up,  which  allowed  for  makeup  flying  days  in  place  of  future  class 
meetings.  Using  more  than  one  airplane  per  flying  period  worked  well 
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during  the  performance  phase,  since  no  special  1 n st r ument at  Ion  was 
required.  In  the  flying  qualities  phase,  however,  with  only  one  set  of 
lnstruaentation,  this  option  was  not  a  good  source  of  flexibility, 
particularly  for  flight  3,  which  required  the  use  of  elevator  and  rudder 
pedal  force  measuring  devices.  Flight  4,  consisting  of  dynamic  maneuvers, 
rolls,  and  stalls,  did  not  require  special  instrumentation  other  than  a 
means  for  measuring  yoke  angle  and  a  stop  watch  for  timing.  The  only  two 
workable  options  for  flights  3  and  4  makeups  were  to  adjust  the  classroom 
lecture  schedule  and  to  arrange  makeup  flights  on  a  volunteer  basis  during 
cadet  free  time.  The  last  option,  however,  was  the  least  desirable  and 
was  used  only  once  during  the  first  15  flights.  Our  prime  objective, 
nonetheless,  was  that  all  cadets  fly  all  the  data  flights  for  which  they 
were  scheduled  . 

To  summarize,  the  course  has  proceeded  smoothly  and  as  scheduled  this 
semester,  weather  and  some  minor  instrumentation  problems  notwithstanding. 
The  scope  and  technical  level  of  the  course  have  been  more  than  adequate 
for  an  undergraduate  level  of  study  and  the  ratio  of  classroom  to  flying 
time  appears  to  be  well-balanced.  On  a  personal  note,  cadet  comments  on 
the  course  have  been  uniformly  positive. 

F.  Potential  Improvement  Areas 

Both  instructors  and  cadets  have  identified  areas  in  which  course 
Improvement  is  needed.  For  example,  the  aircraft  performance  envelope 
that  the  cadets  are  evaluating  needs  to  be  increased.  To  accomplish  this, 
we  would  require  a  bigger,  more  powerful  aircraft;  the  Beech  Bonanza, 
rated  at  285  horsepower,  is  under  consideration.  With  our  current  plan  of 
60  flights  per  year  (40  cadet  sorties  and  20  pilot  proficiency  sorties), 
another  $2,000.00  per  year  added  to  our  current  budget  of  $3,500.00  would 
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cover  the  cost  of  using  the  Bonanza,  possibly  with  three  cadets  per  flight 
Instead  of  the  usual  two. 

Another  advantage  of  the  airborne  laboratory  that  should  be  pursued 
is  the  use  of  this  capability  to  allow  cadets  In  other  courses  (Aero  499, 
Aero  450,  Engineering  420,  etc.)  to  make  actual  engineering  applications 
of  classroom  theory  to  aircraft.  This  would  probably  require  more 
funding,  but  the  potential  seems  enormous. 

III.  Conclusions 

The  Aeronautics  Department  has  the  beginnings  of  an  ongoing  airborne 
laboratory  program.  The  Instructional  method  effectively  links  the 
classroom  theory  to  practical  applications  of  flight  mechanics.  As  a 
laboratory  for  an  academic  course  such  as  Aero  495  and  as  a  future 
research  tool,  the  airplane  is  a  significant  addition  to  our  engineering 
curriculum.  While  the  future  development  of  our  current  program  lies  in 
enrolling  more  cadets,  in  stimulating  research  by  cadets  and  faculty,  and 
even  in  acquiring  our  own  airplane;  for  now,  at  least,  we  have  succeeded 
in  solving  an  old  p  'blera  by  leasing  an  aircraft  and  using  our  own  faculty 
members  as  pilots. 

Appendix 

USAF  Academy,  Department  of  Aeronautics 
Aero  495  -  Test  Plan 

Sierra  C24R  Limited  Performance  Evaluation 
June  >982 


Introduction 

A  limited  performance  evaluation  of  the  Beech  Sierra  C24R  will  be 
conducted  at  the  U.S.  Air  Force  Academy  by  Department  of  Aeronautics 
(DFAN)  faculty  pilots  and  the  students  enrolled  in  Aero  495.  Flight 
testing  will  be  conducted  during  the  Fall  1982  semester  from  the  fifth  to 
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tenth  week  of  classes.  Results  of  the  evaluation  will  be  presented  in  a 
formal  oral  report  given  by  each  of  the  two  student  test  teams. 

Objectives 

The  objectives  of  this  evaluation  are  to  determine  the  Sierra 
C24R's  general  performance  characteristics  and  to  compare  them  to  the 
contractor's  Flight  Manual.  In  addition,  certain  contractual  guarantees 
are  verified.  Specific  objectives  follow. 

A .  Takeoff  Performance 

-  determine  takeoff  power  ground  roll  using  the  Flight  Manual 
takeoff  procedure 

B .  Climb  Performance 

-  determine  the  full  throttle  maximum  rate  of  climb 

-  determine  the  full  throttle  best  angle  of  climb 

-  verify  the  climb  performance  predicted  in  the  contractor's 
Flight  Manual 

C .  Level  Turn  Performance 

-  determine  the  level  sustained  turn  performance  in  cruise  power 
at  2,700  rpm 

-  determine  the  speed  for  optimum  sustained  turn  performance  at 
the  test  altitude 

D .  Cruise  Performance 

-  determine  the  airspeeds  and  rpm  for  maximum  range  and  maximum 
endurance  as  derived  from  test  data 

-  determine  the  aircraft  drag  polar 

-  compare  test  results  with  the  contractor's  Flight  Manual 

E .  Descent  Performance 

-  determine  the  propeller  windmilling  best  no  wind  glide  ratio 

-  determine  the  best  glide  speed  and  minimum  sink  speed  with 
propeller  windmilling 
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-  compare  test  results  with  the  contractor's  Flight  Manual 
maximum  glide  configuration  of  91  knots 

F .  Contractual  Guarantees 

-  maximum  speed  at  sea  level  -  142  knots 

-  cruise  speed  at  75  percent  power,  10,000  feet  -  137  knots 

-  range  with  45  minute  reserve  at  75  percent  power  at  10,000 
feet  -  646  nautical  miles 

-  rate  of  climb  at  sea  level  -  927  fpm 

-  service  ceiling  -  15,385  feet 

-  glide  ratio  -1.7  nautical  miles  per  1000  feet 


Authority 

This  test  program  will  be  conducted  by  Department  of  Aeronautics 
faculty  and  students  as  an  integral  part  of  the  curriculum  for  Aero  495, 
course  in  flight  test  techniques.  The  program  has  the  approval  of  the 
Superintendent,  the  Dean  of  the  Faculty,  the  Head  of  the  Department  of 
Aeronautics,  and  the  Director  of  Flight  Operation  of  Hedrick  Beechcraft 


Test  Team  Organization 

Test  team  organization  shown  in  Figure  1  will  consist  of  two  DFAN 
faculty  pilots  and  two  student  flight  test  engineer  teams.  Each  test  team 
will  be  assigned  to  fly  with  one  faculty  pilot.  A  Test  Director  for  each 
team  will  be  appointed  to  coordinate  the  entire  evaluation  effort.  He 
will  in  turn  appoint  individuals  to  be  in  charge  of  each  test  area  (i.e., 
data  monitors).  It  will  be  the  data  monitor's  responsibility  to  specify 
the  tests  to  be  flown  in  support  of  his  test  area.  Test  areas  to  be 
assigned  are  takeoff,  climb,  turn,  cruise  and  descent  performance. 


Scope /Schedule 


The  evaluation  will  consist  of  sorties  as  specified  in  Table  1. 
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Table  1.  Data  Sorties 

Test 

Sorties 

Per  Test  Team 

Flight  Time 
Per  Sortie 

Flight  #1 

Takeoff  Performance 
Cruise  Performance 

Turn  Performance 

2.5 

1  .0 

Flight  #2 

Takeoff  Performance 
Climbs  and  Descents 

2.5 

1  .0 

*Tot  a  1 

5.0 

*One  sortie  will  he  shared  by  hoth  test  teams. 


Flight  #1  and  It 2  are  scheduled  as  shown  on  the  Integrated  Academics  and 
Flying  Schedule  for  Aero  495.  Mission  time  will  not  exceed  1.0  hour. 


Limitations 

The  following  limitations  will  be  observed  during  this  evaluation. 

A.  The  aircraft  will  be  operated  in  accordance  with  the  Airplane 
Flight  Manual,  FAR  Part  91  and  all  Beech  Aero  Club  Operating  Instructions. 

B.  All  data  sorties  will  be  flown  with  one  DFAN  faculty  pilot  and  two 
cadets. 

C.  Testing  will  only  be  accomplished  under  VFR  daytime  conditions  at 
10.000  ft  MSL  and  below. 

D.  All  testing  will  be  accomplished  within  the  local  flying  area  of 
Colorado  Springs 


Test  Aircraft  Description 

The  Beechcraft  Sierra  C24R,  manufactured  by  Beech  Aircraft 
Corporation,  is  a  six-place,  retractable,  general  aviation  aircraft 
powered  by  one  f ue 1 - i n j ec t ed ,  4-cylinder,  200  HP  Avco  Lycoming  engine. 

The  propeller  Is  a  Hartzell  c on st ant  -  speed ,  two-blade,  aluminum-alloy  prop 
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with  spinner.  See  Figure  2  for  general  dimensions  and  Table  2  for 
Aircraft  Limitations. 


Table  2.  Aircraft  Limitations 


Never  Exceed  Speed  (V  ) 

Maximum  Maneuvering  Speed  (V  ) 

Maximum  Cruising  Speed  In  Turbulent  Air  (V  ) 

1G  Stall  Speed  Gear  and  Flaps  Up  (2,600  lbs) 
(power  idle) 

Maximum  Ramp  Weight 
Maximum  Takeoff  Weight 
Maximum  Landing  Weight 

Flight  Maneuvering  Load  Factor  Flaps  Up 
Flight  Maneuvering  Load  Factor  Flaps  Down 
Maneuver,  Bank  Angles  No  More  Than  60  Degrees 

Service  Celling  15,385  feet 

Test  Plan  Celling  10,000  feet 


IAS 

Knot  s/aph 

168/193 

125/144 

143/165 

65/75 


2,785  lbs 
2,750  lbs 
2,750  lbs 

+3.8  to  -1  . 9G 
+  1  .  9G 


Flight  Test  Instrumentation 

All  test  data  will  be  hand  recorded  using  standard  cockpit 
Instrumentation.  The  only  exceptions  are  the  use  of  an  accelerometer,  a 
stopwatch  and  a  cassette  tape  recorder. 


Weight  and  Balance 

Detailed  weight  and  balance  records  for  each  aircraft  are  available 
at  Hedrick  Beechcraft.  Prior  to  every  data  mission,  student  test 
engineers  will  calculate  aircraft  weight  and  balance  data  for  both  takeoff 
and  landing. 

Test  De sc r 1 p 1 1  on/ Pr oc ed ur e s 

Unless  otherwise  noted,  all  performance  tests  will  be  performed  with 
engine  operating,  landing  gear  and  wing  flaps  retracted.  All  data  will  be 
hand  and  voice  recorded,  and  manually  reduced  to  standard  aircraft  weight 
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and  atmospheric  conditions.  Specific  test  techniques  for  each  area  will 
be  covered  in  classroom  lectures  and  handouts  from  references  2,  3,  and  4. 

A .  Takeoff  Performance 

The  takeoff  ground  roll  will  be  determined  and  hand  recorded  for 
each  sortie. 

All  takeoffs  will  be  made  with  15  degrees  of  flaps  in  accordance 
with  the  Sierra  C24R  Airplane  Flight  Manual. 

B .  Climb  Performance 

Climb  data  at  different  airspeeds  will  be  obtained  using  sawtooth 
climbs.  Data  will  be  obtained  at  2,700  rpm  with  mixture  adjusted  for  best 
power  at  test  altitudes  of  8,000,  8,500,  9,000  and  9,500  feet.  Engine 
operating  limitations  as  specified  in  the  Flight  Manual  will  be  followed. 

C  .  Level  Turn  Performance 

Turn  performance  for  the  Sierra  will  be  determined  from  stabilized 
turns  at  various  altitudes  and  airspeeds.  Data  will  be  obtained  at  2,700 
rpm,  or  full  throttle,  between  7,000  and  10,000  feet. 

D .  Cruise  Performance 

Cruise  performance  will  be  evaluated  using  the  P  versus  V 
test  technique  covered  in  references  2  and  4.  Using  the  backside  trim 
shot  technique ,  the  aircraft  will  be  stabilized  at  several  altitude  and 
airspeed  combinations.  Data  will  be  obtained  at  altitudes  between  7,000 
and  10 , 000  feet . 

E .  Descent  Performance 

Descent  performance  will  be  determined  at  various  airspeeds  at 
test  altitudes  of  8,000,  8,500,  9,000  and  9,500  feet.  The  aircraft  will 
be  operated  with  the  throttle  at  idle  and  propeller  at  high  pitch  for 
gathering  descent  data.  The  Flight  Manual  restriction  concerning 
prolonged  idle  settings  will  be  observed. 
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Tra 1 ni ng 

Both  DFAN  faculty  pilots  will  have  at  least  an  FAA  commercial  pilot 
rating  and  be  current  in  the  Beech  Sierra  C24R  in  accordance  with  FAA  and 
Hedrick  Beechcraft  Aero  Club  standards. 

All  cadets  enrolled  in  Aero  495  will  participate  in  the  flying 
portion  of  the  course  as  passengers  only  and  will  receive  appropriate 
aircraft  orientation  andd  safety  instruction.  All  the  performance  flight 
tea  techniques  required  to  gather  test  data  will  be  covered  during 
classroom  lectures  prior  to  the  flights  for  which  they  will  be  used. 

Crew  Duties 

A .  Pilot 

1.  Check  local  flying  weather. 

2.  Brief  students  on  mission  profile,  and  ground  and  in-flight 

safety . 

3.  Check  maintenance  status  of  aircraft  and  perform  pre-flight. 

4.  Provide  a  stopwatch. 

5.  Provide  the  tachometer  reading  at  which  the  aircraft  was 
refueled  and  the  quantity  of  fuel  and  oil  on  board. 

6.  Act  as  pilot  in  command  of  the  aircraft  and  occupy  the  left 
front  seat  at  all  times. 

B  .  Student  s 

1.  Bring  data  cards  and  a  clipboard. 

2.  Complete  aircraft  weight  and  balance  form. 

3.  Complete  takeoff  data  using  temperature  and  pressure  altitude 
provided  by  the  pilot. 

4.  Provide  cassette  tape  player  for  each  flight,  (optional) 

5.  Record  tachometer  reading  at  which  aircraft  was  refueled  and 
the  quantity  of  fuel  and  oil  on  board. 
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6.  Cadets  will  be  assigned  to  two  man  teams  for  purposes  of  taking 
flight  test  data.  Flight  crew  duties  will  be  rotated  each  flight.  Along 
with  the  pilot,  who  will  be  primarily  concerned  with  precisely  flying  the 
aircraft,  both  cadets  will  act  as  lookouts  and  notify  the  pilot 
Immediately  of  any  aircraft  sighted.  The  cadet  in  the  right  front  seat 
will  act  as  data  observer  and  timekeeper  and  the  cadet  in  the  rear  seat 
will  act  as  iata  recorder. 

Safety 

Flight  personnel  will  adhere  to  the  following  while  on  the  fllghtllne 
and  in  and  around  the  aircraft: 

a.  Smoking  is  prohibited  in  or  near  the  aircraft. 

h.  Seat  belts  will  he  worn  at  all  times. 

c.  Flight  personnel  will  be  seated  in  the  aircraft  prior  to  engine 
start  and  will  remain  seated  until  the  engine  is  stopped. 

d.  Remain  clear  of  the  propeller  area. 

e.  Do  not  stand,  walk,  or  lean  on  the  aircraft  except  in  designated 
areas. 

f.  Do  not  open  aircraft  windows  or  doors  in  flight. 

g.  Advise  the  pilot  immediately  upon  observing  another  aircraft. 

b.  Do  not  manipulate  the  aircraft  flight  controls  or  engine  controls 
unless  told  to  do  so  by  the  pilot. 

I.  Advise  the  pilot  of  impending  airsickness.  Use  the  bag  provided, 
your  hat,  y  our  shoe,  anything  except  the  floor  of  the 

aircraft  . 

J.  Stay  clear  of  taxiing  aircraft  and  other  flightline  vehicles. 
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Command  and  Control 

All  testing  to  be  accomplished  will  be  for  academic  purposes  only  and 
will  be  performed  within  the  restrictions  of  the  Flight  Manual,  Part  91  of 
the  Federal  Aviation  Regulations,  Hedrick  Beechcraft  Aero  Club  Rules  and 
the  limitations  imposed  by  this  test  plan. 

All  Information  with  respect  to  this  test  plan  Is  unclassified. 

Test  Plan  Amendments 

An  amendment  to  this  test  plan  is  required  If  the  flight  test 
envelope  is  expanded  or  if  any  limitations  In  the  test  plan  are  made  less 
restrictive.  An  amendment  to  the  test  plan  must  be  reviewed  and  approved 
by  the  same  authority  who  approved  the  basic  plan. 

References 

1 .  Pilot's  Operating  Handbook  and  FAA  Approved  Airplane  Flight  Manual 
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Aeronautical  History 
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HISTORY  OF  THE  HIGH  ALTITUDE  RECONNAISSANCE  AIRCRAFT 

Ben  R.  Rich* 

Editor ' s  Note 

Mr.  Ben  R.  Rich  has  been  involved  in  aircraft  design  since  1950,  when 
he  joined  Lockheed  Corporation.  At  Lockheed  he  worked  with  the  legendary 
C.L.  "Kelly”  Johnson  and  was  involved  in  the  development  of  all  of  the 
Lockheed  high  altitude  reconnaissance  aircraft  including  the  U-2  and 
SR-71.  Since  1975  he  has  been  the  director  of  Lockheed's  "Skunk  Works" 
organization,  which  is  noted  for  developing  advanced  aircraft.  The 
evolution  of  these  aircraft  is  described  in  the  following  paper  presented 
to  the  faculty  and  cadets  of  the  United  States  Air  Force  Academy  on 
November  16,  1982. 

It  is  Indeed  a  great  pleasure  to  address  the  cadets  and  faculty  of 
the  United  States  Air  Force  Academy.  Today  I  would  like  to  discuss  the 
four  kinds  of  aircraft  development  that  I  have  observed  at  Lockheed's 
Advanced  Development  Projects  Division  or,  as  it  is  more  commonly  known, 
the  "Skunk  Works.” 

First,  there  is  what  I  call  technology  push.  This  is  design 
progression  due  to  technological  advances  —  in  aerodynamics,  propulsion, 
materials,  avionics,  etc.  —  resulting  in  increased  capability  in  speed, 
altitude,  maneuverability,  payload,  persistence,  and  endurance.  The 
evolution  of  the  World  War  I  S pad  into  today's  F-15  is  a  prite  example  of 
technology  push;  at  the  Skunk  Works  I've  seen  the  P-80s  develop  into  the 
F-l 04  s . 

Second,  there  is  progressive  growth  or  the  Improvement  of  a 
particular  model  with  time.  Douglas'  A-4 ,  for  example,  advanced  through 
a  number  of  stages  before  the  final  product,  Model  N,  was  achieved.  We 
had  a  similar  experience  at  the  Skunk  Works  with  the  F-94  series,  from  the 
A  model  to  the  D  model. 

Third,  there  is  what  I  call  requirement  pull:  this  occurs  when  an 

aircraft  is  designed  to  meet  a  special  purpose  or  requirement.  Cood 

•Vice  President  and  General  Manager,  Advanced  Development  Projects, 
Lockheed-Cal if  or nla  Company. 
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exampl e  s  are  the  U-2  and  the  SR-71,  which  were  developed  to  meet 
rec onna 1 s sa nc e  needs. 

Fourth,  there  is  the  breakthrough.  Examples  of  this  include  the  Jet 
engine  and,  more  currently,  the  “low  observable"  or  stealth  aircraft. 

At  the  Skunk  Works  all  four  kinds  of  aircraft  development  occur.  I 
could  tell  you  stories  dating  from  the  time  of  the  XP-38  all  the  way  up  to 
the  present,  but  today  I  will  restrict  my  talk  to  the  evolution  of  the 
different  Blackbirds  —  all  developed  since  I  joined  the  Skunk  Works  in 
1  95 A  . 

When  President  Eisenhower's  “Open  Skies  Policy"  was  rejected  by  the 
Russians  in  1954,  Kelly  Johnson,  the  father  of  the  Skunk  Works  and  my 
predecessor,  made  an  unsolicited  proposal  to  the  Department  of  Defense  to 
develop  a  high-altitude  aircraft:  the  U-2A.  This  aircraft  is  shown  in 
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There  was,  however,  a  mismatch  between  the  aerodynamics  and 
propulsive  thrust  of  this  aircraft.  At  altitudes  above  70,000  feet  the 
airplane  had  a  cruise  speed  envelope  of  four  knots.  If  the  airplane  flew 
too  fast,  it  went  into  Mach  buffet;  if  it  flew  too  slowly,  it  went  into 
stall  buffet.  Even  with  this  shortcoming,  however,  the  aircraft  was  very 
successful.  From  1956  through  the  eventful  May  Day  of  1960,  when  a  U-2 
was  eventually  shot  down  by  the  Russians,  we  overflew  Russia  regularly, 
taking  photographs  at  will.  These  U-2s  were  sometimes  followed  by  as 
many  as  35  Russian  fighter  aircraft  flying  20,000  feet  below.  This 
annoyed  our  pilots,  since  this  aluminum  shield  obstructed  their 
photographs. 

In  late  1965  Kelly  Johnson  decided  to  design  a  wing  for  the  U-2  that 
would  match  the  aircraft's  thrust.  We  looked  at  a  variety  of  fancy- 
wings,  such  as  the  Whitcomb  supercritical  wing,  but  since  none  of  them 
could  match  the  lift-to-drag  (L/D)  ratio  of  the  basic  wing  at  the  high 
values  of  coef f lc lent-of-llf t  required,  we  simply  added  40  percent  more 
area  to  the  existing  U-2  wing.  The  resulting  aircraft  was  the  U-2R  ahown 
in  Figure  3. 


Figure  3.  Lockheed  U-2R 
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By  now  the  U-2s  were  doing  all  types  of  real  time  reconnaissance  (i.e., 
the  results  of  the  reconnaissance  were  immediately  available),  such  as 
photo,  Elint,  and  Comint.  We  added  the  slipper  tanks  and  large  wing  pods 
that  are  still  used  today,  until  the  gross  weight  of  the  aircraft  reached 
40,000  pounds.  Figure  4  shows  a  comparison  of  the  U-2C  to  the  U-2R. 


Figure  4.  A  Side-by-Side  Comparison  of  the  U-2C  and  the  U-2R 

U-2  production  ended  in  1968  --  and  the  25,000  tools  were  hidden  in  five 
different  places  so  that  some  idiot  wouldn't  have  them  cut  up  to  sell  as 
scrap.  This  shortsighted  policy  was  supposed  to  save  money  for  the  United 
States  government. 

In  1977  Kelly  and  I  met  with  General  David  Jones  to  discuss  the 
potential  of  the  U-2  as  a  tactical  reconnaissance  vehicle.  As  a  result, 
Ceneral  Jones  set  up  a  joint  Air  Force/Army  tactical  reconnaissance  study 
which  selected  the  updated  U-2  as  the  best  vehicle  for  this  job.  In  1979 
the  tactical  version  of  the  U-2  (renamed  the  TR-1)  went  back  into 
production.  Since  then,  six  aircraft  have  been  delivered,  five  to  Beale 
Air  Force  Base  and  one  (designated  an  ER-2)  to  NASA  Ames.  Figure  5  is  a 
photograph  of  the  TR-1 . 
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Figure  5.  Lockheed  TR-1 


In  1975  che  Skunk  Works  and  Sperry/Salt  Lake  made  a  proposal  Co  Che 
Air  Force  co  build  a  U-2  RPA .  An  RPA  (remotely  piloted  airplane)  differs 
from  an  RPV  (remotely  piloted  vehicle)  in  that  the  vehicle  can  be  either 
manned  or  unmanned.  The  control  could  be  developed  with  man-ln-loop 
during  early  flights;  during  later  flights  the  pilot  could  operate  the 
vehicle  from  the  ground.  This  system,  called  the  Sperry  Flight  System,  is 
now  used  on  the  F-102  drones.  The  RPA  could  be  highly  effective  In  these 
days  of  large,  expensive  drones  with  very  expensive  black  boxes; 
unfortunately,  we  never  sold  the  concept. 

In  1958,  when  Russian  aircraft  began  to  fire  air-to-air  rockets  at 
the  U-2 s ,  we  began  the  design  of  the  first  of  four  supersonic  Blackbirds, 
the  A-l 1  *  shown  in  Figure  6.  This  was  a  Mach  3+,  ov er-80 , 000-f oo t- 
altitude,  single-place,  reconnaissance  (recce)  aircraft.  Its  first  flight 
occurred  on  April  26,  1962  .  The  A  —  1 1  was  followed  in  1964  by  the 
interceptor  version  shown  in  Figure  7  and  in  1965  by  the  two-place  current 
reconnaissance  aircraft,  the  SR-71,  shown  in  Figure  8.  The  latter  was 
originally  named  the  RS-71  ,  but  when  President  Johnson  goofed  and  called 
it  the  SR-71  we  had  to  change  all  our  drawings. 


•The  term  A-12,  which  is  the  correct  airplane  number,  is  also  used.  The 
aircraft  became  known  as  the  A-l 1  because  of  an  Air  Force  Press  release. 


Figure  6.  Lockheed  A-ll 


Figure  7.  Interceptor  Version  of  Lockheed  A-ll 


Figure  8.  Lockheed  SR-71  (Originally  Named  the  RS-71) 
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These  aircraft  became  operational  In  1969,  and  we  stationed  them  at  Beale 
Air  Force  Base.  Seventeen  of  the  D-21s  are  currently  stored  at 
Da v 1 s- Mont han  Air  Force  Base. 

Not  everything  that  we  do  at  the  Skunk  Works  Is  successful.  Figure 
11  shows  that  we  have  had  some  failures  too.  But  remember,  Jesus  Christ 
was  only  91.5  percent  perfect:  he  missed  on  one  of  his  twelve  apostles. 


SKUNK  WORKS  PROJECTS 


DEVELOPED  IV 
KIUY  JOHNSON 
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X 
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X 

F-94A 
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X 

F-94C 

19.000 
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X 

X 

LOW  CDIT.  MACH  K 

XFV-I 
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X 

BAD  POWER  PLANT 

SATURN 
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X 
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X 
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|  XF-90 

22.500 
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X-7 
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X 
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X 
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T2V 
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X 

XF-104 
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EXCEL 

X 
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X 

U-2 
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X 

JiTSTAR 

17.500 
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X 

CL-400 
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X 

POOR  CONCEPT 

MODEL  A 

SECRET 
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X 

I  YF-I2A 

SECRET 
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X 

8R-7I 

SECRET 
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X 

MODEL  D 

SECRET 

CfCKL 

X 

MODEL  R 

SECRET 

EXCEL. 

X 

•  PARTIAL  SKUNK  WORKS  METHOD 


Figure  11.  A  List  of  Projects  at  the  Skunk  Works 


Let  me  expand  on  one  of  our  boo-boos.  We  called  this  bird,  shown  in 
Figure  12,  the  CL-400.  It  was  a  Mach  2.5,  1 00 , 000-f oo t- a  1 1 1 1 ud e  , 
hydrogen-powered  reconnaissance  aircraft. 
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Figure  12.  The  Proposed  Lockheed  CL-400 


When  the  Russians  brought  their  liquid  hydrogen  (LH2)  expert  Peter 
Kapitsa  out  of  prison  camp  In  1955,  we  thought  they  were  building  a 
hydrogen  airplane  to  shoot  down  the  U-2 .  So  the  following  year  we 
received  a  contract  to  build  16  LH2  recce  aircraft.  It  was  to  be  powered 
by  a  new  type  of  engine:  a  hydrogen  expansion  engine  designed  by  Pratt  4 
Whitney.  Both  Pratt  &  Whitney  and  the  Skunk  Works  learned  how  to  handle 
tons  of  liquid  hydrogen  —  all  before  the  great  space  age.  The  CL-400 
carried  a  quarter  of  an  acre-foot  of  LH2.  Its  stable  or  equilibrium 
temperature  as  a  liquid  Is  20  degrees  C  above  absolute  zero  temperature 
-253  degrees  C  (-420  degrees  F).  One  version  of  the  airplane  was  about 
two-thirds  the  length  of  s  football  field  with  the  fuselage  diameter  of  a 
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widebody  jet  (20  feet)  and  the  vlnga  of  an  F-104.  (In  thoae  days  we 
didn't  know  too  much  about  blended  wing  bodlea.)  The  engine' a  TSFC 
(thruat  specific  fuel  conaumption,  the  engine' a  aaaa  flow  rate  of  fuel 
divided  by  Its  thrust)  was  supposed  to  be  0.6  at  Mach  2.5,  100,000  feet; 

It  turned  out  to  be  0.8.  So  between  the  low  L/D  and  high  TSFC  we  had  an 
airplane  with  a  1900  nautical  mile  range  with  no  possible  future 
lapr oveaent . 

Even  though  Pratt  &  Whitney  had  the  engine  running  well,  and  we  had 
produced  15  allea  of  aetal  extrusion,  we  were  going  to  turn  out  a 
"vide-body "  dog  (failure).  What  would  happen  If  the  airplane  landed  at 
Mt .  Hoae  Air  Force  Base  or  Pease  Air  Force  Base?  Where  would  the  Air 
Force  get  liquid  hydrogen?  If  it  went  operational  overseas,  where  could 
we  find  LH2  In  Turkey,  for  exaaple?  So  Kelly  Johnson  and  Perry  Pratt 
spoke  to  either  General  Irvine  or  General  LeMay  and  recoaaended  that  we 
cancel  the  prograai.  As  It  turned  out,  this  was  a  wise  decision.  If  we 
had  built  the  16  airplanes,  we  would  have  used  one-quarter  of  the  natural 
gas  cooing  into  Los  Angeles  —  right  in  the  middle  of  the  1974  energy 
crisis.  By  the  way,  we  found  out  why  the  Russians  brought  Kapitsa  out  of 
prison  when  they  launched  the  Sputnik  with  a  LH2-fueled  rocket  In  1957. 

Unfortunately,  I  cannot  tell  you  what  we  have  been  doing  for  the  past 
10  years.  We  seem  to  score  a  breakthrough  at  the  Skunk  Works  every 
decade,  so  If  you  Invite  me  back  In  10  years  I  will  be  able  to  tell  you 
what  we  are  doing  now.  I  can  tell  you  about  a  contract  that  we  recently 
received  for  the  vehicle  shown  in  Figure  13.  The  Skunk  Works  has  been 
assigned  the  task  of  getting  "E.T."**  back  home. 

**Extra-Terre str lal  (from  the  movie  of  the  same  name) 
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SECTION  V 

The  Engineer's  Bookshelf 
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THE  ENGINEER'S  BOOKSHELF: 

CLASSICS  OF  CONTEMPORARY  SCIENCE  WRITING  AND  THE  NUCLEAR  DEBATE: 

SOME  SUGGESTED  READING 

James  M.  Kempf* 

In  the  last  issue  of  the  Aeronautic:  Digest  I  asserted  that  the 

best  writing  on  technical  and  scientific  subjects  often  concerns  the 

Impact  of  technology  on  human  society*  The  books  I  am  recommending  in 
this  essay  do  too.  The  writers  I  call  "classics"  of  contemporary 
scientific  literature  all  demonstrate  wide-ranging  knowledge  and  a 
stubborn  allegiance  to  a  singular  belief:  that  no  matter  how  technical  and 
specialized  our  contemporary  culture  has  become,  true  science  confirms  an 
age-old  faith  in  the  interdependence  of  knowledge  as  discovered  by  the 
world's  only  creative  and  thinking  species,  man. 

As  many  writers  have  argued  in  modern  times,  a  major  deficiency  of 

contemporary  education  is  its  increasingly  narrow  training  in  both 

intellectual  disciplines  and  professional  vocations.  The  sad  result  of 
this  narrowed  training  is  the  arrogant  and  ignorant  temperament  It 
produces  in  the  technocratic,  bureaucratic  personalities  that  thrive  in 
modern  institutions.  The  fear  of  complexity,  of  dissent,  of  inquiring 
minds,  and  of  the  questioning  of  values  often  paralyzes  contemporary 
institutions.  This  timid  habit  of  mind  leads  to  the  stifling  of  free 
intellect,  or  more  commonly  to  the  division  of  knowledge  into 
narrowly-focused  concerns,  bo  that  we  lose  our  sense  of  the  whole  and  the 
purpose  of  education  Itself,  which  should  be  the  search  not  just  for  facts 
but  for  wisdom.  By  separating  the  quest  for  knowledge  about  things  from  a 
quest  for  the  values  we  require  to  shape  our  use  of  knowledge,  modern 
education  too  often  is  satisfied  with  routine  mental  exercise,  not  the 
more  difficult  task  of  building  a  just  civilization. 

*Asai8tant  Professor  of  English,  USAF  Academy 
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But  the  concerns  of  the  writers  I  wish  to  bring  to  your  attention  in 
this  essay  and  the  most  widely  discussed  public  issue  of  this  past  year 
concerning  technology  --  the  issue  of  nuclear  weapons  —  Illustrate  the 
impossibility  of  separating  discussion  of  technology  from  the  broad  issues 
of  social  debate.  All  these  writers  argue,  often  from  quite  different 
perspectives,  that  if  we  are  to  accurately  understand  our  technically 
complex  civilization,  we  must  understand  how  science,  philosophy,  art, 
politics,  history,  economics,  moral  values,  and  methodology  as  well  as 
psychology  interact  and  affect  our  perception  of  "reality"  as  well  as  of 
truth.  I  begin  hy  recommending  the  works  of  several  writers  who  are,  or 
were,  distinguished  scientists,  and  who  are  also  distinguished  writers. 

Few  scientists  illustrate  better  than  Lewis  Thomas  how  an  expansive 
mind  is  a  necessary  complement  of  a  technical  researcher  whose  focus,  on 
the  surface,  is  the  most  minute  subject  of  natural  science.  Thomas  is  a 
microbiologist  who  has  spent  his  life  studying  biological  organisms  and 
systems  that  the  human  eye  can't  even  discern.  Vet  his  essays,  in  two 
lucid,  beautifully  written  volumes,  The  Lives  of  a  Cell;  Notes  of  a 
Biology  Watcher  (Viking,  1974)  and  The  Medusa  and  the  Snail:  More  Notes 
of  a  Biology  Watcher  (Viking,  1979),  demonstrate  a  mind  that  ranges  all 
across  the  phenomena  of  human  society.  Thomas  writes  about  everything,  it 
seems,  from  the  ways  in  which  our  advertising  industry  can  misinform  the 
public  about  the  dangers  of  bacterial  organisms  to  how  the  fear  of  genetic 
manipulation  (as  well  as  Ignorance  of  biology)  leads  to  unwarranted 
political  passions  that  fabricate  unnecessary  restrictions  on  recombinant 
DNA  research.  Thomas  is  at  ease  discussing  how  disease,  death,  computers, 
and  public  health  practices  are  fruitfully  analyzed  not  only  by  the 
techniques  of  scientists  but  also  by  the  imaginative  perceptions  of  the 
Argentinian  novelist  Jorge  Luis  Borges.  He  is,  in  short,  a  science  writer 
whose  work  "educates  by  delight,"  an  old  adage  about  good  art  as  well  as 
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good  science.  It  Js  perhaps  indicative  of  his  skill  that  The  Lives  of  a 
Cell  was  honored  by  the  publishing  Industry  in  1975  with  a  National  Book 
Award  —  for  belles-lettres,  not  science. 

Another  equally  Instructive  and  learned  scientist  was  the  late 
paleontologist  Loren  Elseley.  In  a  number  of  books  about  human  evolution 
and  history,  Elseley  brought  to  his  studies  a  poetic  prose  style  and  a 
grand  historical  perspective.  His  essays  range  from  discussion  of  ancient 
paleontology  to  meditations  on  the  "naturalist"  explorations  of  American 
writers  such  as  Henry  David  Thoreau.  The  best  introduction  to  Elseley's 
thought  is  a  volume  of  selected  essays  and  poetry  taken  from  a  lifetime  of 
his  work.  Titled  The  Star  Thrower  (Harcourt  Brace  Jovanovich,  Harvest 
Books,  1978),  this  volume  reveals  Elseley  to  be  a  profound  thinker  and 
scientist,  one  whose  writing  demonstrates  lucid  intelligence  and 
familiarity  with  subjects  as  diverse  as  aesthetic  philosophy  and  the 
"artistic"  skills  which  Cro-Magnon  man  demonstrated  thousands  of  years 
ago.  Elseley's  essay  "The  Illusion  of  Two  Cultures"  might  serve  as  a 
representative  example  of  his  work.  In  this  essay  he  angrily  denounces 
the  "trade  union"  habits  of  much  modern  science  and  meditates  on  his 
discovery  that  the  most  practical  tools  of  ancient  man  were  decorated  or 
designed  as  works  of  art.  He  also  shows  how  the  symbolic  constructs  of 
human  imagination  have  often  become  the  basis  of  later  scientific  theory. 
Elseley  may  well  have  been  an  anachronism,  since  his  commanding  scientific 
concern  for  the  multiple  dimensions  of  human  society  and  his  sensitive 
appreciation  for  the  wonder  of  nature  made  him  turn  to  poetry  as  well  as 
science  to  express  himself.  His  achievement  was  unique,  for  among  his 
many  professional  Involvements  and  awards  was  his  election  to  the  National 
Institute  of  Arts  and  Letters,  an  almost  unheard  of  honor  for  a  scientist. 

Fred  Alan  Wolf  is  a  physicist  and  educator  who  won  the  American  Book 
Award  in  1982  for  the  outstanding  paperback  book  published  the  previous 
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year  In  science.  His  book  describing  the  revolutionary  concepts  of 
twentieth-century  oarticle  physics  and  quantua  mechanics.  Taking  the 
Quantum  Leap,  (Harper  and  Row,  1981),  la  a  pleasurable  excursion  through 
the  history  of  modern  physics.  The  effect  of  that  revolution  on  our 
understanding  of  natural  philosophy  and  the  structure  of  the  universe  Is  a 
collateral  theme  of  Wolf's  book.  By  discussing  concepts  such  as 
Heisenberg's  uncertainty  principle  and  Niels  Bohr's  and  Max  Planck's 
analyses  of  wave  mechanics  and  the  quantum  nature  of  light.  Wolf  connects 
this  history  of  physics  to  the  concurrent  modern  revolutions  in 
psychology,  epistemology,  and  art,  which  all  exploded  during  the  early 
part  of  this  century.  The  virtue  of  Wolf's  book  Is  his  ability  to 
discuss  wllth  ease  and  Insight  the  connection  between  such  apparently 
disparate  subjects  as  wave  mechanics,  perceptual  psychology.  Irrational 
thought  processes,  and  Dadaist  art.  Taking  the  Quantum  Leap  la  a  book 
well  worth  reading. 

Men  like  Thomas  and  Elseley  gently  but  firmly  have  warned  in  several 
essays  that  science  Itself  Is  too  often  endangered  by  rigidity  of  mind, 
hostility  to  new  theory  and  evidence,  and  abuse  of  the  prestige  which 
modern  society  has  bestowed  on  it.  No  better  example  exists  of  these 
dangerous  tendencies  than  Stephen  Jay  Gould's  book,  The  Mlsmeasure  of 
Ma n  (Norton,  1981),  which  won  the  National  Book  Critics  Circle  Award  for 
nonfiction  In  1981.  A  science  writer  of  Increasing  prominence,  Gould 
traces  the  scandalous  ways  In  which.  In  the  name  of  science,  such  quackery 
as  phrenology,  rigged  intelligence  testing,  and  downright  falsified  data 
have  been  used  by  so-called  scientists  over  the  past  century. 

Intelligence  tests,  Gould  demonstrates,  are  often  really  Inspired  by 
political  and  economic  motives  to  suppress  racial  minorities  and 
economically  underprivileged  classes  of  people.  Bis  book  begins  with  an 
epigraph  from  Darwin's  Voyage  of  the  Beagle; 
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tf  the  misery  of  our  poor  be  caused  not  by  the  laws  of 
nature,  but  by  our  Institutions,  great  Is  our  sin. 

It  Is  a  curious  comment,  one  calling  for  moral  judgment,  by  a  scientist 
often  erroneously  cited  by  "social  darwlnlsts”  as  support  for  social  and 
economic  exploitation.  The  ways  in  which  the  "sins  of  science”  have 
fostered  social  oppression  is  the  subject  of  Gould's  marvelous, 
stimulating  history  of  intelligence  testing  and  scientific  fraud. 

The  fact  that  modern  science  has  been  viewed  as  a  subversive  threat 
to  religion  and  has  been  antagonistic  to  religious  and  mythological  world 
views  for  four  centuries  is  a  common  theme  of  Intellectual  history.  In  a 
recent  book  by  one  of  our  finest  living  9c  lent i s t- wr i t e r s ,  Robert  Jastrow, 
the  director  of  NASA's  Goddard  Space  Flight  Center,  this  history  is  shown 
to  have  come  full  circle.  Jastrow's  little  book  God  and  the 
Astronomers  (Norton,  1978)  briefly  describes  the  history  and  biographies 
of  the  men  who  revolutionised  modern  astronomy  and  developed  the  so-called 
"Theory  of  the  Expanding  Universe."  In  short  discussions  of  men  like 
Vesta  Melvin  Slipher,  the  astronomer  who  in  1913  discovered  the 
"red-shift"  of  the  retreating  motion  of  galaxies,  and  Edwin  Hubble,  who 
first  used  Sllpher's  data  and  Einstein's  relativity  theory  to  postulate 
the  "Law  of  the  Expanding  Universe,”  Jastrow  concludes  with  a  startling 
theme.  The  Implications  of  modern  astronomy  have  led  scientists  to 
conclude  that  the  universe  began  with  a  "big  bang"  and  is  expanding  on  a 
scale  that  boggles  the  mind.  The  result  of  this  new  cosmology  is  to  show 
that.  Instead  of  conflicting  with  religious  doctrine,  the  theories  of 
modern  scientific  astronomy  may  in  fact  be  congruent  with  the  religious 
views  of  the  origin  of  the  universe.  Two  other  popular  books  by  Jastrow 
that  show  how  interdependent  the  disciplines  of  both  science  and 
humanistic  study  have  become  are  Until  the  Sun  Dies  (Norton,  1977)  and 
Red  Giants  and  White  Dwarfs  (Warner  Books,  1967/1979). 
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Perhaps  no  recent  public  debate  demonstrates  so  clearly  the 
i apo ss i b i 1 i t y  of  trying  to  Isolate  science  and  technology  in  an  artificial 
vacuum  freed  from  social  implications  as  the  debate  the  past  year  over 
nuclear  weapons  policy  and  technology.  That  discussions  of  nuclear  policy 
are  not  constrained  only  by  technical  questions  is  evident  in  the  Great 
Nuclear  Debate  of  1982,  which  touched  on  a  variety  of  related  issues. 
Including  religious  values  and  morality,  foreign  and  military  policy  and 
strategy,  economic  and  political  values,  and  the  psychological  as  well  as 
biological  and  physical  damage  to  the  western  alliance  caused  by  present 
military  strategy  and  tactics.  For  anyone  who  wishes  to  know  more  about 
this  debate  I  have  listed  below  a  number  of  articles  and  books  that  should 
be  informative. 

Several  works  in  particular  are  useful  starting  points  if  one  wishes 
to  understand  the  context  of  the  current  debate.  Two  English  authors, 

E.P.  Thompson  and  Mary  Kaldor,  have  written  influential  essays  and  books 
that  explain  the  position  and  attitude  of  the  leaders  of  Europe's  nuclear 
disarmament  movement.  Kaldor' s  essay  "Nuclear  Weapons  and  the  Atlantic 
Alliance"  (in  the  January  1982  issue  of  democracy  magazine)  outlines  not 
only  the  European  movement's  desire  to  achieve  a  nuclear  weapons-free 
continent,  but  explains  the  political  thrust  of  the  movement,  which 
Includes  "an  explicit  rejection  of  current  relations  of  power."  Tracing 
the  nuclear  policies  of  the  United  States  and  NATO  to  Cold  War  history, 
Kaldor  argues  that  European  political  and  economic  Interests  have  diverged 
from  the  United  States,  and  the  problem  for  Europe  is  to  free  Itself  from 
its  hostage  relationship  to  both  the  United  States  and  the  Soviets.  In 
another  article,  "Beyond  'No  First  Use':  What  the  Peace  Movement  Really 
Means,"  (The  Nation ,  June  26,  1982),  Kaldor  reaffirms  the  disarmament 
aovesci.t'  s  position  that  its  real  demand  is  for  "changing  political 
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relations'  so  that  Europe  does  not  become  the  theatre  for  limited  nuclear 

war  . 

Thompson’s  ideas  are  summarized  In  his  article  "letter  to  America," 
The  Nation  (January  24,  1981),  and  in  two  books:  Beyond  the  Cold  War;  A 
New  Approach  to  the  Arms  Race  and  Nuclear  Annihilation  (Pantheon  1981), 
and  Protest  and  Survive,  co-edited  with  Dan  Smith  (Monthly  Review  Press, 
1981).  A  good  analysis  of  the  views  of  these  leaders  of  Europe’s  nuclear 
disarmament  movement  Is  Alan  Wolfe’s  essay  "Europe  in  Search  of 
Autonomy,"  The  Nation  (February  27,  1982). 

The  reasons  behind  America's  nuclear  freeze  movement  are  most 
clearly  articulated  by  Randall  Forsberg  in  her  essay  "A  Bilateral 
Nuc 1 ea r-Wea pon  Freeze,"  Scientific  American  (November,  1982).  Forsberg 
wrote  the  April  1980  document,  “Call  to  Halt  the  Nuclear  Arms  Race,"  which 
led  to  the  June  1982  anti-nuclear  rally  in  New  York’s  Central  Park  and, 
eventually,  to  the  numerous  freeze  resolutions  on  the  ballot  in  many 
states  during  the  November  1982  elections.  The  article  cited  above 
outlines  the  rationale  of  the  freeze  campaign,  particularly  its  attempt  to 
prevent  production  of  new  " c oun t e r f or c e  weapons”  and  to  freeze  current 
force  levels  at  a  relative  position  of  "parity."  The  article  is 
accompanied  by  excellent  charts  of  the  comparative  levels  of  weapon 
systems  and  warheads  possessed  by  the  United  States  and  the  U.S.S.R. 

An  equally  influential  work  of  the  American  anti-nuclear  movement  is 
Jonathan  Schell's  book  The  Fate  of  the  Earth  (Avon,  1982).  Schell  has 
been  widely  criticized,  as  well  as  praised,  for  his  idealistic  view  that 
we  must  undergo  a  change  of  consciousness  to  save  the  earth  from  what  he 
believes  Is  an  Increasing  risk  of  nuclear  holocaust.  Though  strongly 
criticised  for  his  unrealistic  Idealism,  the  virtue  of  his  book  Is  Its 
delineation  of  the  known  facts  about  nuclear  detonations.  The  dramatic 
scene  he  faints  Is  culled  from  unclassified  government  reports,  studies  of 
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the  National  Academy  of  Sciences,  and  research  data  compiled  by 
biologists,  physicists,  and  medical  personnel  from  experiments  conducted 
to  test  nuclear  radiation,  blast,  and  heat  phenomena,  Hl6  book  has  flaws, 
but  it  has  succeeded  in  forcing  the  general  public  to  think  conceptually 
about  the  real  consequences  of  nuclear  warfare.  That  alone  is  no  small 
achievement . 

Finally,  a  response  to  the  anti-nuclear  movement  written  in  the 
spring  of  1982  by  four  members  of  America's  defense  and  foreign  policy 
establishment  generated  wide  discussion.  The  essay  "Nuclear  Weapons  and 
the  Atlantic  Alliance”  by  McGeorge  Bundy,  George  F.  Kennan,  Robert  S. 
McNamara,  and  Gerard  Smith  in  Foreign  Affairs  (Spring  1982)  describes 
these  men's  proposal  for  changing  America's  military  strategy  to  a  "no 
first  use"  doctrine.  They  view  their  tactic  as  a  way  of  both  allaying 
growing  public  unrest  over  nuclear  policy  and  reuniting,  under  a  common 
policy,  the  fragmenting  forces  within  NATO.  A  no  first  use  policy  would 
"neutralize  the  highly  disruptive"  debate  and  fear  in  Europe  that  the 
United  States  wants  to  fight  a  limited  nuclear  war  in  Europe. 

In  response  to  thesi  various  positions  which  I  have  outlined,  scores 
of  articles  and  hooks  have  been  pouring  forth  from  the  presses  to  weigh  in 
on  the  debate  over  western  military  and  nuclear  weapons  policy.  1  have 
merely  tried  to  list  below  a  number  of  these  essays  and  books.  They 
represent  the  opinions  of  a  large  number  of  influential  policy  makers, 
scholars,  and  military  strategists.  They  also  represent  diverse  political 
viewpoints  and,  consequently,  quite  divergent  responses  to  the  issue  of 
nuclear  strategy.  Many  of  these  authors  respond  to  the  ideas  of  the 
authors  I  have  briefly  summarized  in  the  preceding  paragraphs. 

Leon  Wieseltier,  "Nuclear  War,  Nuclear  Peace:  The  Great  Nuclear  Debate," 
The  New  Republic  (January  10  and  17,  1983);  Albert  Gore,  Jr.,  “The  Fork 
in  the  Road:  A  New  Plan  for  Nuclear  Peace,”  The  New  Republic  (May  5, 
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1982);  Charles  Krauthammer,  "The  Real  Way  to  Prevent  Nuclear  War,”  The 
New  Republic  (April  28,  1982); 

William  J.  Perry  and  Cynthia  A.  Roberts,  "Winning  Through  Sophistication: 
How  to  Meet  the  Soviet  Military  Threat,"  Technology  Review  (July  1982); 
Lynn  R.  Sykes  and  Jack  F.  Everden,  "The  Verification  of  a  Comprehensive 
Nuclear  Test  Ban,"  Scientific  American  (October  1982); 

Edward  N.  Luttwak,  "How  to  Think  about  Nuclear  War,"  Commentary  (August 
1982);  and  Edward  Jay  Epstein,  "Disinformation:  Or  Why  the  CIA  Cannot 
Verify  an  Arms-Control  Agreement,"  Commentary  (July  1982); 

George  Kennan,  "Cease  This  Madness:  Our  Collision  Course  with  the 
U.S.S.R.,”  The  Atlantic  (January  1981)  and  "America's  Unstable  Soviet 
Policy,"  The  Atlantic  (November  1982),  and  The  Nuclear  Delusion 
(1982); 

"The  Nuclear  Debate,"  Part  I  by  Fred  Kaplan,  "Russian  and  American 
Intentions,"  Part  II  by  Jerome  B.  Wlesner,  "Russian  and  American 
Capabilities,"  The  Atlantic  (July  1982);  and  Peter  H.  Stone,  "The  Bomb: 

The  Last  Epidemic,"  The  Atlantic  (February  1982). 

Two  Influential  foreign  policy  journals  have  published  a  continuing 

series  of  essays  on  the  nuclear-military  policy  debate.  Foreign  Affairs 

has  published  the  following  essays  during  the  past  year: 

"Nuclear  Weapons  in  the  1980s”  (Winter  1981/82)  included  essays  by 
Spurgeon  M.  Ke*>ny,  Jr.  and  Wolfgang  K.H.  Panofsky,  "Mad  vs.  Nuts:  The 
Mutual  Hostage  Relationship  of  the  Superpowers";  Christopher  Bertram,  "The 
Implications  of  Theatre  Nuclear  Weapons  in  Europe";  and  Stanley  Hoffman, 
"NATO  and  Nuclear  Weapons." 

"Military  and  Political  Policy"  (Summer  1982)  included  essays  by  Robert  W. 
Komer,  "Maritime  Strategy  vs  Coalition  Defense";  General  Bernard  Rogers, 
"The  Atlantic  Alliance";  and  Karl  Kaiser  et  a  1  .  ,  "Nuclear  Weapons  and 
the  Preservation  of  Peace:  A  German  Response  to  No  First  Use." 

The  Fall  1982  issue  of  Foreign  Affairs  included  the  essays,  "Toward  an 
Overall  Western  Strategy  for  Peace,  Freedom  and  Progress,"  by 
Hans-Dietr Ich  Genscher;  Stansfleld  Turner  and  George  Thibault.  "Preparing 
for  the  Unexpected:  The  Need  for  a  New  Military  Strategy";  and  Francois  de 
Rose,  "Inflexible  Response." 

A  similar  series  of  essays  has  been  run  by  Foreign  Policy  magazine 
and  Includes  the  following  essays: 

"Future  of  Arms  Control"  (Fall  1981)  includes  essays  by  Leon  V.  Slgal, 
"Kennan1 a  Cuts";  Peter  D.  Zimmerman,  "Quota  Testing";  and  David  A. 
Andelman,  "Space  Wars." 

"The  Great  Strategic  Debate"  (Winter  1981-82)  Includes  essays  by  John 
Steinbrunner ,  "Nuclear  Decapitation”;  Robert  J.  Elnhorn,  "Treaty 
Compliance";  and  Lawrence  Freedman,  "NATO  Myths." 
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"The  Nuclear  Jumble"  (Fall  1982)  includes  essays  by  Jonathan  Dean,  "Beyond 
First  Use";  Leon  V.  Slgal,  "Warming  to  the  Freeze";  Jan  M.  Lodal, 
"Finishing  Start";  and  L.  Bruce  van  Voorst,  "The  Critical  Masses." 


Two  military  journals  have  run  continuing  analyses  of  the  nuclear 
debate.  The  Naval  War  College  Review  has  recently  published  the 
following  articles: 

Peter  D.  Zimmerman  and  G.  Allen  Greb,  "The  Bottom  Rung  of  the  Ladder: 
Battlefield  Nuclear  Weapons  in  Europe"  ( November-December  1981); 

Donald  M.  Snow,  "Strategic  Uncertainty  and  Nuclear  Deterrence," 
(November-December  1981); 

Lt .-Commander  T.  Wood  Parker,  "Theater  Nuclear  Warfare  and  the  U.S.  Navy" 
( Janua r y-Febr ua r y  1982). 

Several  essays  in  the  Naval  War  College  Review  issue  of 
September-Oc tober  1982  are  especially  useful: 

Laurence  Martin,  "National  Security  in  an  Insecure  Age”;  Michael 
Mandelbaum,  "The  Future  of  Nuclear  Weapons";  Robert  Komer,  "Is 
Conventional  Defense  of  Europe  Feasible." 

The  book  considered  the  "classic"  study  of  military  strategy  in  the 

age  of  nuclear  technology  is  Bernard  Brodie's  Strategy  in  the  Missile 

Age  (Princeton  University  Press,  1959).  A  good  discussion  of  this 

book's  relevance  to  current  debate  and  American  military  strategy  is 

included  in  Major  Leslie  J.  Hamblin's  "Deterrence:  After  the  Golden  Age,” 

Air  University  Review  ( J a nua r y-Febr ua r y  1982). 

The  Air  University  Review  also  recently  published  these  useful 

e  ssays : 

Stephen  M.  Millett,  "Soviet  Perceptions  of  Nuclear  Strategy  and 
Implications  for  U.S.  Deterrence"  (March-April  1982); 

Group  Captain  R.  A.  Mason,  "Western  Deterrence:  Posture  and  Rationale"; 
John  Borowski,  "Theater  Nuclear  Arms  Control  and  Forward-Based  Systems" 
(May-June  1982). 

There  are  scores  of  other  articles  and  books  on  this  issue,  and  the 
bibliography  expands  each  week.  But  the  list  above  should  be  sufficient 
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to  get  anyone  who  is  interested  started  on 
issue  of  technology  and  policy  confronting 
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